Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,207,747 A→T 7.5% L25M (TTG→ATG) ‡ 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,7470AT 7.5% 189.9 / 4.3 66L25M (TTG→ATG) ‡PP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base T (3/2);  total (37/29)
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.78e-01

CGCCATACCTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGATGGC  >  NC_002947/4207599‑4207895
                                                                                                                                                    |                                                                                                                                                    
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:861692/150‑1 (MQ=255)
                        tCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTa                                                                                                                             >  1:1662441/1‑150 (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:967009/1‑150 (MQ=255)
                                                         aCGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGa                                                                                            >  1:948854/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                                    <  1:1335948/150‑1 (MQ=255)
                                                                                                                 gTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAAt                                    >  2:351343/1‑150 (MQ=255)
                                                                                                                 gTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAAt                                    >  2:1335948/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                                                             >  1:1241850/1‑113 (MQ=255)
                                                                                                                             ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTc                                                             <  2:1241850/113‑1 (MQ=255)
                                                                                                                             ggCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCaga                        >  2:814586/1‑150 (MQ=255)
                                                                                                                              gCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCagag                       >  1:1133846/1‑150 (MQ=255)
                                                                                                                                   tgtcgttttCGCCACCAAGGTGGCGATGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGccc                  <  1:857418/144‑1 (MQ=255)
                                                                                                                                    ttCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAgg                                                <  2:1700941/119‑1 (MQ=255)
                                                                                                                                    ttCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAgg                                                >  1:1700941/1‑119 (MQ=255)
                                                                                                                                       aacctggCCACCTTGGTGGCCAGGGTCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCctct              >  2:594899/7‑150 (MQ=255)
                                                                                                                                        gTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCtt             >  2:607826/1‑150 (MQ=255)
                                                                                                                                        gTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCtt             >  1:1599031/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)
                                                                                                                                                    |                                                                                                                                                    
CGCCATACCTTGGGTTTGAAAACCTCGGGGTCGGAGGTCAAGGTCAGTGTTGCATTGACGAAGCGGTCTTCGGTCTCGTTCGCGATCCGAAGTTCGCGAATTACGGGAATGGCGTTCTGGAAGTCGGCAAGGTTCAGTTTCGCCACCAAGGTGGCCAGGATCTTGACTTCTTTAGGCGGGATTTCTTCTTGAGCTTTGAGGGGCTGAACGATGTCGTCCATGCGCAACCTTCCTTGTCGTAATGCCCAAGGGAGTGCCGGAATTGGCCGGCCAGAGGGCCCCTCTTGAGCTGATGGC  >  NC_002947/4207599‑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: ‑