Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,861 T→A 10.8% intergenic (+78/+74) 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,8610TA10.8% 736.4 / 6.7 272intergenic (+78/+74)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (42/175);  new base A (34/2);  total (95/177)
Fisher's exact test for biased strand distribution p-value = 1.12e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

CGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATC  >  NC_000913/3504713‑3505010
                                                                                                                                                    |                                                                                                                                                      
 ggTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGTATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCa                                                                                                                                                      <  1:754513/150‑1 (MQ=255)
  gTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCacac                                                                                                                                                     >  1:1616055/1‑150 (MQ=255)
  gTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAc                                                                                                                                                     <  1:1592462/150‑1 (MQ=255)
  gTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAc                                                                                                                                                     <  1:5025800/150‑1 (MQ=255)
   tGCACGTCTCCGAACTCTTCTGCCGCGCCAATCTCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAcc                                                                                                                                                    <  2:4289839/150‑1 (MQ=255)
   tGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTACCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCAcc                                                                                                                                                    <  1:6266223/150‑1 (MQ=255)
     cACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCt                                                                                                                                                  <  2:2487461/150‑1 (MQ=255)
         tctcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTTATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCTCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAcc                                                                                                                                              <  1:4873010/150‑1 (MQ=255)
             cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCctct                                                                                                                                          <  2:1265646/150‑1 (MQ=255)
             cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGACCCCAAACCcact                                                                                                                                          >  1:4042434/1‑150 (MQ=39)
                aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACCACCTAACCCACAccc                                                                                                                                       >  1:5673074/1‑150 (MQ=39)
                aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTccc                                                                                                                                       <  2:760970/150‑1 (MQ=255)
                aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTccc                                                                                                                                       <  2:6068599/150‑1 (MQ=255)
                  tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCCCCCCAACccccacccca                                                                                                                                     >  2:5553334/1‑150 (MQ=39)
                  tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCTAAACCACACCCCa                                                                                                                                     >  1:5098191/1‑150 (MQ=39)
                   cttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGACCACACCCAAAACCAAACCCCaa                                                                                                                                    >  1:3097307/1‑150 (MQ=38)
                      ctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACACCCCAAAgg                                                                                                                                 >  1:5263458/1‑150 (MQ=255)
                         ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCGCACCCCCAAGGGGc                                                                                                                              >  2:6181676/1‑150 (MQ=39)
                            cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                           <  1:4383702/150‑1 (MQ=255)
                             gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                          <  2:6018759/150‑1 (MQ=39)
                              ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                         <  2:3128001/150‑1 (MQ=255)
                              ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                         <  2:4713722/150‑1 (MQ=255)
                              ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAggg                                                                                                                         <  1:598957/150‑1 (MQ=39)
                                  tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCACCCGAAACCACAACCCAAAAGGGCCAGGGGccc                                                                                                                    >  2:4961078/1‑150 (MQ=38)
                                  tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGAGAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgggttag                                                                                                                     <  2:4045556/150‑5 (MQ=38)
                                      aTCACCTTCACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTc                                                                                                                <  2:6816948/150‑1 (MQ=255)
                                      aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTc                                                                                                                <  1:5439652/150‑1 (MQ=255)
                                       tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                               <  2:5602308/150‑1 (MQ=255)
                                          ccTTAACTATTGATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                            <  2:4628358/150‑1 (MQ=255)
                                          ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                            >  2:5886799/1‑150 (MQ=255)
                                            ttAACTATTGATAATACAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGCCCCGGACCCTTACCCCCTCCCCAAAGGGGCGGGGGGCCCGTTCCGGACg                                                                                                          >  2:5059547/1‑150 (MQ=38)
                                            ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACACCCCAAAGGGGCGAGGGGGCCGGTAAGTACg                                                                                                          >  2:6687863/1‑150 (MQ=38)
                                            ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTACCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGTCCGTGAAGTAcc                                                                                                          >  1:3633726/1‑149 (MQ=37)
                                              aaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACCCCCAAACCCACACCCCAAAACGGGGAGGGGCCCGGTCCGTACGGa                                                                                                        >  1:3414427/1‑150 (MQ=37)
                                                 tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCGAACCCACAACCCAAAGGGGCGAGGGGGCCGGTAAGTAAgggaat                                                                                                     >  1:5770614/1‑146 (MQ=37)
                                                   ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                   <  1:7121431/150‑1 (MQ=255)
                                                   ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                   <  1:2747703/150‑1 (MQ=255)
                                                   ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                   <  2:1221020/150‑1 (MQ=255)
                                                   ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCAAACCCAAACCCACAACCCCAACGGGCGCGGGGACCGGTCAGTAAGGCATTTa                                                                                                   >  1:3483017/1‑150 (MQ=37)
                                                      ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAACCCACTCCCCCAACGGGCGAGGGGGCCCGTCAGTACGGCATTTATCt                                                                                                >  1:3806155/1‑150 (MQ=37)
                                                           aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACg                                                                                           <  1:6750132/150‑1 (MQ=255)
                                                            aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGTGGACCGTTCAGTACGGAATTTAGCTTGACgg                                                                                          >  2:5028695/1‑150 (MQ=255)
                                                             aaGACCTTAACCGTAGGCCGGCTAAGATGCGCCAGCATCGCATCCGGCGATGCTCGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                         <  2:3962773/150‑1 (MQ=39)
                                                              aGACATGAACCGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                        <  2:4546434/150‑1 (MQ=255)
                                                                   ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACACCCCAAAGGGGCGACGGGACCGTTCCGTACGGAATTTAGCATGACGGTGCTATc                                                                                   >  2:6554163/1‑150 (MQ=38)
                                                                   ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCAACCGGCGATGCTGGCGCGTTGAATTTTAAATCCCGTACGTTCCCCACACCCAAAACCACACCCCCAACGGGCGAGGGGGCCGTTCAGTACGGAAATTAGCGTTACGGTGCGAta                                                                                   >  1:7213941/1‑149 (MQ=37)
                                                                    tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                  <  2:2481093/150‑1 (MQ=255)
                                                                    tAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCAAACCCACAACCCCAACGGGCGGGGGGACCCGTACGTAACGGAATTAGCATTGCGGTGCTAtac                                                                                  >  1:4786795/1‑148 (MQ=19)
                                                                     aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                 <  1:4961078/150‑1 (MQ=255)
                                                                      aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                <  2:3436312/150‑1 (MQ=255)
                                                                      aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                <  2:1944697/150‑1 (MQ=255)
                                                                             gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                         <  2:2082650/150‑1 (MQ=255)
                                                                              ccGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGGCCGTTCCTTAAGGCAATTAGCTTGGAGGTGCTATCCAGCATGAAAt                                                                        >  2:5194763/1‑150 (MQ=38)
                                                                                 gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCACTCCCCCAACGGGCGAGGGGCCCGTTCAGTACCGAATTTAGCGTGACGGGGCTATCCCGCAATAAAAAca                                                                     >  2:6440115/1‑149 (MQ=255)
                                                                                  aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                    <  1:6819473/150‑1 (MQ=255)
                                                                                  aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                    <  1:1396524/150‑1 (MQ=255)
                                                                                  aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                    <  1:4917357/150‑1 (MQ=255)
                                                                                   tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                   <  2:1747868/150‑1 (MQ=255)
                                                                                   tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                   <  2:1384864/150‑1 (MQ=255)
                                                                                   tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                   <  1:6707999/150‑1 (MQ=255)
                                                                                   tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCGCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                   <  2:2320487/150‑1 (MQ=255)
                                                                                    gtccatcgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGc                                                                  <  2:7013918/144‑1 (MQ=255)
                                                                                      gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCTGCATTAAATACTGGCaa                                                                <  2:5453625/150‑1 (MQ=255)
                                                                                      gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAAACCTCTACCCCAATGGGCGGGGGGACCGGTCAGTACGGGCTTTATCATTACGGTGCTCTCCAGCATTAAAAACTGGCaa                                                                >  2:5624124/1‑150 (MQ=255)
                                                                                          cgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCGAACCCGCTACCCAAAGGGGCGAGGGGCCCGGTCAGTACGGAAATTAGCTTGTCGGTGCTAACCAGCATTAAATACTGGCAATaaa                                                            >  2:6253448/1‑148 (MQ=255)
                                                                                           gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                           <  2:220494/150‑1 (MQ=255)
                                                                                           gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                           <  2:612489/150‑1 (MQ=255)
                                                                                              aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                        <  1:1093768/150‑1 (MQ=255)
                                                                                                 atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCTCTCCCCACAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCACTACTCTCACGc                                                     >  2:884745/1‑150 (MQ=255)
                                                                                                 atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGc                                                     <  1:6155050/150‑1 (MQ=255)
                                                                                                        ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                              <  1:4670958/150‑1 (MQ=255)
                                                                                                          ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCAACCAAAACCAAAACCCAAATGGGGGGTGGGTCCGTTGATTAAGGAATTTATGTTTTAGGTTGTATACAGCATTTAATAAATGGAATAATATAAAGCAAAATACgc                                            >  2:2747105/1‑150 (MQ=255)
                                                                                                          ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                            >  2:4823404/1‑150 (MQ=255)
                                                                                                               tGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                       <  2:200585/150‑1 (MQ=255)
                                                                                                                    gcgcGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtct                                  <  2:5702520/150‑1 (MQ=255)
                                                                                                                    gcgcGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTCACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtct                                  <  2:1905140/150‑1 (MQ=255)
                                                                                                                         ttGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTACCCAAAGGGGCGCGGGGACCGTTCAGTACCGAAATTAGCTTTACGGTGCTAACCACCAATAAAAACAGGCAATACACTCACGCAAAAACCGCACCACCGTCTCTTCa                             >  2:6222937/1‑150 (MQ=37)
                                                                                                                          tGAATTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                            <  1:5562288/150‑1 (MQ=255)
                                                                                                                              ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                        <  1:2341562/150‑1 (MQ=255)
                                                                                                                               tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                       <  2:3321762/150‑1 (MQ=255)
                                                                                                                               tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                       <  1:1909344/150‑1 (MQ=255)
                                                                                                                                ttACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGc                      <  1:1638368/150‑1 (MQ=255)
                                                                                                                                ttACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGc                      <  1:4606174/150‑1 (MQ=255)
                                                                                                                                     tCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                 <  1:6937739/150‑1 (MQ=255)
                                                                                                                                     tCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                 <  2:5374508/150‑1 (MQ=255)
                                                                                                                                       ccGTACGTTCCCCACACCCGAAACCTCTCCCCAAAGGGGCGCGGGGCCCGTTCAGTACGGAATTTAGCATGACGGTGCTATCCAGCATTAAAAACTGGCAATACTCTCACCCAAAATCCGCAGCACCGTCTCTTCACCGCAGCGATTAgg               >  2:2981532/1‑150 (MQ=255)
                                                                                                                                         gTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt             >  2:4644165/1‑150 (MQ=255)
                                                                                                                                         gTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt             <  1:930421/150‑1 (MQ=255)
                                                                                                                                          tACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGtt            <  1:1208524/150‑1 (MQ=255)
                                                                                                                                            cGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTATCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa          <  1:4795958/150‑1 (MQ=255)
                                                                                                                                             ggtCCCCTCACCCGAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt         <  1:6181676/148‑1 (MQ=255)
                                                                                                                                               tCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc       <  2:1616055/150‑1 (MQ=255)
                                                                                                                                                 cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa     <  1:85459/150‑1 (MQ=255)
                                                                                                                                                  ccTCACCCTAACCCTCTCCCCCACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa    <  1:6020155/150‑1 (MQ=255)
                                                                                                                                                  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa    >  1:4221327/1‑150 (MQ=255)
                                                                                                                                                  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa    >  1:3570648/1‑150 (MQ=255)
                                                                                                                                                  ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa    <  1:6605962/150‑1 (MQ=255)
                                                                                                                                                   cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt   >  1:1707915/1‑150 (MQ=255)
                                                                                                                                                    tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc  <  1:2226737/150‑1 (MQ=255)
                                                                                                                                                    tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc  >  2:3694296/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATC  >  NC_000913/3504713‑3505010

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: