Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,008,361 C→T 12.2% S259S (TCG→TCA PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,008,3610CT12.2% 212.0 / 17.1 82S259S (TCG→TCAPP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (6/4);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGA  >  NC_002947/3008218‑3008503
                                                                                                                                               |                                                                                                                                              
cTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCgg                                                                                                                                          >  1:330404/1‑150 (MQ=255)
 tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCTGCCCGGTCGTCGATGACCggg                                                                                                                                         >  1:1214869/1‑150 (MQ=255)
 tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCggg                                                                                                                                         >  1:1171263/1‑150 (MQ=255)
 tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCAGCCCGGTCGTCGATGACCggg                                                                                                                                         >  1:787385/1‑150 (MQ=255)
     cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTCCCATGGCTGCGCGCGCAAAGGGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                     >  1:357136/1‑150 (MQ=255)
      aGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATATTCCATGGCGGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGGCAGAAACCCGCCCGGTCATCGACGACCGGGCATTg                                                                                                                                    >  2:1249215/1‑150 (MQ=255)
        gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                  <  2:1171263/150‑1 (MQ=255)
        gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                  <  2:433668/150‑1 (MQ=255)
        gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                  <  1:1427970/150‑1 (MQ=255)
                ggCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCt                                                                                                                          <  1:784436/150‑1 (MQ=255)
                 gCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTg                                                                                                                         >  2:904127/1‑150 (MQ=255)
                       aCCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGGAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAc                                                                                                                   <  2:1327432/150‑1 (MQ=255)
                        ccGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGCGCTCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAcc                                                                                                                  <  1:1192213/150‑1 (MQ=255)
                           gTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTg                                                                                                               <  1:622175/150‑1 (MQ=255)
                            tCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGc                                                                                                              <  1:906939/150‑1 (MQ=255)
                            tCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATCGGCCTGCTCCTGCCCGACCTTGc                                                                                                              <  2:1972555/150‑1 (MQ=255)
                                cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                  <  1:2287419/126‑1 (MQ=255)
                                cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                  >  2:2287419/1‑126 (MQ=255)
                                cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAg                                                                                                          >  1:2113325/1‑150 (MQ=255)
                                 ttGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGc                                                                                                         >  1:1992156/1‑150 (MQ=255)
                                  tGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGcc                                                                                                        <  2:2171918/150‑1 (MQ=255)
                                    aaaTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTgt                                                                                                                                   >  1:567495/1‑120 (MQ=255)
                                    aaaTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTgg                                                                                                                                   <  2:567495/121‑1 (MQ=255)
                                    aaaTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgg                                                                                                      >  2:2196495/1‑150 (MQ=255)
                                      aTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgggg                                                                                                    <  2:1273476/150‑1 (MQ=255)
                                         ttCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGAc                                                                                                 <  2:1358960/150‑1 (MQ=255)
                                         ttCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGATATAACCCTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGAc                                                                                                 <  1:869591/150‑1 (MQ=255)
                                             aTGGCTGCGGTGGCAAAGTGCGCTTCGCGGGTGGAGTGCGAAAACAGTGTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCGTCGACAGCCGGGCATTGGCCTGCTGCTGCCCGACCGGGCCCAGCCGGGGTACGaaa                                                                                             >  2:678634/1‑150 (MQ=255)
                                             aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             <  2:1090937/150‑1 (MQ=255)
                                             aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             >  1:1427900/1‑150 (MQ=255)
                                             aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             <  2:2113325/150‑1 (MQ=255)
                                                gCTGCGCGCGCAAAGTGCGCTTCGCGGTCGTAGTGTGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATAc                                                                                          <  2:482934/150‑1 (MQ=255)
                                                       gcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATgg                                                                                   >  2:1800901/1‑150 (MQ=255)
                                                         gcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGcc                                                                                 >  1:200967/1‑150 (MQ=255)
                                                             aGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGc                                                                             <  1:904127/150‑1 (MQ=255)
                                                             aGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGc                                                                             <  2:745811/150‑1 (MQ=255)
                                                             aGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAACACAGCATGGCCAAGc                                                                             <  2:357136/150‑1 (MQ=255)
                                                                 cgcTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         <  1:2196495/150‑1 (MQ=255)
                                                                         ggTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         <  1:1109261/142‑1 (MQ=255)
                                                                         ggTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         >  2:1109261/1‑142 (MQ=255)
                                                                                 tGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCtt                                                                                                                <  1:2020425/95‑1 (MQ=255)
                                                                                 tGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCtt                                                                                                                >  2:2020425/1‑95 (MQ=255)
                                                                                    gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGcc                                                                            >  1:1980634/1‑128 (MQ=255)
                                                                                    gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGcc                                                                            <  2:1980634/128‑1 (MQ=255)
                                                                                     aaaaCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACAGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCCCGGCCAAGCCGCGCATCACGACCCAGGCAATGa                                                     >  1:121177/1‑150 (MQ=255)
                                                                                           gTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             <  1:1127092/104‑1 (MQ=255)
                                                                                           gTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             >  2:1127092/1‑104 (MQ=255)
                                                                                            ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAATTACAGCATgg                                                                                   >  1:576897/1‑113 (MQ=255)
                                                                                            ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAATTACAGCATgg                                                                                   <  2:576897/113‑1 (MQ=255)
                                                                                                aCCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATc                                          >  1:929431/1‑150 (MQ=255)
                                                                                                    tggtgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGACATCAAGCATCCTCAATCCTAACATCACTACCCAGGCAATGAAGGGCCGCATCAGtt                                      >  2:2135438/1‑150 (MQ=255)
                                                                                                    tggtgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGtt                                      >  1:539510/1‑150 (MQ=255)
                                                                                                    tggtgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGtt                                      >  2:1652060/1‑150 (MQ=255)
                                                                                                        gATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             <  2:1936455/91‑1 (MQ=255)
                                                                                                        gATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                             >  1:1936455/1‑91 (MQ=255)
                                                                                                         aTGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                 <  2:1427900/150‑1 (MQ=255)
                                                                                                          tGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCTACGACCGGTTATTGTCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgat                                <  1:1249215/150‑3 (MQ=255)
                                                                                                           gTTGGTGTAAGCGTCAGAAAACCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGc                               >  2:43811/1‑150 (MQ=255)
                                                                                                              ggTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCaaa                            >  2:369826/1‑150 (MQ=255)
                                                                                                              ggTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCaaa                            >  1:998593/1‑150 (MQ=255)
                                                                                                                     gCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCGATGAAGGGCAGCATCAGTTTCAGGGCAAAGccgccg                     >  1:1906342/1‑150 (MQ=255)
                                                                                                                           gAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAgg               <  2:929431/150‑1 (MQ=255)
                                                                                                                               cccgcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCaa           >  2:2067223/1‑150 (MQ=255)
                                                                                                                               accgcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTCCCCGCCATTGCCCAGCCGGGGCACGAAATACAGCCTGGCCATGCCGCACCGGAAAACTGATGCAAAGAACCACAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCaa           <  1:43811/149‑1 (MQ=255)
                                                                                                                               aatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         <  2:1913111/85‑1 (MQ=39)
                                                                                                                               aatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         <  2:751596/85‑1 (MQ=38)
                                                                                                                               aatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         >  1:751596/4‑88 (MQ=38)
                                                                                                                               aatgcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                         >  1:1913111/4‑88 (MQ=39)
                                                                                                                                ccgcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGtt                                      <  1:2142003/122‑1 (MQ=255)
                                                                                                                                ccgcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGtt                                      >  2:2142003/1‑122 (MQ=255)
                                                                                                                                ccgcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAAt          >  1:228650/1‑150 (MQ=255)
                                                                                                                                  gcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGAGGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGCCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGa        >  2:251171/1‑150 (MQ=255)
                                                                                                                                  gcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCAGCCGGCCAGGGCAATGa        >  2:261768/1‑150 (MQ=255)
                                                                                                                                   cccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAAt          <  2:1232057/147‑1 (MQ=255)
                                                                                                                                   cccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAAt          >  1:1232057/1‑147 (MQ=255)
                                                                                                                                   cccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGAt       >  2:1953893/1‑150 (MQ=255)
                                                                                                                                    ccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATc      <  1:1953893/150‑1 (MQ=255)
                                                                                                                                     cgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATcg     <  1:2067223/150‑1 (MQ=255)
                                                                                                                                      ggTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGCAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATcgc    <  2:330404/150‑1 (MQ=255)
                                                                                                                                      ggTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATcgc    >  1:1662464/1‑150 (MQ=255)
                                                                                                                                        tCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGa  <  1:1800901/150‑1 (MQ=255)
                                                                                                                                        tCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGa  <  1:1652060/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGA  >  NC_002947/3008218‑3008503

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: