Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,435,428 C→A 5.7% L64L (CTG→CTT uspF ← stress‑induced protein, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,4280CA 5.7% 229.1 / 3.6 70L64L (CTG→CTTuspFstress‑induced protein, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base A (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

AATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTT  >  NC_000913/1435281‑1435567
                                                                                                                                                   |                                                                                                                                           
aaTGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGGTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCa                                                                                                                                            <  4:237050/149‑1 (MQ=255)
 atgatgATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAg                                                                                                                                           <  6:227907/149‑1 (MQ=255)
      gatCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCAt                                                                                                                                      >  2:173326/1‑149 (MQ=255)
                   gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGGGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTaa                                                                                                                         <  1:225658/149‑1 (MQ=255)
                      aTCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTc                                                                                                                      >  7:9310/1‑149 (MQ=255)
                        cttcttCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTg                                                                                                                    >  3:44555/1‑149 (MQ=255)
                           cttcGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGc                                                                                                                   <  1:6046/147‑1 (MQ=255)
                              cGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGa                                                                                                                <  8:238531/147‑1 (MQ=255)
                              cGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGGCATCCATTGCCGGTAATTCTGCgg                                                                                                                 <  5:180888/146‑1 (MQ=255)
                              cGCCAATTCCAGAATGCGGTATTTGGGCGAGCCTTCCTCAACATGCACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTTATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCgg                                                                                                                 <  5:160764/146‑1 (MQ=255)
                               gCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATa                                                                                                             >  1:184736/1‑149 (MQ=255)
                               gCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCACTTGCGATTTGGCATCCGCTATCAGGTCATCCATTGCCGGTAATTCTCCGGAAt                                                                                                              <  1:33469/148‑1 (MQ=255)
                                ccAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATAc                                                                                                            >  8:176476/1‑149 (MQ=255)
                                  aaTTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGc                                                                                                          >  8:78468/1‑149 (MQ=255)
                                       cAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCa                                                                                                                                   <  6:53436/119‑1 (MQ=255)
                                       cAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCa                                                                                                                                   >  5:53436/1‑119 (MQ=255)
                                          aaTGCGGTCTTTGGGCGCGCCTTCCTCAACATGGACATGCACTCAGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAAcc                                                                                                    >  6:89093/1‑147 (MQ=255)
                                            tGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCaga                                                                                                <  8:102797/149‑1 (MQ=255)
                                            tGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCa                                                                                                  >  7:133346/1‑147 (MQ=255)
                                              cGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCagaga                                                                                              <  6:344153/149‑1 (MQ=255)
                                               ggTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCagagag                                                                                             >  6:291838/1‑149 (MQ=255)
                                                       ggCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAtagta                                                                                     <  5:115156/149‑1 (MQ=255)
                                                          gAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAggg                                                                                  >  7:104557/1‑149 (MQ=255)
                                                          gAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAggg                                                                                  >  1:187022/1‑149 (MQ=255)
                                                            gCCTTCCTCAACATGGACATGGACTCTGTCGGGTGGCATTTTAAATTTTTTAATGATCTCTTCCAGTTGCGTTTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCa                                                                                <  8:9310/149‑1 (MQ=255)
                                                            gCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCa                                                                                <  2:187022/149‑1 (MQ=255)
                                                               ttCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCTGAATACGCTAAACACAGAGAGGCATAGTAGGGCAGt                                                                              >  7:264918/1‑148 (MQ=255)
                                                               ttCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGt                                                                              >  1:25612/1‑148 (MQ=255)
                                                               ttCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTAAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGt                                                                              >  2:182706/1‑148 (MQ=255)
                                                               ttCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTAAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGt                                                                              >  2:188309/1‑148 (MQ=255)
                                                                   tCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAgg                                                                         >  8:349245/1‑149 (MQ=255)
                                                                   tCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTAGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAgg                                                                         >  3:98990/1‑149 (MQ=255)
                                                                     aaCATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAACCCCAGAGAGGCATAGTAGGGCAGTGAAGGTa                                                                       >  4:299539/1‑149 (MQ=255)
                                                                     aaCATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGt                                                                        <  4:140953/148‑1 (MQ=255)
                                                                      aCATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTAt                                                                      >  4:122606/1‑149 (MQ=255)
                                                                        aTGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTa                                                                    <  8:133346/149‑1 (MQ=255)
                                                                        aTGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTa                                                                    <  3:122606/149‑1 (MQ=255)
                                                                           gACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCg                                                                 >  1:202400/1‑149 (MQ=255)
                                                                            aCATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGt                                                                <  5:337402/149‑1 (MQ=255)
                                                                              aTGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCa                                                              <  7:349245/149‑1 (MQ=255)
                                                                                gCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAg                                                             >  7:106173/1‑148 (MQ=255)
                                                                                 cACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAgg                                                            <  2:202400/148‑1 (MQ=255)
                                                                                   ctctGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGaaa                                                         <  7:78468/149‑1 (MQ=255)
                                                                                    tctGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAAt                                                        >  1:221887/1‑149 (MQ=255)
                                                                                        tCGGTTTGAACTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTAATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGATAGGCATAGTAGCGCAGTGAAGGTATTACCGTCAGGAAATGaa                                                     <  5:161266/148‑1 (MQ=255)
                                                                                        tCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAtagta                                                                                     <  6:153059/116‑1 (MQ=255)
                                                                                        tCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAtagta                                                                                     >  5:153059/1‑116 (MQ=255)
                                                                                        tCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAAc                                                    <  7:176476/149‑1 (MQ=255)
                                                                                                  gTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGcatca                                          <  5:321540/149‑1 (MQ=255)
                                                                                                   tttAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGcataat                                         >  3:226257/1‑149 (MQ=255)
                                                                                                      aaaTTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAg                                                             >  8:227697/1‑126 (MQ=255)
                                                                                                      aaaTTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAg                                                             <  7:227697/126‑1 (MQ=255)
                                                                                                      aaaTTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGAGGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGGAGTGAAGGTATTACCGTCAGGAAATGAACGTCTGCATGATCaa                                      >  4:9354/1‑149 (MQ=255)
                                                                                                       aaTTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCATGTCATCCATTGCCGGTAATTCTGCGTAATACGCTAAACCCAGAGAGTCATATTAGGGCAGTGAAGGTATTACCTTCAGGAAATGAACCTCTGAATCATCaa                                      >  7:251101/1‑148 (MQ=255)
                                                                                                           ttttAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCtt                                  <  8:106173/148‑1 (MQ=255)
                                                                                                             ttAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGc                               >  2:56095/1‑149 (MQ=255)
                                                                                                              tAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGc                               <  2:349502/148‑1 (MQ=255)
                                                                                                                aTGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCctc                            >  5:339617/1‑149 (MQ=255)
                                                                                                                 tGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCctct                           >  8:182774/1‑149 (MQ=255)
                                                                                                                    tctcTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTcc                        >  5:340923/1‑149 (MQ=255)
                                                                                                                           cAGTTGCGATTTGGCTTCCGCTTTAAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGt                 <  1:188309/149‑1 (MQ=255)
                                                                                                                           cAGTTGCGATTTGGCTTCCGCTTTAAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGt                 <  1:182706/149‑1 (MQ=255)
                                                                                                                             gTTGCTATTTGGCTTCCTCTTTCAGGGAATCCATTGCCGGTAATTGTGCGGAATACGCTAAACCCATAGATGCATAGTAGGGCAGTGAATGTATTACCGTCATGAAATGATACTCTGCATAATCAATCTTTGCCTCTTCCTCTACGTgg               <  8:264918/149‑1 (MQ=255)
                                                                                                                             gTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTgg               <  6:340923/149‑1 (MQ=255)
                                                                                                                             gTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGGGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTgg               <  5:242324/149‑1 (MQ=255)
                                                                                                                             gTTGCGATTTGGCTTCAGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTgg               >  8:149878/1‑149 (MQ=255)
                                                                                                                                  gaattGGCATCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCTTCAGGAAATGAACCACTGCATCATCAATCTATTCCACTTCCTCAACGTGGCTAAt          >  4:56846/4‑149 (MQ=255)
                                                                                                                                  gATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTaa           >  8:89437/1‑148 (MQ=255)
                                                                                                                                          ttCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGtt  >  4:224782/1‑149 (MQ=255)
                                                                                                                                            ccGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGtt  <  4:44555/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
AATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTT  >  NC_000913/1435281‑1435567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: