Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,826,672 G→A 6.2% L1521L (CTC→CTT pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,826,6720GA 6.2% 420.7 / 10.7 140L1521L (CTC→CTTpvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (73/58);  new base A (4/5);  total (77/63)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACGGGCCAGG  >  NC_002947/4826526‑4826819
                                                                                                                                                  |                                                                                                                                                   
atctcACTGTAATAGCCGCCCAGTAGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGACCGTTCTCCGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGCTCGAgg                                                                                                                                                  <  1:1574276/146‑1 (MQ=255)
   ccACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               <  1:823952/150‑1 (MQ=255)
   ccACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCAACCAAGCAGCCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               >  1:1694228/1‑150 (MQ=255)
        gTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAg                                                                                                                                          <  1:1747037/150‑1 (MQ=255)
           aTAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGCCCTGAACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTa                                                                                                                                       >  2:661950/1‑150 (MQ=255)
                cGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCa                                                                                                                                  >  1:620960/1‑150 (MQ=255)
                 gACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAg                                                                                                                                 <  2:1035704/150‑1 (MQ=255)
                 gACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAg                                                                                                                                 >  2:1150091/1‑150 (MQ=255)
                  aCCAGTTGCGCACCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAgg                                                                                                                                <  1:1063061/150‑1 (MQ=255)
                           gcgcCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTCGACCAGCAGGTCAAggt                                                                                                                                                 <  2:258153/124‑2 (MQ=255)
                              ccAATGGCGTCCTTGTGCATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGa                                                                                                                    >  2:1190301/1‑150 (MQ=255)
                               cAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               >  1:620217/1‑122 (MQ=255)
                               cAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               <  2:620217/122‑1 (MQ=255)
                               cAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGaa                                                                                                                   >  1:2312325/1‑150 (MQ=255)
                                aaTGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAAc                                                                                                                  <  1:1150091/150‑1 (MQ=255)
                                         cTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCg                                                                                                         >  2:1810939/1‑150 (MQ=255)
                                          ttGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCgg                                                                                                        >  2:1698415/1‑150 (MQ=255)
                                                aTCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCGAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGa                                                                                                  <  2:1694228/150‑1 (MQ=255)
                                                aTCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGa                                                                                                  <  2:183911/150‑1 (MQ=255)
                                                 tCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGAt                                                                                                 >  1:1825431/1‑150 (MQ=255)
                                                  cAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGgcg                                                                                                     >  2:1940450/1‑145 (MQ=255)
                                                  cAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGgcg                                                                                                     <  1:1940450/145‑1 (MQ=255)
                                                  cAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATc                                                                                                <  2:1900687/150‑1 (MQ=255)
                                                    gCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAg                                                                                                                                          <  1:1690581/106‑1 (MQ=255)
                                                    gCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAg                                                                                                                                          >  2:1690581/1‑106 (MQ=255)
                                                    gCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGAgcgc                                                                                              >  1:347129/1‑150 (MQ=255)
                                                    gCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCCGCCGAACCTAGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGCCGTCGAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGc                                                                                              >  1:2176137/1‑150 (MQ=255)
                                                      aCCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCa                                                                                            >  1:1031792/1‑150 (MQ=255)
                                                      aCCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCa                                                                                            <  2:1747064/150‑1 (MQ=255)
                                                      aCCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCa                                                                                            >  1:1060727/1‑150 (MQ=255)
                                                        cAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACt                                                                                          >  1:2048944/1‑150 (MQ=255)
                                                        cAGCGCCTGGTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACt                                                                                          >  2:2144429/1‑150 (MQ=255)
                                                         aGCGCCTGTTCCACACCCGCCTGTTCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTACTGGTCGAGGCGGCCAAGGTATTTCAGGGCCTCCAGCTGAACCTGCCAGAGGGCGCGATCGCAACTg                                                                                         >  2:1070379/1‑150 (MQ=255)
                                                         aGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTCGACCAGCAGGTCAAggt                                                                                                                                                 >  1:1011397/1‑93 (MQ=255)
                                                         aGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTCGACCAGCAGGTCAAggt                                                                                                                                                 <  2:1011397/94‑2 (MQ=255)
                                                          gcgcCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGc                                                                                        >  2:377339/1‑150 (MQ=255)
                                                          gcgcCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGc                                                                                        >  2:1524615/1‑150 (MQ=255)
                                                             ccTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGt                                                                                     <  1:1689226/150‑1 (MQ=255)
                                                                gTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACa                                                                                  >  1:1633729/1‑150 (MQ=255)
                                                                   ccACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCACGAGCTATGAGGGCCTCCAGCTGAACCTGCCGTCGTGCGCGGTCGCCACTGCGGTACAGAc                                                                               >  1:1778878/1‑150 (MQ=255)
                                                                   ccACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGAc                                                                               <  2:1060727/150‑1 (MQ=255)
                                                                     acacCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGc                                                                             >  1:378308/1‑150 (MQ=255)
                                                                         ccGCCTGTGCCAGCAGGCAGGGCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGAc                                                                         <  1:1360080/150‑1 (MQ=255)
                                                                         ccGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGAc                                                                         >  2:1242070/1‑150 (MQ=255)
                                                                         ccGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCCCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGAc                                                                         <  2:1288911/150‑1 (MQ=255)
                                                                         ccGCCTCTGCCAGCAAAAAAGCCTGCACATCCTCGAGCTCGACCCGGCAACCGCGCAACTTGACCTGCTGGTCGAGACGGCAAAGGTATTCCAGGGCCTACAGCTGAACCTGACAGCGGACACGATCGCCACTGCGGTACAGAcaatgaa                                                                         >  2:482701/1‑144 (MQ=255)
                                                                             cTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCAtt                                                                     <  2:1839734/150‑1 (MQ=255)
                                                                              tgtgCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               >  2:2187322/1‑75 (MQ=255)
                                                                              tgtgCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAggcgg                                                                                                                                               <  1:2187322/75‑1 (MQ=255)
                                                                              tgtgCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTg                                                                    >  2:271927/1‑150 (MQ=255)
                                                                               gtgCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGc                                                                   <  2:1031792/150‑1 (MQ=255)
                                                                                     gcaggcagGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTa                                                                                    >  2:1644780/1‑127 (MQ=255)
                                                                                     gcaggcagGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTa                                                                                    <  1:1644780/127‑1 (MQ=255)
                                                                                         gcaggcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAg                                                                                                                                          <  1:910314/69‑1 (MQ=255)
                                                                                         gcaggcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAg                                                                                                                                          >  2:973284/1‑69 (MQ=255)
                                                                                         gcaggcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTg                                                         >  1:1471606/1‑150 (MQ=255)
                                                                                           aggcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTggg                                                       <  1:1475311/150‑1 (MQ=255)
                                                                                           aggcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTggg                                                       <  2:887764/150‑1 (MQ=255)
                                                                                             gcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGcc                                                                                                                             <  2:1136896/78‑1 (MQ=255)
                                                                                             gcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGcc                                                                                                                             >  1:1136896/1‑78 (MQ=255)
                                                                                             gcCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGTCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGCCGATGACCATTGCCGTCGGCCTCGGGt                                                     >  2:449709/1‑150 (MQ=255)
                                                                                                tGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACg                                                  >  1:1980281/1‑150 (MQ=255)
                                                                                                 gCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTgg                                                        >  1:1391662/1‑143 (MQ=255)
                                                                                                 gCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTgg                                                        <  2:1391662/143‑1 (MQ=255)
                                                                                                  cACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTc                                                                                                                           <  2:739655/75‑1 (MQ=255)
                                                                                                   aCTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCt                                                                                                                      <  1:426438/79‑1 (MQ=255)
                                                                                                   aCTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAg                                               >  2:723483/1‑150 (MQ=255)
                                                                                                    cTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAgc                                              <  1:1764582/150‑1 (MQ=255)
                                                                                                        cTCGGGCTCTACGCGGAAACCGCGCACCTTGACCCGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTc                                          <  2:1588088/150‑1 (MQ=255)
                                                                                                         tCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCg                                         >  1:767442/1‑150 (MQ=255)
                                                                                                          cGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCgg                                        <  1:1038791/150‑1 (MQ=255)
                                                                                                            ggCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCa                                      <  2:887488/150‑1 (MQ=255)
                                                                                                                cGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGa                                  <  2:1980281/150‑1 (MQ=255)
                                                                                                                    gcgGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCg                                                              >  2:951789/1‑118 (MQ=255)
                                                                                                                    gcgGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCg                                                              <  1:951789/118‑1 (MQ=255)
                                                                                                                        aaaCCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCgggg                          <  2:1471606/150‑1 (MQ=255)
                                                                                                                          aCCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCg                                                              <  1:1493018/112‑1 (MQ=255)
                                                                                                                          aCCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCg                                                              >  2:1493018/1‑112 (MQ=255)
                                                                                                                          aCCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGggcg                        <  2:620960/150‑1 (MQ=255)
                                                                                                                            cgcgCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGc                                                                             >  1:720058/1‑95 (MQ=255)
                                                                                                                            cgcgCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGTCCTGTGGTACGAAGCGCTCGGCAGTGAGGCCGGggcggc                      <  1:449709/150‑1 (MQ=255)
                                                                                                                             gcgcACCTTGACCTGCTGGTCGAGGCGACCAAGGTATTCAAGGGCCTCCAGCTGAACCTGCCAGCGGGCGAGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGGCTGGGGTACGAAGCGCTCGGCAGAGAGGGGGGgggggcg                     >  1:1129486/1‑149 (MQ=255)
                                                                                                                                cACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCaa                  <  1:723483/150‑1 (MQ=255)
                                                                                                                                   cTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTa               <  1:1518104/150‑1 (MQ=255)
                                                                                                                                   cTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTa               >  1:812639/1‑150 (MQ=255)
                                                                                                                                   cTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTa               >  2:780823/1‑150 (MQ=255)
                                                                                                                                   cTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTa               >  2:1411217/1‑150 (MQ=255)
                                                                                                                                      gACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAAcc            >  1:49137/1‑150 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGa                                                                          >  2:1358613/1‑87 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGa                                                                          <  1:1358613/87‑1 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCa           <  2:1488914/150‑1 (MQ=255)
                                                                                                                                         ctgctgGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTg                                                                    <  1:603754/91‑1 (MQ=255)
                                                                                                                                         ctgctgGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACg         <  1:21930/150‑1 (MQ=255)
                                                                                                                                          tgctgGGCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACgg        >  2:2144228/1‑150 (MQ=255)
                                                                                                                                           gcaggTCAAGGTGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAg                                 <  2:1250091/121‑1 (MQ=255)
                                                                                                                                           gcaggTCAAGGTGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAg                                 >  1:1250091/4‑124 (MQ=255)
                                                                                                                                            ctgGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGTACCCGGGGCGGCCCAAGTAACCACGGGc      <  1:1230357/150‑1 (MQ=255)
                                                                                                                                               gTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACGGGCCAg   >  1:2258607/1‑150 (MQ=255)
                                                                                                                                               gTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACGGGCCAg   <  2:2312325/150‑1 (MQ=255)
                                                                                                                                                tCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACGGGCCAgg  >  2:2682/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCCACTGTAATAGCCGACCAGTTGCGCGCCAATGGCGTCCTTGTGTATCAGCACCAGCGCCTGTTCCACACCCGCCTGTGCCAGCAGGCAGGCCTGCACTTCCTCGGGCTCGACGCGGAAACCGCGCACCTTGACCTGCTGGTCGAGGCGGCCAAGGTATTCCAGGGCCTCCAGCTGAACCTGCCAGCGGGCGCGATCGCCACTGCGGTACAGACGCTGACCATTGCCGTCGGCCTGGGGTACGAAGCGCTCGGCAGTGAGGCCGGGGCGGCCCAAGTAACCACGGGCCAGG  >  NC_002947/4826526‑4826819

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: