Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,207,746 A→T 7.7% L25* (TTG→TAG) ‡ PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,207,7460AT 7.7% 189.2 / 5.4 65L25* (TTG→TAG) ‡PP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base T (3/2);  total (37/28)
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.12e-01

TCCGCCATACCTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGATGGC  >  NC_002947/4207597‑4207895
                                                                                                                                                     |                                                                                                                                                     
tCCGCCATGCCCTAGGTGGGAAGACCACGTGGGCGGAGGTCAATGTCAGGGTTGCAGTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGccacca                                                                                                                                                       >  2:1552505/1‑150 (MQ=255)
  cGCCATACCTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAg                                                                                                                                                     <  2:167320/150‑1 (MQ=255)
          cTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCTTGGTGGcgaa                                                                                                                                             <  2:198982/150‑4 (MQ=255)
                   gAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGAc                                                                                                                                    >  1:1072388/1‑150 (MQ=255)
                    aaaaCCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTCACGCAGAGGCCTACGAACTCATTCGCCAACCTGTGCACGCACAGTACATAAAGGGCTTTGTGGTGTTCTGCGAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGAct                                                                                                                                   >  2:1161103/1‑150 (MQ=255)
                      aaCCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACGAATGTGGCCAGGATCTTGActtc                                                                                                                                 >  2:439110/1‑150 (MQ=255)
                        ccTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGCATGACGACATAATATTATGCATTTTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGActtctt                                                                                                                               >  1:1869501/1‑150 (MQ=255)
                         cTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGGCTTACGGGAATGGCGTTCTGGAAGTCGGCGAGGGTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCttt                                                                                                                              >  1:602602/1‑150 (MQ=255)
                         cTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCttt                                                                                                                              >  2:1128621/1‑150 (MQ=255)
                          tCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTa                                                                                                                             >  1:1662441/1‑150 (MQ=255)
                          tCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTa                                                                                                                             <  1:861692/150‑1 (MQ=255)
                            ggggTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGACGGTGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAgg                                                                                                                           <  1:1079363/150‑1 (MQ=255)
                                 cGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGa                                                                                                                      <  1:1330942/150‑1 (MQ=255)
                                  ggAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGAt                                                                                                                     <  2:1623247/150‑1 (MQ=255)
                                     ggtcaaggtcaGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATttc                                                                                                                  <  2:1662441/150‑1 (MQ=255)
                                      gtcaaggtcaGTGTTGCATTGACGAAGCGGTGTTAAGAGCCGTCCGCAATCCGTGGTTCGAGCATTACGGTAAAGGCGTTCTGGAAGACGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTCCTTTAGGCGGGATttct                                                                                                                 >  1:1471558/1‑150 (MQ=255)
                                       tcaaggtcaGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATttctt                                                                                                                <  2:262329/150‑1 (MQ=255)
                                            gtcaGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCGACCAAGGTGGCCAGGACCTTGACTTCTTTAGGCGGGATTTCTTCTTGa                                                                                                           <  2:375393/150‑1 (MQ=255)
                                               aGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCt                                                                                                        >  1:345856/1‑150 (MQ=255)
                                                      cATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCAAGGACACTCACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAggg                                                                                                 >  1:422537/1‑150 (MQ=255)
                                                           aCGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGa                                                                                            >  1:948854/1‑150 (MQ=255)
                                                           aCGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGa                                                                                            >  1:967009/1‑150 (MQ=255)
                                                           aCGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAATTGGCGAATTACGGGCATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGa                                                                                            >  1:152405/1‑150 (MQ=255)
                                                            cGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGaa                                                                                           <  2:948854/150‑1 (MQ=255)
                                                                gCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGAt                                                                                       <  1:1145023/150‑1 (MQ=255)
                                                                      ttCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATgtcgtc                                                                                 >  1:606143/1‑150 (MQ=255)
                                                                       tCGGTCTCGTTGGCGATCCGAAGTTCGCGACCCATAAAAGTGGCTTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTcc                                                                                <  2:912296/150‑1 (MQ=255)
                                                                          gTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATg                                                                             >  1:1395912/1‑150 (MQ=255)
                                                                                ttCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAAc                                                                       <  1:1163088/150‑1 (MQ=255)
                                                                                   gcgATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAAcctt                                                                    <  2:1471558/150‑1 (MQ=255)
                                                                                     gATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCCAccttcc                                                                  <  2:967009/150‑1 (MQ=255)
                                                                                       tCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAAccttcctt                                                                <  1:1112337/150‑1 (MQ=255)
                                                                                         cGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGt                                                              >  1:1832281/1‑150 (MQ=255)
                                                                                              ttCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTaa                                                         <  1:680570/150‑1 (MQ=255)
                                                                                                 gcgAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCATTGTCGTAATGc                                                      <  2:1395912/150‑1 (MQ=255)
                                                                                                      ttACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAg                                                 >  2:89831/1‑150 (MQ=255)
                                                                                                       tACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAgg                                                <  2:1832281/150‑1 (MQ=255)
                                                                                                            gAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTg                                           >  1:1598914/1‑150 (MQ=255)
                                                                                                            gAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTg                                           >  2:1281161/1‑150 (MQ=255)
                                                                                                             aaTGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGc                                          <  2:345856/150‑1 (MQ=255)
                                                                                                                   gTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAAt                                    >  2:1335948/1‑150 (MQ=255)
                                                                                                                   gTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAAt                                    <  1:1335948/150‑1 (MQ=255)
                                                                                                                   gTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAAt                                    >  2:351343/1‑150 (MQ=255)
                                                                                                                      cTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGCTCATGTTGGCTTTAGGCGGGACTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTgg                                 <  1:1552505/150‑1 (MQ=255)
                                                                                                                       tGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTAGTCGTAATGCCCAAGGGAGTGCCGGAATTggc                                <  2:1072388/150‑1 (MQ=255)
                                                                                                                        ggAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTggcc                               >  2:372445/1‑150 (MQ=255)
                                                                                                                         gAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTggccg                              >  1:1028805/1‑150 (MQ=255)
                                                                                                                          aaGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTggccgg                             >  2:1925087/1‑150 (MQ=255)
                                                                                                                               ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAAc                                                                                          <  1:1198909/84‑1 (MQ=255)
                                                                                                                               ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAAc                                                                                          >  2:1198909/1‑84 (MQ=255)
                                                                                                                               ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTc                                                             <  2:1241850/113‑1 (MQ=255)
                                                                                                                               ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTc                                                             >  1:1241850/1‑113 (MQ=255)
                                                                                                                               ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCaga                        >  2:814586/1‑150 (MQ=255)
                                                                                                                                gCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCagag                       >  1:1133846/1‑150 (MQ=255)
                                                                                                                                     tgtcgttttCGCCACCAAGGTGGCGATGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGccc                  <  1:857418/144‑1 (MQ=255)
                                                                                                                                      ttCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAgg                                                >  1:1700941/1‑119 (MQ=255)
                                                                                                                                      ttCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAgg                                                <  2:1700941/119‑1 (MQ=255)
                                                                                                                                         aacctggCCACCTTGGTGGCCAGGGTCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCctct              >  2:594899/7‑150 (MQ=255)
                                                                                                                                          gTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCtt             >  1:1599031/1‑150 (MQ=255)
                                                                                                                                          gTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCtt             >  2:607826/1‑150 (MQ=255)
                                                                                                                                          gTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCACAAGGGAATGCCGTAATTGGCCGGCCAGAGGGCCCCTCtt             >  2:310735/1‑150 (MQ=255)
                                                                                                                                           tttCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCAAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTg            >  1:422793/1‑150 (MQ=255)
                                                                                                                                            ctggCCACCTTGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGa           >  1:1587257/4‑150 (MQ=255)
                                                                                                                                             tggCCACCTTGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCaga                        <  1:1285122/134‑1 (MQ=255)
                                                                                                                                             tggCCACCTTGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAg          >  1:642161/3‑150 (MQ=255)
                                                                                                                                               gCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCt        >  1:1742284/1‑150 (MQ=255)
                                                                                                                                                  accaAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGAt     <  2:1869501/150‑1 (MQ=255)
                                                                                                                                                     aaGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGATGGc  <  2:683796/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGCCATACCTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGATGGC  >  NC_002947/4207597‑4207895

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