Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,414,515 T→C 5.5% S218P (TCC→CCC)  acrE → multidrug efflux pump membrane fusion lipoprotein AcrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,414,5150TC 5.5% 222.0 / 2.5 72S218P (TCC→CCC) acrEmultidrug efflux pump membrane fusion lipoprotein AcrE
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base C (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AATCTTGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGAT  >  NC_000913/3414368‑3414664
                                                                                                                                                   |                                                                                                                                                     
aaTCTTGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATcc                                                                                                                                                     >  2:823184/1‑150 (MQ=255)
aaTCTTGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATcc                                                                                                                                                     <  2:1177904/150‑1 (MQ=255)
    ttGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCa                                                                                                                                                 <  1:1001919/150‑1 (MQ=255)
    ttGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCa                                                                                                                                                 <  2:720663/150‑1 (MQ=255)
      gCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAAc                                                                                                                                               >  2:876971/1‑150 (MQ=255)
       cTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACg                                                                                                                                              >  2:1105019/1‑150 (MQ=255)
           taCCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGAttt                                                                                                                                          >  2:663774/1‑150 (MQ=255)
             ccAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGAGCCAATCCAGGAACg                                                                                                                                              <  2:1057648/144‑1 (MQ=255)
             ccAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACg                                                                                                                                              >  1:1057648/1‑144 (MQ=255)
              cAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAACCCAGCAACGATTTTAt                                                                                                                                       >  2:1044524/1‑150 (MQ=255)
              cAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAACCCAGCAACGATTTTAt                                                                                                                                       >  2:1039496/1‑150 (MQ=255)
                 aGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAg                                                                                                                                    <  2:539512/150‑1 (MQ=255)
                      cTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGCCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGa                                                                                                                               >  2:114273/1‑150 (MQ=255)
                        gcgcCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAg                                                                                                                             <  1:896778/150‑1 (MQ=255)
                          gcCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGAAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCa                                                                                                                           >  2:1189399/1‑150 (MQ=255)
                            cAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATcc                                                                                                                                                     >  2:952784/1‑122 (MQ=255)
                            cAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATcc                                                                                                                                                     <  1:952784/122‑1 (MQ=255)
                               ttAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCg                                                                                                                      <  2:519745/150‑1 (MQ=255)
                                 aGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTa                                                                                                                    >  2:362876/1‑150 (MQ=255)
                                 aGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTa                                                                                                                    >  2:919774/1‑150 (MQ=255)
                                   cGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTaga                                                                                                                  >  2:1119759/1‑150 (MQ=255)
                                    ggACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTagag                                                                                                                 <  1:823184/150‑1 (MQ=255)
                                    ggACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTagag                                                                                                                 <  1:876971/150‑1 (MQ=255)
                                      aCGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCa                                                                                                               <  2:599110/150‑1 (MQ=255)
                                            cGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGaaa                                                                                                         <  1:492563/150‑1 (MQ=255)
                                            cGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGaaa                                                                                                         <  1:497347/150‑1 (MQ=255)
                                                aaaTCGACTGTGACCGAAGGCGCTCTTATCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTg                                                                                                     <  1:377215/150‑1 (MQ=255)
                                                 aaTCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGc                                                                                                    >  1:1056527/1‑150 (MQ=255)
                                                 aaTCGACTGTGACCGAAGGCGCTCTTGTCACTAAAGGGCAAACGACTGAACTCGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGAGACCCAATACAGCAACGATTTTATGAGGCTGAAGAAATCCGTAGAGCAAGGAAATTAGc                                                                                                    >  1:567787/1‑150 (MQ=255)
                                                  aTCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCa                                                                                                   <  1:1418370/150‑1 (MQ=255)
                                                   tCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCAt                                                                                                  <  2:611164/150‑1 (MQ=255)
                                                       cTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAgg                                                                                              >  2:968767/1‑150 (MQ=255)
                                                          tgACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGaaa                                                                                           <  1:1105019/150‑1 (MQ=255)
                                                           gACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATAATATGAGCCTGAAGCAATACGTAGAGCAAGGAAATTTGAATAAGCaaaa                                                                                          >  2:634529/1‑150 (MQ=39)
                                                              cGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGc                                                                                       >  2:675948/1‑150 (MQ=255)
                                                              cGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGc                                                                                       >  1:728763/1‑150 (MQ=255)
                                                                    cgcTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAg                                                                                 >  1:666386/1‑150 (MQ=255)
                                                                        ctTGTCACTAATGGGCAAACGACTGAACTGGCCACTGTCCAGCAGCACCATCCTAACAACGTTGATCTGCCCCAATCCAGCAACGATATTATGAGCCTCAAGCAATCCGTACAGCAAGGAAAAATGCAAAACGAAAACGCCACCAGCAAc                                                                             >  2:1258541/1‑150 (MQ=255)
                                                                           gTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTa                                                                          <  2:1105172/150‑1 (MQ=255)
                                                                            tCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTag                                                                         <  1:481084/150‑1 (MQ=255)
                                                                                 aaTGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACg                                                                            <  2:12721/142‑1 (MQ=255)
                                                                                 aaTGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACg                                                                            >  1:12721/1‑142 (MQ=255)
                                                                                      gCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCAt                                                               >  2:526710/1‑150 (MQ=255)
                                                                                        aaaCGACTTAACTGGCGAATGTCCAGCAGCTCGATCCTATGTACTTTTATGTGACCCAATCCATCATCGATTTTATGAGGCTGAAGCAATCCGTAGAGGAAGGAAATTTTCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATgg                                                             <  2:567787/150‑1 (MQ=255)
                                                                                         aaCGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGa                                                            >  1:734813/1‑150 (MQ=255)
                                                                                         aaCGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGa                                                            <  2:1056527/150‑1 (MQ=255)
                                                                                                 aaCTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTc                                                    <  1:968767/150‑1 (MQ=255)
                                                                                                     ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTagag                                                                       <  1:1412563/127‑1 (MQ=255)
                                                                                                     ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTagag                                                                       >  2:1412563/1‑127 (MQ=255)
                                                                                                     ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAAc                                                >  1:1177156/1‑150 (MQ=255)
                                                                                                       cGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAAGGTAGAGTTGGTCATGGAAAACGGTCAAACCt                                              >  2:692646/1‑150 (MQ=255)
                                                                                                         aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAACCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTAt                                            <  1:1039496/150‑1 (MQ=255)
                                                                                                         aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAACCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTAt                                            <  1:1044524/150‑1 (MQ=255)
                                                                                                          cTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATc                                           <  2:989153/150‑1 (MQ=255)
                                                                                                             tCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATcccc                                        >  2:1296871/1‑150 (MQ=255)
                                                                                                              ccAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCt                                       >  1:459133/1‑150 (MQ=255)
                                                                                                                  cagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAACTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGaaa                                   <  2:1203255/150‑1 (MQ=255)
                                                                                                                  cagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGaaa                                   <  1:1262012/150‑1 (MQ=255)
                                                                                                                       cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTAc                              <  2:828959/150‑1 (MQ=255)
                                                                                                                            cTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGc                         <  2:197306/150‑1 (MQ=255)
                                                                                                                             tATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGt                                <  1:288637/142‑1 (MQ=255)
                                                                                                                             tATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGt                                >  2:288637/1‑142 (MQ=255)
                                                                                                                               tCTACGTTGATGTGACCCAATCCAGCAACTATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAt                      <  1:1142419/150‑1 (MQ=255)
                                                                                                                               tCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCATCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAt                      >  1:923570/1‑150 (MQ=255)
                                                                                                                                  aCGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtct                   >  2:653371/1‑150 (MQ=255)
                                                                                                                                    gTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                 <  1:227383/150‑1 (MQ=255)
                                                                                                                                    gTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                 <  1:902065/150‑1 (MQ=255)
                                                                                                                                        aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg             <  1:362876/150‑1 (MQ=255)
                                                                                                                                        aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg             <  1:919774/150‑1 (MQ=255)
                                                                                                                                        aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg             <  1:1258541/150‑1 (MQ=255)
                                                                                                                                         tgtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgt            <  1:1119759/150‑1 (MQ=255)
                                                                                                                                           tgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGa          >  2:571034/1‑150 (MQ=255)
                                                                                                                                                   tCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGATCAAGGAAATTTGCATAAGGAAAACGCCCCCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACTCTGCAATTCTCCGATGTGACCGTtgat  <  2:591264/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATCTTGCTTATACCAAAGTCACTGCGCCAATTAGCGGACGTATCGGCAAATCGACTGTGACCGAAGGCGCTCTTGTCACTAATGGGCAAACGACTGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGAT  >  NC_000913/3414368‑3414664

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