Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,824,197 C→T 6.1% E2346E (GAG→GAA pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,824,1970CT 6.1% 355.8 / 6.7 123E2346E (GAG→GAApvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (47/67);  new base T (5/3);  total (53/70)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 128 total aligned reads displayed.

TGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCT  >  NC_002947/4824048‑4824346
                                                                                                                                                     |                                                                                                                                                     
tGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAc                                                                                                                                                       >  2:1614196/1‑150 (MQ=255)
tGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAc                                                                                                                                                       >  1:1408356/1‑150 (MQ=255)
    cgcgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                   >  1:69486/1‑150 (MQ=255)
       gcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCTGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCAt                                                                                                                                                <  2:1601126/150‑1 (MQ=255)
       gcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCAt                                                                                                                                                <  2:1601094/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)
                       accacACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                   >  1:41568/1‑131 (MQ=255)
                       accacACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                   <  2:41568/131‑1 (MQ=255)
                         cacacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCt                                                                                                                              >  2:394077/1‑150 (MQ=255)
                          acacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCtt                                                                                                                             >  2:385077/1‑150 (MQ=255)
                           cacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTg                                                                                                                            <  2:1670099/150‑1 (MQ=255)
                            accacaGCGCCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTgg                                                                                                                           <  1:2099612/147‑1 (MQ=255)
                            acaacaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTgg                                                                                                                           >  1:2428259/1‑150 (MQ=255)
                              aacaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGt                                                                                                                         >  1:2209175/1‑150 (MQ=255)
                               acaGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGtt                                                                                                                        >  1:301924/1‑150 (MQ=255)
                                caGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGttt                                                                                                                       <  1:690635/150‑1 (MQ=255)
                                  gCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGc                                                                                                                     <  2:2209175/150‑1 (MQ=255)
                                   cGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTNGcc                                                                                                                    <  2:1520145/150‑1 (MQ=255)
                                       aCCCACAGGCGCTCGCCGGCACCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACGCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTg                                                                                                                <  2:778034/150‑1 (MQ=255)
                                        cccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTgg                                                                                                               >  2:2210791/1‑150 (MQ=255)
                                         ccACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCa                                                                                                                                                 <  2:2015211/115‑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)
                                           acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                            <  1:2080243/150‑1 (MQ=255)
                                           acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                            >  1:2279325/1‑150 (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)
                                                   tcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGGGCCTGAGgtg                                                                                                    <  2:1553335/150‑1 (MQ=255)
                                                    cgctGGCAGCGTCGATCTTGAACGAATAGAAATTCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                   <  2:1097456/150‑1 (MQ=255)
                                                         gCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATc                                                                                              >  2:43003/1‑150 (MQ=255)
                                                          cAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                             <  2:776510/150‑1 (MQ=255)
                                                          cAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                             <  2:39188/150‑1 (MQ=255)
                                                            gCGTCGAAATTGATAGAATACAAATGCAGCTCGCAATCAGCGCTGCGCATGCCGAATGCGGCGATCAGCGCCTGGCAGTGCAGGGCGAACCCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGCAGCCTGAGGTGTAGATCAgc                                                                                           <  2:356638/150‑2 (MQ=255)
                                                               tCGAAATTGATCGAATAGAAATGCAGCTCGCAATCAGCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTaa                                                                                        <  1:1081296/150‑1 (MQ=255)
                                                                 gAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGc                                                                                      <  1:1095785/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)
                                                                    aTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAc                                                                                                                                                       <  2:1284448/82‑1 (MQ=255)
                                                                        aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCg                                                                                                                   <  1:2189525/114‑1 (MQ=255)
                                                                        aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCg                                                                                                                   >  2:2189525/1‑114 (MQ=255)
                                                                        aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgc                                                                               >  2:837736/1‑150 (MQ=255)
                                                                         tCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgct                                                                              <  1:2210791/150‑1 (MQ=255)
                                                                            aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                       >  2:1823709/1‑122 (MQ=255)
                                                                            aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                       <  1:1823709/122‑1 (MQ=255)
                                                                                gAAATGCAGCTCGCAATCATCGCTGCGCATGCCGACTGTGGCGCTCACCGCCTGGCAGTGCATGGCAGACGGGCCATGGTTGTGCATCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGGAGCTCAGGTAAGCCTGGTGCTGCGGCAg                                                                       >  1:1023842/1‑150 (MQ=255)
                                                                                   aTGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCt                                                                    >  1:1276168/1‑150 (MQ=255)
                                                                                      cAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                                 <  2:69486/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                                                                                   <  1:1591940/128‑1 (MQ=255)
                                                                                          tCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTg                                                             >  1:672742/1‑150 (MQ=255)
                                                                                          tCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTg                                                             >  2:2113710/1‑150 (MQ=255)
                                                                                           cGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGt                                                            <  1:1135802/150‑1 (MQ=255)
                                                                                            gCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTc                                                           <  2:672742/150‑1 (MQ=255)
                                                                                             cAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgctgc                                                                            >  1:1019233/1‑132 (MQ=255)
                                                                                              aaTCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAAGTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                             >  2:1283508/1‑114 (MQ=255)
                                                                                                     gCTGCGCATGCCGAATGCGGCGATCAGCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGc                                                  <  1:1064695/150‑1 (MQ=255)
                                                                                                         cgcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCg                                              >  2:694639/1‑150 (MQ=255)
                                                                                                          gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGc                                             <  2:76054/150‑1 (MQ=255)
                                                                                                             tGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGt                                                                                 >  2:900858/1‑111 (MQ=255)
                                                                                                             tGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGt                                                                                 <  1:900858/111‑1 (MQ=255)
                                                                                                              gCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGTGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTa                                         <  2:1276168/150‑1 (MQ=255)
                                                                                                              gCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTa                                         <  2:2398191/150‑1 (MQ=255)
                                                                                                               ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTg                                                                                   <  1:2401603/107‑1 (MQ=255)
                                                                                                                cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                   >  1:810262/1‑90 (MQ=255)
                                                                                                                cGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAgg                                       >  2:1713599/1‑150 (MQ=255)
                                                                                                                 gAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGACTGAGGTTGTCCAGCGGTGCATCGCTGTAGgc                                      >  2:2001659/1‑150 (MQ=255)
                                                                                                                  aaTGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGgcc                                     <  2:2428259/150‑1 (MQ=255)
                                                                                                                      cggcgATCCCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCa                                 <  1:2202119/150‑1 (MQ=255)
                                                                                                                       ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGGCCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAg                                >  1:1959990/1‑150 (MQ=255)
                                                                                                                       ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAg                                >  2:14446/1‑150 (MQ=255)
                                                                                                                       ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAg                                <  1:2400518/150‑1 (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)
                                                                                                                                           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)
                                                                                                                                             aTGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTgg                           <  2:1196248/133‑1 (MQ=255)
                                                                                                                                             aTGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTgg                           >  1:1196248/1‑133 (MQ=255)
                                                                                                                                                gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc       <  1:1713599/150‑1 (MQ=255)
                                                                                                                                                gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTc       <  2:1067338/150‑1 (MQ=255)
                                                                                                                                                  gAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGTCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCa     >  1:1761929/1‑150 (MQ=255)
                                                                                                                                                  gAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCa     >  2:210456/1‑150 (MQ=255)
                                                                                                                                                   aaCTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGTTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAg    >  2:1934828/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCGCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAGCT  >  NC_002947/4824048‑4824346

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: