Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,920,796 T→C 17.5% A31A (GCA→GCG 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,7960TC17.5% 200.1 / 33.1 80A31A (GCA→GCGatpFF0 sector of membrane‑bound ATP synthase, subunit b
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (7/7);  total (40/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 = 9.53e-01

GATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAG  >  NC_000913/3920664‑3920940
                                                                                                                                    |                                                                                                                                                
gatTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcggc                                                                                                                                   >  1:947005/1‑148 (MQ=255)
             tCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCCGCCATACGTACTTCa                                                                                                                      >  1:118697/1‑148 (MQ=255)
                   ttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcaga                                                                                                                <  4:1056960/148‑1 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTaa                                                                                                                                          <  2:14770/116‑1 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTaa                                                                                                                                          >  1:14770/1‑116 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCa                                                                                                                                 <  4:446691/125‑1 (MQ=255)
                         tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGCCCGCCATTAATGGCAGCCa                                                                                                                                 >  3:446691/1‑125 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACaaa                                                                                             <  2:761457/148‑1 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACaaa                                                                                             >  1:761457/1‑148 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggca                                                                                                                                     <  4:959304/108‑2 (MQ=255)
                                      cggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggca                                                                                                                                     >  3:959304/1‑107 (MQ=255)
                                        gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                                                                       <  3:1038523/136‑1 (MQ=255)
                                        gtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                                                                       >  4:1038523/1‑136 (MQ=255)
                                              cttgCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtgg                                                                                     <  2:118697/145‑1 (MQ=255)
                                                    tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcgg                                                                                                                                    <  2:556093/95‑1 (MQ=255)
                                                    tttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcgg                                                                                                                                    >  1:556093/1‑95 (MQ=255)
                                                     ttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    <  2:1069513/94‑3 (MQ=255)
                                                     ttGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    >  1:1069513/1‑92 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTa                                                                                                                                           >  4:644281/1‑84 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTa                                                                                                                                           <  3:644281/84‑1 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctg                                                                                 >  1:751216/1‑142 (MQ=255)
                                                        caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctg                                                                                 <  2:751216/142‑1 (MQ=255)
                                                             tcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                        <  3:1086168/98‑1 (MQ=255)
                                                             tcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                        >  4:1086168/1‑98 (MQ=255)
                                                                  gtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaa                                                                                                      <  3:341502/111‑1 (MQ=255)
                                                                  gtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaa                                                                                                      >  4:341502/1‑111 (MQ=255)
                                                                   tcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATa                                                                                                                               <  4:264875/85‑1 (MQ=255)
                                                                   tcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATa                                                                                                                               >  3:264875/1‑85 (MQ=255)
                                                                   tcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    >  4:772686/1‑78 (MQ=255)
                                                                   tcCATATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    <  3:772686/80‑3 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                      <  2:721770/138‑1 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                      >  1:1052536/1‑138 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                      <  2:1052536/138‑1 (MQ=255)
                                                                       tATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                      >  1:721770/1‑138 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaac                                                                                                              <  3:580067/94‑1 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaac                                                                                                              >  4:580067/1‑94 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    <  4:174605/72‑3 (MQ=255)
                                                                           tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATggcag                                                                                                                                    >  3:174605/1‑70 (MQ=255)
                                                                             cTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                    <  4:907926/118‑1 (MQ=255)
                                                                             cTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGAAGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                    >  3:907926/1‑118 (MQ=255)
                                                                              tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                   >  2:479059/1‑134 (MQ=255)
                                                                              tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgttg                                                                   <  1:479059/134‑1 (MQ=255)
                                                                                    tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                    <  3:983688/95‑1 (MQ=255)
                                                                                    tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                    >  4:983688/1‑95 (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                                                  <  4:951052/121‑1 (MQ=255)
                                                                                                            ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGca                                                  >  3:951052/1‑121 (MQ=255)
                                                                                                              cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                    <  1:302209/85‑1 (MQ=255)
                                                                                                              cTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggc                                                                                    >  2:302209/1‑85 (MQ=255)
                                                                                                               tttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGata                    <  2:947005/148‑1 (MQ=255)
                                                                                                                ttttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGatat                   >  1:871375/1‑148 (MQ=255)
                                                                                                                  ttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAt                                                          <  3:382416/107‑1 (MQ=255)
                                                                                                                  ttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAt                                                          >  4:382416/1‑107 (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                                                                     >  3:66568/1‑90 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgt                                                                     <  4:66568/90‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGt                                                                >  2:248085/1‑95 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGt                                                                <  1:248085/95‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTa                               <  1:473502/128‑1 (MQ=255)
                                                                                                                        tttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTa                               >  2:473502/1‑128 (MQ=255)
                                                                                                                               aTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTa                                       <  1:289994/113‑1 (MQ=255)
                                                                                                                               aTGGCCGCCATTAATGGCAGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTa                                       >  2:289994/1‑113 (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)
                                                                                                                                    |                                                                                                                                                
GATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAG  >  NC_000913/3920664‑3920940

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: