Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,873 T→A 18.1% intergenic (+90/+62) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8730TA18.1% 421.7 / 31.2 177intergenic (+90/+62)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (33/104);  new base A (27/2);  total (70/107)
Fisher's exact test for biased strand distribution p-value = 2.47e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022
                                                                                                                                                     |                                                                                                                                                     
tcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGACTTTTACATCCCGTACGTTCCCCTCACCCAAACCct                                                                                                                                                       >  2:3275876/1‑150 (MQ=255)
tcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCct                                                                                                                                                       <  1:2585644/150‑1 (MQ=255)
 ccGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTTGCCTGACCCTAACCcta                                                                                                                                                      >  1:2037428/1‑149 (MQ=255)
 ccGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTTACCCTAACCctc                                                                                                                                                      >  2:3531916/1‑150 (MQ=255)
 ccGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCctc                                                                                                                                                      >  1:609609/1‑150 (MQ=255)
    aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcc                                                                                                                                                   <  1:325968/150‑1 (MQ=255)
      ctctTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcccc                                                                                                                                                 <  1:2765591/150‑1 (MQ=255)
         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCACAACCCaaa                                                                                                                                              >  2:1447008/1‑150 (MQ=39)
            tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                                           <  1:2470252/150‑1 (MQ=255)
            tGCCGCGCCAATCGCACCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                                           <  1:1072575/150‑1 (MQ=255)
             gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                                          <  1:3053603/150‑1 (MQ=255)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCTCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       <  2:2112388/150‑1 (MQ=255)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGa                                                                                                                                       >  1:2859571/1‑150 (MQ=255)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCCCACCCCAAAGGGGCgg                                                                                                                                       >  2:1682296/1‑149 (MQ=255)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAAACCACACCCCCAACGGGCgc                                                                                                                                       >  1:2273866/1‑149 (MQ=38)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAAACCACACCCCAAAGGGGCGa                                                                                                                                       >  1:364535/1‑150 (MQ=39)
                   ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCTCACCCCCAAGGGGCGCggg                                                                                                                                    >  2:3820434/1‑150 (MQ=255)
                    cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGGAGGCCGGATAATATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTACATCCCGTACGTTCCCCGCCCCCACACCCACTCCCCAAACGGGCGAgggg                                                                                                                                   >  2:2297316/1‑150 (MQ=255)
                       tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAcc                                                                                                                                <  2:4201660/150‑1 (MQ=255)
                          cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGtt                                                                                                                             <  2:3124786/150‑1 (MQ=255)
                           aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTc                                                                                                                            <  2:1213319/150‑1 (MQ=255)
                           aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCGAACCCGCTCCCCAAAGGGGGGAGGGGGCCGTTc                                                                                                                            >  1:2962049/1‑150 (MQ=39)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  1:4023618/150‑1 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  1:4023479/150‑1 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  1:3501187/150‑1 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                        <  2:1200841/150‑1 (MQ=255)
                                cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                       <  1:1633994/150‑1 (MQ=255)
                                 ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCAACGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAACCCCCTACCCCAACGGGCGCGGGGACCCTTCAgaaag                                                                                                                      >  1:301841/1‑146 (MQ=37)
                                  tatcTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTCCATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCTGTACgg                                                                                                                     <  2:680706/147‑1 (MQ=255)
                                  tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACgg                                                                                                                     <  1:2627515/150‑1 (MQ=255)
                                    aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCACTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                   >  2:698549/1‑150 (MQ=255)
                                    aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCACCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                   <  1:1035711/150‑1 (MQ=255)
                                      tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCACACCCCCAAAGGGCGAGGGGACCGGTCAGTAAGGaaat                                                                                                                 >  1:2194533/1‑148 (MQ=38)
                                         ggATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGACTTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  2:2064256/149‑1 (MQ=255)
                                          gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                             <  1:120895/150‑1 (MQ=255)
                                           ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACCCCCACACCCACACCCCCACCGGGGGAGGGGCCCGGGCAGTACGGCCTATAGCt                                                                                                            >  2:4131915/1‑150 (MQ=25)
                                           ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCTAAACCTCTCCCCCAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCt                                                                                                            >  1:4210260/1‑150 (MQ=39)
                                             aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACACCCCCAAGGGGCGCGGGGGCCGGTAAGTACGGAATTTAGCTTg                                                                                                          >  1:1802340/1‑150 (MQ=38)
                                              ataaAAGACATTACCCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGa                                                                                                         <  2:740271/150‑1 (MQ=39)
                                              ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGa                                                                                                         <  1:1667575/150‑1 (MQ=255)
                                               taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAACCCAAGCCCCAAAGGGGCGGGGGGGCCCTTCAGTACGGCATTTAGCtttaa                                                                                                        >  1:1739822/1‑147 (MQ=37)
                                               taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                        <  1:2393039/150‑1 (MQ=255)
                                                 aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGACCCTTCAGTACGGAAATTAGCTTGAAgg                                                                                                      >  1:3183852/1‑150 (MQ=39)
                                                 aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCCAAACCCAACCCCAAACGGGGGAGGGGGCCGGTCAGTAAGGGAATTATAGTGACgg                                                                                                      >  1:2004266/1‑150 (MQ=25)
                                                 aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAAATTTAAATCCCGTACGTTCCCCTCCCCCAAACCCACACCCCAAACGGGCGGGGGGCCCCTGCCGTAAGGAAGTTAGCATGaacg                                                                                                      >  1:4272761/1‑147 (MQ=25)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  1:555303/150‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  1:4067771/150‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  1:2071973/150‑1 (MQ=255)
                                                   aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                    <  1:4005539/150‑1 (MQ=255)
                                                       aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAACCTAACCCGCTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGGCGGTGCTAt                                                                                                >  2:732395/1‑150 (MQ=255)
                                                         tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACACCCCAAAGGGGCGCGGGGACCGGGCAGTACCGAAGTTAGCGTGACGGTgggctcc                                                                                              >  2:3356316/1‑144 (MQ=37)
                                                         tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACACTCTACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                              >  2:3963690/1‑150 (MQ=255)
                                                         tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCGAACCCACGCCCCCAAGGGGCGGGGGGCCCCGTCCGTAAGGCATTTAGATTGGCGGGGCTAtac                                                                                              >  2:4094060/1‑148 (MQ=25)
                                                          aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                             <  2:4230893/150‑1 (MQ=255)
                                                            ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                           <  2:2830175/150‑1 (MQ=255)
                                                            ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                           <  1:634106/150‑1 (MQ=255)
                                                            ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTCTCCCCCACGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                           >  2:2743394/1‑150 (MQ=255)
                                                                        tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                               <  1:3221440/150‑1 (MQ=255)
                                                                          aGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACCGAATTTAGCTTGACGGTGCTATCCCGCATTAAATACTGGCa                                                                             >  2:3583007/1‑150 (MQ=255)
                                                                                 ccAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGGACGGTCCCCTCCCCCTCCCCCGCGCCCCCCAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGCCGGTGCTATCCACCATTAAATACTGGCAATAcact                                                                      <  2:3954481/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACCTCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  1:2999035/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  1:3628028/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  2:2553505/150‑1 (MQ=255)
                                                                                        cgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTCACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTTCGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaa                                                               <  1:1493490/150‑1 (MQ=255)
                                                                                        cgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCACTCCCCACAGGGGCGCGGGGCCCCGTCGGGACGGCCTTTGGCGTGGCGGGGGTCTCCACCCGTAAATACGGGCCAGACGCGCCCGCaa                                                               >  2:4078765/1‑150 (MQ=255)
                                                                                         gcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGATGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaa                                                              <  2:108469/150‑1 (MQ=255)
                                                                                         gcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaa                                                              <  2:2836464/150‑1 (MQ=255)
                                                                                          catcCGGCGATGCTGGCGATTTGAATTTTACATCCAGTACGTTCCCCTCACCATAACCCTCTCCCCACAGGGGCGAGTGGACAGTTCATTACGGAATTTAGCTTGACTGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCaaaa                                                             <  1:1120804/150‑1 (MQ=255)
                                                                                           atcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAAt                                                            <  2:1863095/150‑1 (MQ=255)
                                                                                            tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAAACCACTCCCCAAATGGGCGCGGGGACCCGTCAGTACGGACTTTAGCATGACGGGGCTATCCCGCATTAAATAAAGGAAAAACACTCACGCAAAATc                                                           >  2:811860/1‑150 (MQ=255)
                                                                                            tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCTCTACCCAAAGGGGCGAGGGGACCCTTCAGTACGGAAATTAGCTTGACGGTGCTATCCAGCATTAAAAACTGGCAATACTCTCAAGCaaaaaa                                                           >  2:464707/1‑148 (MQ=255)
                                                                                             ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                          <  2:609609/150‑1 (MQ=255)
                                                                                               ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                        <  1:3856530/150‑1 (MQ=255)
                                                                                                gCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgca                                                       >  2:2921773/1‑150 (MQ=255)
                                                                                                    tGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCTACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAAGGTGCTATCCAGCATTAAATACTGGCAATACTATCACGCAAAATCCGCAGCAc                                                   >  2:1818825/1‑150 (MQ=255)
                                                                                                                  aTTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc                                     <  1:714958/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  1:1447008/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:919315/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  1:466716/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  1:285753/150‑1 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCACTTCACCGCTGCGATTCg                            >  2:4315360/1‑150 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGATGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCg                            <  1:2180109/150‑1 (MQ=255)
                                                                                                                              gTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                         <  1:737673/150‑1 (MQ=255)
                                                                                                                                   gtCCCCTTACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCTGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                    <  1:3531916/149‑1 (MQ=255)
                                                                                                                                    tCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                   <  1:898890/150‑1 (MQ=255)
                                                                                                                                      cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                 >  2:2607164/1‑150 (MQ=255)
                                                                                                                                       cctcaccctcacCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAACATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCgta                <  1:539756/150‑3 (MQ=255)
                                                                                                                                       ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                >  2:3867705/1‑150 (MQ=255)
                                                                                                                                       ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAGCCAGCATTAAATCCTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGac                >  2:239695/1‑149 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt               >  1:1409792/1‑150 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt               <  2:4082854/150‑1 (MQ=255)
                                                                                                                                         tCACCCTACCCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc              <  2:3662196/150‑1 (MQ=255)
                                                                                                                                          cACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCg             >  1:1007184/1‑150 (MQ=255)
                                                                                                                                            cccgaccccTCTCCCCACAGGGGCGAGGGGACCGTTCAGTAAGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAGCGTCTCTTCACCTCTGCGTTTCGGGTTAATGCGAATCGca           <  1:4176019/144‑1 (MQ=255)
                                                                                                                                            cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca           <  2:3054165/150‑1 (MQ=255)
                                                                                                                                             ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac          >  1:3853436/1‑150 (MQ=255)
                                                                                                                                              cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca         >  1:3832028/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGACTCGCACACTGt     >  2:402896/1‑150 (MQ=255)
                                                                                                                                                   ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt    >  1:1795504/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc   >  1:1590521/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022

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