Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,409,076 A→G 6.7% F428L (TTC→CTC)  PP_3893 ← DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,409,0760AG 6.7% 290.5 / 8.5 104F428L (TTC→CTC) PP_3893DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (56/41);  new base G (3/4);  total (59/45)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 105 total aligned reads displayed.

TTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGCTCGATCTCGCCCG  >  NC_002947/4408934‑4409225
                                                                                                                                              |                                                                                                                                                     
ttCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGATGGCAAACCGGGCGAACTCGccc                                                                                                                                                >  2:2359654/1‑150 (MQ=255)
   ggTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGGGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAGTTCGcccggt                                                                                                                                             <  1:1537775/150‑4 (MQ=255)
         ttCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGAtgctgc                                                                                                                                       >  2:647910/1‑150 (MQ=255)
         ttCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGAtgctgc                                                                                                                                       >  1:1802178/1‑150 (MQ=255)
               cggcAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAgcagg                                                                                                                                 <  2:1461651/150‑1 (MQ=255)
               cggcAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAgcagg                                                                                                                                 <  2:2252104/150‑1 (MQ=255)
                          tCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGAc                                                                                                                      >  1:187771/1‑150 (MQ=255)
                          tCAGAACCTCTCACCTGGATGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGAc                                                                                                                      >  2:825600/1‑150 (MQ=255)
                             gAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGa                                                                                                                       >  2:809466/1‑146 (MQ=255)
                             gAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAgg                                                                                                                   >  2:869192/1‑150 (MQ=255)
                              aaCCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGACGACAGgc                                                                                                                  <  2:1802178/150‑1 (MQ=255)
                                ccTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGgcgc                                                                                                                <  1:587425/150‑1 (MQ=255)
                                 ctctCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGcgcg                                                                                                               >  2:435263/1‑150 (MQ=255)
                                   ctcCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGAGTg                                                                                                             >  1:409502/1‑150 (MQ=255)
                                    tcCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGa                                                                                                            <  2:550638/150‑1 (MQ=255)
                                     ccccTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTg                                                                                                                                              <  1:1041120/115‑1 (MQ=255)
                                     ccccTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTg                                                                                                                                              >  2:1041120/1‑115 (MQ=255)
                                      cccTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACg                                                                                                          <  1:97879/150‑1 (MQ=255)
                                      cccTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACg                                                                                                          >  2:598065/1‑150 (MQ=255)
                                      cccTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACg                                                                                                          >  2:1695309/1‑150 (MQ=255)
                                       ccTGGCTGGAAGTTTTTCAGCATCGAGCGGGAGGACCTGCGAGCCGCTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACgg                                                                                                         <  1:264765/150‑1 (MQ=255)
                                        ctggctggAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCCGCCTGCTGCTGCTCGCGAGCGCCGTCGTAGTTGGCAAAACGGGCGAACAGGACATAATCCTTCACCAGGCAGTGGCCGACAGGCGCGTGACggg                                                                                                        >  1:540706/1‑149 (MQ=255)
                                            ctggAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCtt                                                                                                    <  1:598065/150‑1 (MQ=255)
                                               gAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGAc                                                                                                 <  2:413547/150‑1 (MQ=255)
                                               gAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGAc                                                                                                 >  1:928263/1‑150 (MQ=255)
                                               gAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGAc                                                                                                 >  2:1628645/1‑150 (MQ=255)
                                               gAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGAc                                                                                                 <  1:924772/150‑1 (MQ=255)
                                                  gttgttCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACa                                                                                                                     >  1:608046/1‑127 (MQ=255)
                                                  gttgttCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACa                                                                                                                     <  2:608046/127‑1 (MQ=255)
                                                     gttCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATc                                                                                           >  1:1066723/1‑150 (MQ=255)
                                                         aGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAGTTCGcccggt                                                                                                                                             >  1:125963/1‑93 (MQ=255)
                                                         aGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAGTTCGcccggt                                                                                                                                             <  2:125963/96‑4 (MQ=255)
                                                          gCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGAtct                                                                                      >  2:2316287/1‑150 (MQ=255)
                                                           cATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGAt                                                                                                                                            <  2:115598/95‑1 (MQ=255)
                                                           cATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGAt                                                                                                                                            >  1:115598/1‑95 (MQ=255)
                                                            aTCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCg                                                                                    >  2:448951/1‑150 (MQ=255)
                                                               gAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCGCGGCGATCTCGGTa                                                                                 >  1:496982/1‑150 (MQ=255)
                                                                     gaggaCCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAAcgccg                                                                           >  2:1349045/1‑150 (MQ=255)
                                                                     gaggaCCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACACCGATCTCGGTCAcgccg                                                                           >  2:382814/1‑150 (MQ=255)
                                                                         aCCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGctct                                                                       >  2:2339079/1‑150 (MQ=255)
                                                                          ccTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACa                                                                                                                     >  1:1021761/1‑103 (MQ=255)
                                                                          ccTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACa                                                                                                                     <  2:1021761/103‑1 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAGTTCGcccggt                                                                                                                                             <  1:1154440/78‑4 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAGTTCGcccggt                                                                                                                                             >  2:1154440/1‑75 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGCAGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGATTGACTACCTGGAGCCCGGTAACGCCGCTCTgg                                                                     >  1:1021906/1‑150 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCgg                                                                                   >  1:1609309/1‑136 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCgg                                                                                   <  2:1609309/136‑1 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCg                                                                                    >  1:863566/1‑135 (MQ=255)
                                                                           cTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCg                                                                                    <  2:863566/135‑1 (MQ=255)
                                                                            tGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGc                                                                    >  2:1003286/1‑150 (MQ=255)
                                                                            tGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGc                                                                    >  1:1939669/1‑150 (MQ=255)
                                                                              cGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGccc                                                                  >  2:449574/1‑150 (MQ=255)
                                                                                    ggTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGANGGTGCTTGACCACATCGATNTCGGTAACGCCGCTCTGGCCCAAGTcc                                                            <  1:1628645/150‑1 (MQ=255)
                                                                                       tGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGaa                                                         <  2:1066723/150‑1 (MQ=255)
                                                                                           ctgctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCAAGAGGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCgg                                                     >  1:636116/1‑150 (MQ=255)
                                                                                             gctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTgg                                                                     <  2:1633640/132‑1 (MQ=255)
                                                                                             gctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTgg                                                                     >  1:1633640/1‑132 (MQ=255)
                                                                                             gctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCa                                                   >  2:2005411/1‑150 (MQ=255)
                                                                                              ctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACAt                                                                                            >  2:2177325/1‑108 (MQ=255)
                                                                                              ctgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACAt                                                                                            <  1:2177325/108‑1 (MQ=255)
                                                                                               tgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGAtc                                                                                       <  2:190430/112‑1 (MQ=255)
                                                                                               tgcctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGAtc                                                                                       >  1:190430/1‑112 (MQ=255)
                                                                                                 cctgctgcTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATcgc                                               >  1:1628307/1‑150 (MQ=255)
                                                                                                   tgctgctgcTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCa                                                   >  2:1984584/1‑144 (MQ=255)
                                                                                                   tgctgctgcTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCa                                                   <  1:1984584/144‑1 (MQ=255)
                                                                                                    gctgctgctCGCGGTCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCCGATGCTGCAGCAGGCAGTGGCCGATAGGCGCGTGACGGTGCTCGACCCCATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGt                                            <  1:841844/150‑1 (MQ=255)
                                                                                                    gctgctgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGt                                            <  1:1349045/150‑1 (MQ=255)
                                                                                                    gctgctgctCGCGAGCTCCGGCGGAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGt                                            <  2:540706/150‑1 (MQ=255)
                                                                                                     ctgctgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGccc                                                                                                                                                <  2:1125004/49‑1 (MQ=255)
                                                                                                     ctgctgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGccc                                                                                                                                                >  1:1125004/1‑49 (MQ=255)
                                                                                                      tgctgctCGCGAGCTCCGGCGTAGTTGGCAAGCCGGGCGAACTCGCTCCCAGTCTGCAGTAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTgg                                          <  2:65081/150‑1 (MQ=255)
                                                                                                      tgctgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTgg                                          <  2:2293791/150‑1 (MQ=255)
                                                                                                      ggcggttCGCGAGCTGCGGCGTAGTTGGCAGACCGGGCGAACTCGCCCTGATGGTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGATTGACCACATCGATCTCGGTAACGCCGCTCTGGGCCAAGTCCGAATCGGCATCGCGGTgg                                          <  2:745734/144‑1 (MQ=255)
                                                                                                         tgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGaa                                                         >  2:727253/1‑132 (MQ=255)
                                                                                                         tgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGaa                                                         <  1:727253/132‑1 (MQ=255)
                                                                                                         tgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCg                                       >  2:197533/1‑150 (MQ=255)
                                                                                                         tgctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCg                                       >  1:1166904/1‑150 (MQ=255)
                                                                                                          gctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGa                                      <  1:1276490/150‑1 (MQ=255)
                                                                                                           ctCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGTGTGACGGTGCTTGACCACATCGATCTCGTTAACGCCGCCCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGAt                                     >  1:37216/1‑150 (MQ=255)
                                                                                                            tCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATc                                    >  2:2357340/1‑150 (MQ=255)
                                                                                                                    tCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGgatgat                            >  1:1112860/1‑150 (MQ=255)
                                                                                                                       ggCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTAGGCGATCAGGATGATGTc                         >  1:124826/1‑150 (MQ=255)
                                                                                                                        gCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACTACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCg                                       <  1:86816/135‑1 (MQ=255)
                                                                                                                        gCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACTACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCg                                       >  2:86816/1‑135 (MQ=255)
                                                                                                                        gCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCaa                                                                <  2:1012247/110‑1 (MQ=255)
                                                                                                                        gCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCaa                                                                >  1:1012247/1‑110 (MQ=255)
                                                                                                                           tAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCg                     <  2:1939669/150‑1 (MQ=255)
                                                                                                                           tAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCg                     >  2:566235/1‑150 (MQ=255)
                                                                                                                            aGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGt                    <  1:869192/150‑1 (MQ=255)
                                                                                                                             gTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTc                   >  2:325050/1‑150 (MQ=255)
                                                                                                                              ttGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCgg                                                                                   <  2:1786500/85‑1 (MQ=255)
                                                                                                                              ttGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCt                  >  1:292391/1‑150 (MQ=255)
                                                                                                                                ggCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTg                <  1:825600/150‑1 (MQ=255)
                                                                                                                                 gCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGc               >  1:1640629/1‑150 (MQ=255)
                                                                                                                                 gCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGc               <  1:435263/150‑1 (MQ=255)
                                                                                                                                  gcatcagggCGAGTTCGCCCGGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGc               >  1:1902731/7‑149 (MQ=255)
                                                                                                                                  gcatcagggCGAGTTCGCCCGGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGc               <  2:1902731/143‑1 (MQ=255)
                                                                                                                                   aaaCCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGCTc             <  1:449574/150‑1 (MQ=255)
                                                                                                                                      ccGGGCGAACTCGCCCGGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGCTCGAt          <  2:1085997/150‑1 (MQ=255)
                                                                                                                                              aCTCGCCCTGATGCTGCAGCAGGCCGGGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCCCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGCTCGATCTCGCCCg  >  1:933500/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGGCGAACTCGCCCTGATGCTGCAGCAGGCAGTGGCCGACAGGCGCGTGACGGTGCTTGACCACATCGATCTCGGTAACGCCGCTCTGGCCCAAGTCCGAATCGGCATCGCGGTGGGCGATCAGGATGATGTCGGCGTCTTGCTCGATCTCGCCCG  >  NC_002947/4408934‑4409225

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