Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,824,196 A→G 6.4% F2347L (TTC→CTC)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,824,1960AG 6.4% 365.8 / 9.7 126F2347L (TTC→CTC) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (49/68);  new base G (5/3);  total (55/71)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 130 total aligned reads displayed.

CTGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAG  >  NC_002947/4824047‑4824344
                                                                                                                                                     |                                                                                                                                                    
cTGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGaa                                                                                                                                                      >  1:2398191/1‑150 (MQ=255)
cTGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGaa                                                                                                                                                      >  1:658137/1‑150 (MQ=255)
cTGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGaa                                                                                                                                                      <  2:2045661/150‑1 (MQ=255)
     cgcgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                 >  1:69486/1‑150 (MQ=255)
        gcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCTGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCAt                                                                                                                                              <  2:1601126/150‑1 (MQ=255)
        gcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCAt                                                                                                                                              <  2:1601094/150‑1 (MQ=255)
         cAGTACCACGCGCGCACCACACAAGAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAAGCATCGCTGCGCATGCCGAATGCGGGGATCTCCGCCTGGCAGTGCATGGCGAGTTCGCCATg                                                                                                                                             <  2:704059/150‑1 (MQ=255)
          aGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATgg                                                                                                                                            <  2:2059926/150‑1 (MQ=255)
           gTACCACGCGCGCACCACACGACAGTGGCACCCACAGGCGCTCGCTGGGAGCGACGAAATTGATCGAATAGAAATGCAGCTCGCAATCACCGCTGCGCAGGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcc                                                                                                                                           <  1:1542675/150‑3 (MQ=255)
             aCCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTg                                                                                                                                         <  1:1049930/150‑1 (MQ=255)
             aCCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTg                                                                                                                                         <  1:300568/150‑1 (MQ=255)
                aCGCGCGCACCACACAACAGCGGCACCCACAGGCGCACACTCGCCCGGTCCACAAGGATCGATAATAAATGGAGCCCACAATCATCGCTACGCATGCCGAAAGCGGCGATCACCGCCTGGCAGTGGATGGCGAACTCGCCATGGCTGacc                                                                                                                                      >  1:1225800/1‑150 (MQ=255)
                       caccacACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCGCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGc                                                                                                                               <  2:18714/150‑1 (MQ=255)
                          cacacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCt                                                                                                                            >  2:394077/1‑150 (MQ=255)
                            cacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTg                                                                                                                          <  2:1670099/150‑1 (MQ=255)
                             acaacaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTgg                                                                                                                         >  1:2428259/1‑150 (MQ=255)
                               aacaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGt                                                                                                                       >  1:2209175/1‑150 (MQ=255)
                                acaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGtt                                                                                                                      >  1:301924/1‑150 (MQ=255)
                                   gCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGc                                                                                                                   <  2:2209175/150‑1 (MQ=255)
                                    cGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTNGcc                                                                                                                  <  2:1520145/150‑1 (MQ=255)
                                       cACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGt                                                                                                               >  2:1414780/1‑150 (MQ=255)
                                        aCCCACAGGCGCTCGCCGGCACCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACGCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                              <  2:778034/150‑1 (MQ=255)
                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCa                                                                                                                                               >  1:2015211/1‑115 (MQ=255)
                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCCCGCCCTTGGGTTTGCCGGTGGa                                                                                                            <  1:1253845/150‑1 (MQ=255)
                                          ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGa                                                                                                            >  2:1266282/1‑150 (MQ=255)
                                           cacaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAg                                                                                                           <  1:1815607/150‑1 (MQ=255)
                                            acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGGATACAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGAGGCGGTCTCCGCCTTTCAGTGTAGGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGatc                                                                                                          <  2:398349/150‑3 (MQ=255)
                                            acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                          >  1:2279325/1‑150 (MQ=255)
                                            acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                          <  1:2080243/150‑1 (MQ=255)
                                            acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCAAAAAATTTATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                          <  1:1332604/150‑1 (MQ=38)
                                                gcgcTCGCTGGCATCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGa                                                                                                      <  2:301924/150‑1 (MQ=255)
                                                 cgctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCCCCAAGTCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                     <  1:1023857/150‑1 (MQ=255)
                                                 cgctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                     <  2:2279325/150‑1 (MQ=255)
                                                   ctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgt                                                                                                   <  1:1617858/150‑1 (MQ=255)
                                                     cgctGGCAGCGTCGATCTTGAACGAATAGAAATTCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                 <  2:1097456/150‑1 (MQ=255)
                                                          gCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATc                                                                                            >  2:43003/1‑150 (MQ=255)
                                                           cAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                           <  2:39188/150‑1 (MQ=255)
                                                           cAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                           <  2:776510/150‑1 (MQ=255)
                                                            aGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg                                                                                          <  1:2357906/150‑1 (MQ=255)
                                                             gCGTCGAAATTGATAGAATACAAATGCAGCTCGCAATCAGCGCTGCGCATGCCGAATGCGGCGATCAGCGCCTGGCAGTGCAGGGCGAACCCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGCAGCCTGAGGTGTAGATCAgc                                                                                         <  2:356638/150‑2 (MQ=255)
                                                                tCGAAATTGATCGAATAGAAATGCAGCTCGCAATCAGCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTaa                                                                                      <  1:1081296/150‑1 (MQ=255)
                                                                  gAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGc                                                                                    <  1:1095784/150‑1 (MQ=255)
                                                                    aaTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCt                                                                                  <  1:2451413/150‑1 (MQ=255)
                                                                     aTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAc                                                                                                                                                     >  1:1284448/1‑82 (MQ=255)
                                                                         aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCg                                                                                                                 <  1:2189525/114‑1 (MQ=255)
                                                                         aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCg                                                                                                                 >  2:2189525/1‑114 (MQ=255)
                                                                         aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgc                                                                             >  2:837736/1‑150 (MQ=255)
                                                                          tCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgct                                                                            >  2:361212/1‑150 (MQ=255)
                                                                             aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcac                                                                                                                                          >  2:2062524/1‑82 (MQ=255)
                                                                             aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcac                                                                                                                                          <  1:2062524/85‑4 (MQ=255)
                                                                             aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                     <  1:1823709/122‑1 (MQ=255)
                                                                              aTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCgg                                                                        >  2:2207429/1‑150 (MQ=255)
                                                                                 gAAATGCAGCTCGCAATCATCGCTGCGCATGCCGACTGTGGCGCTCACCGCCTGGCAGTGCATGGCAGACGGGCCATGGTTGTGCATCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGGAGCTCAGGTAAGCCTGGTGCTGCGGCAg                                                                     >  1:1023842/1‑150 (MQ=255)
                                                                                       cAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               <  2:69486/150‑1 (MQ=255)
                                                                                       cAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               <  1:1199091/150‑1 (MQ=255)
                                                                                       cAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               <  1:1195993/150‑1 (MQ=255)
                                                                                        aGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAgg                                                              <  1:1266282/150‑1 (MQ=255)
                                                                                           tCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCAGCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTg                                                           >  2:1064695/1‑150 (MQ=255)
                                                                                           tCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTg                                                                                 >  2:1591940/1‑128 (MQ=255)
                                                                                           tCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTg                                                           >  1:672742/1‑150 (MQ=255)
                                                                                            cGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGt                                                          <  1:1135802/150‑1 (MQ=255)
                                                                                             gCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTc                                                         <  2:672742/150‑1 (MQ=255)
                                                                                              cAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgctgc                                                                          <  2:1019233/132‑1 (MQ=255)
                                                                                              cAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgctgc                                                                          >  1:1019233/1‑132 (MQ=255)
                                                                                               aaTCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAGTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                           >  2:1283508/1‑114 (MQ=255)
                                                                                               aaTCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                           <  1:1283508/114‑1 (MQ=255)
                                                                                                      gCTGCGCATGCCGAATGCGGCGATCAGCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGc                                                <  1:1064695/150‑1 (MQ=255)
                                                                                                      gCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGAGACCTCGGCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGc                                                >  1:1637086/1‑150 (MQ=255)
                                                                                                          cgcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCg                                            >  2:694639/1‑150 (MQ=255)
                                                                                                             aTGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTg                                         <  1:1414780/150‑1 (MQ=255)
                                                                                                              tGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGt                                                                               >  2:900858/1‑111 (MQ=255)
                                                                                                               gCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGTGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTa                                       <  2:1276168/150‑1 (MQ=255)
                                                                                                               gCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTa                                       <  2:2398191/150‑1 (MQ=255)
                                                                                                                ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTg                                                                                 >  2:2401603/1‑107 (MQ=255)
                                                                                                                ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTg                                                                                 <  1:2401603/107‑1 (MQ=255)
                                                                                                                 cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                 >  1:810262/1‑90 (MQ=255)
                                                                                                                 cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAgg                                     >  2:1713599/1‑150 (MQ=255)
                                                                                                                 cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAgg                                     >  1:1762167/1‑150 (MQ=255)
                                                                                                                  gAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGACTGAGGTTGTCCAGCGGTGCATCGCTGTAGgc                                    >  2:2001659/1‑150 (MQ=255)
                                                                                                                       cggcgATCCCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCa                               <  1:2202119/150‑1 (MQ=255)
                                                                                                                        ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAg                              >  2:14446/1‑150 (MQ=255)
                                                                                                                          cgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGc                            <  2:303419/150‑1 (MQ=255)
                                                                                                                           gATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCAt                                              <  1:550575/131‑1 (MQ=255)
                                                                                                                           gATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCAt                                              >  2:550575/1‑131 (MQ=255)
                                                                                                                               aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcg                       >  2:877293/1‑150 (MQ=255)
                                                                                                                               aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcg                       >  1:1480935/1‑150 (MQ=255)
                                                                                                                                ccGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcgg                      >  2:791496/1‑150 (MQ=255)
                                                                                                                                      ggCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGt                <  1:2207429/150‑1 (MQ=255)
                                                                                                                                        cAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGCCTggcggcctacagc              <  1:694639/150‑8 (MQ=255)
                                                                                                                                            ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa          >  1:347873/2‑150 (MQ=255)
                                                                                                                                            ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCACGGAGGCCGGGTGCTCCGGCAAGCAGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa          >  2:457417/2‑150 (MQ=37)
                                                                                                                                            ccATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCCGAGGTGTCGATGGGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGGCCAGCGGGGCATCGCAGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCa          >  1:297607/2‑150 (MQ=37)
                                                                                                                                             cATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCAACGTCTTCCAg         >  2:16069/1‑150 (MQ=255)
                                                                                                                                              aTGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg                                                                                          >  1:1669675/1‑68 (MQ=255)
                                                                                                                                              aTGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg                                                                                          <  2:1669675/68‑1 (MQ=255)
                                                                                                                                               tGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTcagatc       >  1:1107178/1‑145 (MQ=255)
                                                                                                                                                 gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc     <  1:1713599/150‑1 (MQ=255)
                                                                                                                                                 gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc     <  2:1067338/150‑1 (MQ=255)
                                                                                                                                                   gAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCa   >  2:210456/1‑150 (MQ=255)
                                                                                                                                                    aaCTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGTTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAg  >  2:1934828/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAG  >  NC_002947/4824047‑4824344

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