Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,321,656 T→A 5.3% V249D (GTC→GAC)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,321,6560TA 5.3% 1211.2 / 35.0 350V249D (GTC→GAC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (121/209);  new base A (14/5);  total (136/214)
Fisher's exact test for biased strand distribution p-value = 2.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 354 total aligned reads displayed.

TCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGAA  >  NC_002947/4321509‑4321803
                                                                                                                                                    |                                                                                                                                                   
tCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                                                   <  2:35025/150‑1 (MQ=255)
  aCGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCAcc                                                                                                                                                 >  1:2204229/1‑150 (MQ=255)
    ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                               <  2:2777124/150‑1 (MQ=255)
      cccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGc                                                                                                                                             <  1:2449181/150‑1 (MQ=255)
       ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                                                            <  1:3288713/150‑1 (MQ=255)
       ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                                                            >  2:1418366/1‑150 (MQ=255)
       ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                                                            >  1:3266109/1‑150 (MQ=255)
           aaaaCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGGGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTCCCTGTCGGCCTGCCTCGCTGAGCCCGGTGTCCCCTTGGtcctc                                                                                                                                        >  2:3810530/1‑150 (MQ=255)
             aaCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagcga                                                                                                                                      <  2:3426868/150‑5 (MQ=255)
               cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATgg                                                                                                                                    <  2:2600832/150‑1 (MQ=255)
               cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATgg                                                                                                                                    <  2:4176632/150‑1 (MQ=255)
                ccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGt                                                                                                                                   <  2:6057819/150‑1 (MQ=255)
                 cAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTa                                                                                                                                  >  2:2216622/1‑150 (MQ=255)
                   ggCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGATGAGCATGGtgac                                                                                                                                <  1:2775350/150‑4 (MQ=255)
                     ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCa                                                                                                                              <  1:809271/150‑1 (MQ=255)
                     ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCa                                                                                                                              <  2:5861960/150‑1 (MQ=255)
                     ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCa                                                                                                                              <  2:3439700/150‑1 (MQ=255)
                     ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGATGAGCATGGtgacac                                                                                                                              <  1:2936140/150‑6 (MQ=255)
                      cTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAAGTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCTTGCTCATCATGGTACGCa                                                                                                                              <  2:1107814/150‑1 (MQ=255)
                         gTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCa                                                                                                                          >  1:6029544/1‑150 (MQ=255)
                             gACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTg                                                                                                                      >  2:5561233/1‑150 (MQ=255)
                             gACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTg                                                                                                                      <  1:4736553/150‑1 (MQ=255)
                              aCCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgc                                                                                                                     <  2:3926924/150‑1 (MQ=255)
                               ccccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgca                                                                                                                    <  2:5459084/150‑1 (MQ=255)
                                cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                                   >  2:810453/1‑150 (MQ=255)
                                cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                                   >  1:2546753/1‑150 (MQ=255)
                                cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCCCCGGCGCCGTGATCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCAGTCTCCCTGCAACGCTGAGCAAGGTGACACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                                   >  2:3167932/1‑150 (MQ=255)
                                 ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagc                                                                                                                  <  1:1300612/150‑1 (MQ=255)
                                  cTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagca                                                                                                                 >  2:1712633/1‑150 (MQ=255)
                                  cTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagca                                                                                                                 >  1:6394131/1‑150 (MQ=255)
                                  cTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagca                                                                                                                 >  1:5415831/1‑150 (MQ=255)
                                  cTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTAACTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGATCATAATGGTACACAGCCACCTgcagca                                                                                                                 >  1:6116990/1‑150 (MQ=255)
                                   tCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCAt                                                                                                                <  1:3989052/150‑1 (MQ=255)
                                       tCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCAtcttc                                                                                                            <  1:2823419/150‑1 (MQ=255)
                                         cAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTtctc                                                                                                          <  2:6799200/150‑1 (MQ=255)
                                          aggaggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCg                                                                                                         >  2:2092578/1‑150 (MQ=255)
                                                ccGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGcc                                                                                                   >  2:1967372/1‑150 (MQ=255)
                                                     ctgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGca                                                                                              >  2:3752147/1‑150 (MQ=255)
                                                       gcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaat                                                                                            <  1:3571509/150‑1 (MQ=255)
                                                       gcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaat                                                                                            >  2:3309745/1‑150 (MQ=255)
                                                       gcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaat                                                                                            >  2:5167260/1‑150 (MQ=255)
                                                           ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaa                                                                                         >  1:45549/1‑150 (MQ=255)
                                                           ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaa                                                                                         >  1:70209/1‑150 (MQ=255)
                                                            tttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                        <  1:5951388/150‑1 (MQ=255)
                                                            tttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                        <  2:6149372/150‑1 (MQ=255)
                                                            tttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                        <  1:1485919/150‑1 (MQ=255)
                                                               ttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                        >  2:4282805/1‑147 (MQ=255)
                                                               ttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                        <  1:4282805/147‑1 (MQ=255)
                                                               ttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcc                                                                                     <  1:4420961/150‑1 (MQ=255)
                                                                    cgccgGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCa                                                                                >  1:1851626/1‑150 (MQ=255)
                                                                        ggCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCAc                                                                            <  1:6664556/150‑1 (MQ=255)
                                                                              gtgtTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAgg                                                                      <  2:6288501/150‑1 (MQ=255)
                                                                                 ttCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTTTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCat                                                                   >  1:1680347/1‑150 (MQ=255)
                                                                                  tCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCata                                                                  >  1:3516196/1‑150 (MQ=255)
                                                                                   cGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCatat                                                                 >  1:5547610/1‑150 (MQ=255)
                                                                                       ccTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGtt                                                             <  1:2730055/150‑1 (MQ=255)
                                                                                             cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                                                       >  1:4983748/1‑150 (MQ=255)
                                                                                             cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                                                       >  1:6638024/1‑150 (MQ=255)
                                                                                             cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                                                       >  1:6545583/1‑150 (MQ=255)
                                                                                                       tCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                             >  2:1281792/1‑150 (MQ=255)
                                                                                                        cAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTg                                            <  1:4120198/150‑1 (MQ=255)
                                                                                                        cAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTg                                            <  2:17897/150‑1 (MQ=255)
                                                                                                         aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                           <  2:1154895/150‑1 (MQ=255)
                                                                                                         aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                           <  2:1046358/150‑1 (MQ=255)
                                                                                                         aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                           <  2:4965457/150‑1 (MQ=255)
                                                                                                          aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                          <  1:2200806/150‑1 (MQ=255)
                                                                                                          aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                          <  2:2470529/150‑1 (MQ=255)
                                                                                                                ccAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttg                                    >  2:915392/1‑150 (MQ=255)
                                                                                                                ccAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCACCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttg                                    >  1:6296849/1‑150 (MQ=255)
                                                                                                                 cAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttga                                   >  2:2308735/1‑150 (MQ=255)
                                                                                                                 cAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttga                                   >  1:3810482/1‑150 (MQ=255)
                                                                                                                 cAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttga                                   <  2:5955237/150‑1 (MQ=255)
                                                                                                                  ccttcccTGTCGGCCAGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgag                                  <  1:744317/145‑1 (MQ=255)
                                                                                                                  aaTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgag                                  >  2:2415184/1‑150 (MQ=255)
                                                                                                                   aTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgagt                                 <  2:6379893/150‑1 (MQ=255)
                                                                                                                   aTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgagt                                 <  1:2162652/150‑1 (MQ=255)
                                                                                                                    ttACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgagtt                                >  2:1505519/1‑150 (MQ=255)
                                                                                                                    ttACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgagtt                                <  1:2283401/150‑1 (MQ=255)
                                                                                                                       ccTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGAt                             <  2:6638024/150‑1 (MQ=255)
                                                                                                                         tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCa                           <  1:1461932/150‑1 (MQ=255)
                                                                                                                         tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCa                           <  1:5708838/150‑1 (MQ=255)
                                                                                                                           tCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt                         <  1:1034542/150‑1 (MQ=255)
                                                                                                                           tCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt                         <  1:1551957/150‑1 (MQ=255)
                                                                                                                           tCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt                         <  1:4860861/150‑1 (MQ=255)
                                                                                                                              gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTgcg                      >  2:6424794/1‑150 (MQ=255)
                                                                                                                                cTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCt                    <  1:3443366/150‑1 (MQ=255)
                                                                                                                                 tGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCtt                   <  1:4899873/150‑1 (MQ=255)
                                                                                                                                  gCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtc                  >  1:3684522/1‑150 (MQ=255)
                                                                                                                                  gCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtc                  <  1:1806751/150‑1 (MQ=255)
                                                                                                                                   cATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtca                 <  1:4477654/150‑1 (MQ=255)
                                                                                                                                      tgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcatca              >  1:2617958/4‑150 (MQ=255)
                                                                                                                                      tgatGAGCATGGTGACACCGTGCTCAGCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcatca              <  2:4845063/147‑1 (MQ=255)
                                                                                                                                         tGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCCCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgcg           >  2:4345217/1‑150 (MQ=255)
                                                                                                                                         tGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGc                                                                     >  2:4575664/1‑92 (MQ=255)
                                                                                                                                          gAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgcgc          >  1:3341128/1‑150 (MQ=255)
                                                                                                                                           aGCATGGTGGCACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCt         >  2:5373/1‑150 (MQ=255)
                                                                                                                                              aCGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGTCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGc      <  1:6796450/150‑1 (MQ=255)
                                                                                                                                              aCGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGc      >  1:4625441/1‑150 (MQ=255)
                                                                                                                                                 gtgtCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGa   >  1:3179164/1‑150 (MQ=255)
                                                                                                                                                  tgtCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGaa  <  2:2054465/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGAA  >  NC_002947/4321509‑4321803

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