Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,920,809 G→A 10.3% P27L (CCG→CTG)  atpF ← F0 sector of membrane‑bound ATP synthase, subunit b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,920,8090GA10.3% 236.2 / 16.1 79P27L (CCG→CTG) atpFF0 sector of membrane‑bound ATP synthase, subunit b
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base A (4/4);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

GATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACG  >  NC_000913/3920664‑3920954
                                                                                                                                                 |                                                                                                                                                 
gatTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcggc                                                                                                                                                 >  1:947005/1‑148 (MQ=255)
             tCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCCGCCATACGTACTTCa                                                                                                                                    >  1:118697/1‑148 (MQ=255)
                   ttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcaga                                                                                                                              <  4:1056960/148‑1 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCa                                                                                                                                               <  4:446691/125‑1 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGCCCGCCATTAATGGCAGCCa                                                                                                                                               >  3:446691/1‑125 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACaaa                                                                                                           >  1:761457/1‑148 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACaaa                                                                                                           <  2:761457/148‑1 (MQ=255)
                                        gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                                                                                     >  4:1038523/1‑136 (MQ=255)
                                        gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                                                                                     <  3:1038523/136‑1 (MQ=255)
                                              cttgCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtgg                                                                                                   <  2:118697/145‑1 (MQ=255)
                                                    tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcgg                                                                                                                                                  <  2:556093/95‑1 (MQ=255)
                                                    tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcgg                                                                                                                                                  >  1:556093/1‑95 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctg                                                                                               >  1:751216/1‑142 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctg                                                                                               <  2:751216/142‑1 (MQ=255)
                                                             tcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                                      >  4:1086168/1‑98 (MQ=255)
                                                             tcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                                      <  3:1086168/98‑1 (MQ=255)
                                                                  gtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaa                                                                                                                    <  3:341502/111‑1 (MQ=255)
                                                                  gtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaa                                                                                                                    >  4:341502/1‑111 (MQ=255)
                                                                   tcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATa                                                                                                                                             >  3:264875/1‑85 (MQ=255)
                                                                   tcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATa                                                                                                                                             <  4:264875/85‑1 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                                    <  2:721770/138‑1 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                                    <  2:1052536/138‑1 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                                    >  1:1052536/1‑138 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                                    >  1:721770/1‑138 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaac                                                                                                                            <  3:580067/94‑1 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaac                                                                                                                            >  4:580067/1‑94 (MQ=255)
                                                                             cTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                                  <  4:907926/118‑1 (MQ=255)
                                                                             cTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGAAGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                                  >  3:907926/1‑118 (MQ=255)
                                                                              tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                                 <  1:479059/134‑1 (MQ=255)
                                                                              tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                                 >  2:479059/1‑134 (MQ=255)
                                                                                    tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                                  >  4:983688/1‑95 (MQ=255)
                                                                                    tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                                  <  3:983688/95‑1 (MQ=255)
                                                                                         aaGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGa                                                                                                      >  2:1047/1‑102 (MQ=255)
                                                                                         aaGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGa                                                                                                      <  1:1047/102‑1 (MQ=255)
                                                                                            cAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                                  >  2:625460/1‑87 (MQ=255)
                                                                                            cAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                                  <  1:625460/87‑1 (MQ=255)
                                                                                              aGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTAtt                                                   >  1:386932/1‑148 (MQ=255)
                                                                                               ggCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                                 >  3:239401/1‑117 (MQ=255)
                                                                                               ggCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                                 <  4:239401/117‑1 (MQ=255)
                                                                                                 ccGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGGCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAg                                                >  3:629838/1‑148 (MQ=255)
                                                                                                          aTTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGGCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTc                                       <  4:629838/148‑1 (MQ=255)
                                                                                                            ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                                                                                         >  4:1099969/1‑96 (MQ=255)
                                                                                                            ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                                                                                         <  3:1099969/96‑1 (MQ=255)
                                                                                                            ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGca                                                                >  3:951052/1‑121 (MQ=255)
                                                                                                            ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGca                                                                <  4:951052/121‑1 (MQ=255)
                                                                                                              cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                                  >  2:302209/1‑85 (MQ=255)
                                                                                                              cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                                  <  1:302209/85‑1 (MQ=255)
                                                                                                               tttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGata                                  <  2:947005/148‑1 (MQ=255)
                                                                                                                ttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGatat                                 >  1:871375/1‑148 (MQ=255)
                                                                                                                  ttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAt                                                                        >  4:382416/1‑107 (MQ=255)
                                                                                                                  ttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAt                                                                        <  3:382416/107‑1 (MQ=255)
                                                                                                                      cGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAt                                                                                      >  2:256470/1‑89 (MQ=255)
                                                                                                                      cGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAt                                                                                      <  1:256470/89‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcct                                                                                                <  1:341583/77‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcct                                                                                                >  2:341583/1‑77 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgt                                                                                   <  4:66568/90‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgt                                                                                   >  3:66568/1‑90 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGt                                                                              >  2:248085/1‑95 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGt                                                                              <  1:248085/95‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTa                                             >  2:473502/1‑128 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTa                                             <  1:473502/128‑1 (MQ=255)
                                                                                                                               aTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTa                                                     <  1:289994/113‑1 (MQ=255)
                                                                                                                               aTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTa                                                     >  2:289994/1‑113 (MQ=255)
                                                                                                                                 ggccgCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAAt                                                            >  4:265487/5‑104 (MQ=255)
                                                                                                                                 ggccgCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGa                                    <  4:608226/124‑1 (MQ=255)
                                                                                                                                 ggCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAg                <  2:386932/148‑1 (MQ=255)
                                                                                                                                 ggCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAAt                                                            <  3:265487/104‑1 (MQ=255)
                                                                                                                                 ggCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGa                                    >  3:608226/1‑128 (MQ=255)
                                                                                                                                   cTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                                 >  1:274347/1‑81 (MQ=255)
                                                                                                                                   cTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgctg                                                                                 <  2:274347/81‑1 (MQ=255)
                                                                                                                                     gCCATTAATGGCGGCCATACGTACTTCATGGAGAACAGAACGAACAGGACAAACGCGa                                                                                                      <  3:946844/58‑1 (MQ=255)
                                                                                                                                     gCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGa                                                                                                      >  4:946844/1‑58 (MQ=255)
                                                                                                                                     gCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTa                                                 >  2:242250/1‑111 (MQ=255)
                                                                                                                                     gCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTa                                                 <  1:242250/111‑1 (MQ=255)
                                                                                                                                      ccATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGatat                                 <  1:899015/126‑1 (MQ=255)
                                                                                                                                      ccATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGatat                                 >  2:899015/1‑126 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCctc                                                       <  4:844605/101‑1 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCctc                                                       >  3:844605/1‑101 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTg                               >  2:654003/1‑125 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTg                               <  1:654003/125‑1 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGct                             <  4:965494/127‑1 (MQ=255)
                                                                                                                                         ttAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGct                             >  3:965494/1‑127 (MQ=255)
                                                                                                                                               gcggcCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTAcg  <  1:1068587/148‑1 (MQ=255)
                                                                                                                                               gcggcCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTAcg  >  2:1068587/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACG  >  NC_000913/3920664‑3920954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: