Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,453 T→A 5.5% intergenic (+9/+149) fhuB → / ← hemL iron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913173,4530TA 5.5% 1205.6 / 9.1 412intergenic (+9/+149)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base T (184/191);  new base A (34/2);  total (219/193)
Fisher's exact test for biased strand distribution p-value = 2.47e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 428 total aligned reads displayed.

CTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  NC_000913/173304‑173602
                                                                                                                                                     |                                                                                                                                                     
     gggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAc                                                                                                                                                  >  2:4813883/1‑150 (MQ=255)
        tggtTTAGTGCGGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGc                                                                                                                                               <  1:6663057/150‑1 (MQ=255)
         ggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGcc                                                                                                                                              >  2:2246763/1‑150 (MQ=255)
           tttAGTGGTGGTGTTGGGTGAGTGGTGTGGGGGGTTGGTGGTGTTTCCTTTCCTGATCCCGGGGGGGCTGCTGTCAACCTTTATCGGCGCGCAATATTTTATCTTTTTGTTGAGAAAGCAGTGCGGTTAATTTTTCCATATAACGGCCCt                                                                                                                                            <  1:4649019/150‑1 (MQ=255)
           tttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCt                                                                                                                                            >  1:1077939/1‑150 (MQ=255)
           tttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCt                                                                                                                                            >  1:2853194/1‑150 (MQ=255)
            ttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTg                                                                                                                                           >  2:1648399/1‑150 (MQ=255)
                  ctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGcc                                                                                                                                     >  1:4903236/1‑150 (MQ=255)
                  ctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAAATAACGGCCCTGTACGcc                                                                                                                                     >  1:2164289/1‑150 (MQ=255)
                  ctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTAAGcc                                                                                                                                     >  2:5133596/1‑150 (MQ=255)
                  ctgGTGTTCACTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGcc                                                                                                                                     >  2:1444824/1‑150 (MQ=255)
                    ggTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTg                                                                                                                                   >  2:471948/1‑150 (MQ=255)
                      tgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGa                                                                                                                                 <  2:2644821/150‑1 (MQ=255)
                      tgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGa                                                                                                                                 <  2:2535600/150‑1 (MQ=255)
                           gCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCGGTACGCCTGGACGGAt                                                                                                                            >  1:898773/1‑150 (MQ=255)
                            cTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATa                                                                                                                           <  1:5540699/150‑1 (MQ=255)
                            cTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCAGTAAGCCTGGACGGATa                                                                                                                           >  1:2986515/1‑150 (MQ=255)
                            cTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTATCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATa                                                                                                                           >  2:3463952/1‑150 (MQ=255)
                             tGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATAAAAAGGCCCAGTACCCCACGACGGaaaa                                                                                                                          >  2:4164373/1‑150 (MQ=255)
                              gACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa                                                                                                                         <  2:1365972/150‑1 (MQ=255)
                              gACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATAAAACGGCCCAGTACGCCAGGACGGATaaa                                                                                                                         >  1:2768636/1‑150 (MQ=255)
                                cTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc                                                                                                                       <  1:6515965/150‑1 (MQ=255)
                                 tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCg                                                                                                                      >  2:5872109/1‑150 (MQ=255)
                                 tGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATAAAAAGGCCCAGTAACCCAGGACGGGTAAATCg                                                                                                                      >  1:6262136/1‑150 (MQ=255)
                                    tgtgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAAAAAAGGCCCAGTAACCCAGGACGGATAAAGAGttt                                                                                                                   >  2:608501/1‑150 (MQ=255)
                                      tgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTGATTTTTAAATATAAAGGCCCAGTAAGCCTGGGCGGATAAATCGGTTAc                                                                                                                 >  1:1777062/1‑150 (MQ=255)
                                      tgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                                                                 <  2:211728/150‑1 (MQ=255)
                                      tgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                                                                 <  2:3793292/150‑1 (MQ=255)
                                        ggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc                                                                                                               <  1:2249362/150‑1 (MQ=255)
                                        ggCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACCGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc                                                                                                               >  2:5459530/1‑150 (MQ=255)
                                          cGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTAcgcgc                                                                                                             >  2:1075910/1‑149 (MQ=255)
                                            gATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCGGTACGCCTGGACGGATAAATCGTTTACGCGACa                                                                                                           >  2:6474570/1‑150 (MQ=255)
                                             aTGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACaa                                                                                                          <  2:2300010/150‑1 (MQ=255)
                                                 tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGa                                                                                                      >  2:1231851/1‑150 (MQ=255)
                                                 tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCTGTAAGCCTGGACGGATAAAGCGTTTACGCGACAACCGa                                                                                                      >  2:3000301/1‑150 (MQ=255)
                                                 tGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGCCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGa                                                                                                      <  1:4071239/150‑1 (MQ=255)
                                                   cTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCATATAAAGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATg                                                                                                    >  2:817407/1‑150 (MQ=255)
                                                      tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTCGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATTccc                                                                                                 >  1:3839018/1‑150 (MQ=255)
                                                        tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGa                                                                                               <  2:4517598/150‑1 (MQ=255)
                                                        tccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAATCGTTTACGCGACAACGGATGCCCGa                                                                                               >  1:4300167/1‑150 (MQ=255)
                                                        tccatgccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTATTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGa                                                                                               <  1:2303373/150‑1 (MQ=255)
                                                         ccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                                                              <  1:2603660/150‑1 (MQ=255)
                                                         ccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                                                              <  1:3188821/150‑1 (MQ=255)
                                                         ccattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCGGTACGCCAGGACGGATAAAGCGGTTACGCGACAACGGATGCCCCAt                                                                                              >  1:4742140/1‑150 (MQ=255)
                                                          cattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGACAAAGGATGCCCGATg                                                                                             >  2:5204310/1‑150 (MQ=255)
                                                          cattccaGATCCCGGCGGGGCAGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTCCCTAAAACGGCCCTGGGCGCCGGGCCGGCAAAACCGGTTACGCGACAACCGGTGCCCGagg                                                                                             >  2:6244425/1‑148 (MQ=255)
                                                             tccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGTATTTTTGCATATAAAGGCCCGGTACGCCGGGGCGGATAAAGCGGTTTAGCGCCAACCGATGCCCCCTGCgc                                                                                          >  2:2873378/1‑149 (MQ=255)
                                                              ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                                                         <  1:4980306/150‑1 (MQ=255)
                                                                aGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGACAACGGATGCCCGATGCGACGc                                                                                       >  2:5601731/1‑150 (MQ=255)
                                                                  aTCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTg                                                                                     <  2:6542422/150‑1 (MQ=255)
                                                                  aTCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTACGCCAGGACGGATAAAGCGTTTAAGCGACAACGGATGCCCGATGCGACGCTg                                                                                     >  2:1939954/1‑150 (MQ=255)
                                                                   tCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTACGCCAGGTCGGATAAAGCGTTTTCGCGACAACGGATGCCCGAGGCGCCGCggg                                                                                    >  1:3467354/1‑150 (MQ=255)
                                                                     ccGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAAGCGACGCTGgcg                                                                                  >  2:4299933/1‑150 (MQ=255)
                                                                     ccGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCATTACGCCATGACGGATAAAGCGTTTACGCGAAAACCGATGCCCGATTCGACTCTGgcg                                                                                  >  2:3714568/1‑150 (MQ=255)
                                                                      cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcgc                                                                                 <  2:5726628/150‑1 (MQ=255)
                                                                      cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcgc                                                                                 <  2:5510874/150‑1 (MQ=255)
                                                                      cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCAGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATGCCCGATGCGACGCTGgcgc                                                                                 >  2:3432218/1‑150 (MQ=255)
                                                                      cggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCTGTACCCCAGGACGGATAAATCGTTTTCGCGACAACGGATGCCCGATGCGACGCAGgcgc                                                                                 >  2:595497/1‑150 (MQ=255)
                                                                        gcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTAGTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGt                                                                               <  1:770368/150‑1 (MQ=255)
                                                                          ggggCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCt                                                                             <  1:447381/150‑1 (MQ=255)
                                                                              ctgctgTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATc                                                                         <  1:1648399/150‑1 (MQ=255)
                                                                                  tgTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCATATAAAGGCCCAGTAAGCCAGGACGGATAAAACGTTTAACACACAACCGATTCCCCAAGCGGCGCTGGCGCGTCTTAAAGGGc                                                                     >  1:674716/1‑150 (MQ=255)
                                                                                    tCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTACATATAAAGGCCCAGTACGCCAGGACGGATAAAGAGTTTACGCGACCAACGATTCCCCATTCGACGCAGGCCCGTCTTATAGGGcca                                                                   >  1:3923959/1‑149 (MQ=255)
                                                                                      aaCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAc                                                                 <  1:1427867/150‑1 (MQ=255)
                                                                                      aaCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAc                                                                 <  1:2047325/150‑1 (MQ=255)
                                                                                      aaCCTTTATCGGCGCGCCATATTTTATCTATGTGTTGAGACAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAc                                                                 <  1:4164373/150‑1 (MQ=255)
                                                                                       aCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACa                                                                >  2:5375723/1‑150 (MQ=255)
                                                                                         cTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCAGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaaa                                                              >  2:3164892/1‑150 (MQ=255)
                                                                                         cTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTAAAAATAAAGGCCCGGTACCCCTGGACGGCTAAAGCGTTTACGCGACAACGGCTGCCCCATGCGACGCAGGCGCGTCTTATCGGGCCGACaaa                                                              >  1:6027739/1‑150 (MQ=255)
                                                                                          tttATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAAc                                                             <  1:115667/150‑1 (MQ=255)
                                                                                           ttATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTAAAAAAAAAGGCCCAGGACGCCAGGGCGGGTAAAACCTTTAAGCGACAACGGATGCCCCCAGCGACCCGGGCGCGTCGTATAGGGCCAACAAAcc                                                            >  1:6116927/1‑149 (MQ=255)
                                                                                            tATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACgg                                                           <  2:2380130/150‑1 (MQ=255)
                                                                                             aTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCAACAAACGGc                                                          <  2:5080574/150‑1 (MQ=255)
                                                                                                    cgcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCGGTACGCCTGGACGGATAAAGCGGTTACGCGACAACCGATGCCCGATGCGACGCTGGCGCGTCTTATAGGGCCTACAAACCGCCCCGAAt                                                   >  2:687381/1‑150 (MQ=255)
                                                                                                     gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATc                                                  <  2:4635277/150‑1 (MQ=255)
                                                                                                     gcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCTTTACGCCTGGAAGGATAAAGCGTTTAAGCGACAACGGATTCCCGATTCGACGCTTGCGCGTATTATAGGGCCTACAAAACGCCCCCaaac                                                  >  2:1517191/1‑148 (MQ=255)
                                                                                                        atatTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTa                                               <  2:1812078/150‑1 (MQ=255)
                                                                                                           ttttATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTAAGCCTGGAAGGATAAAGCCTTTACGCGACAACGGATTACCCATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAACGTAGGt                                            >  2:891636/1‑150 (MQ=255)
                                                                                                               aTCTATTTGTTGAGAAAGCAGAGCCGTTAATTTGTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGa                                        <  1:3718619/150‑1 (MQ=255)
                                                                                                                tCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATAGGGCCTAAAAACGGCCCCGAATCGTAGGTCGGAt                                       >  2:1616975/1‑150 (MQ=255)
                                                                                                                   aTTTGTTGAGAAAGCTTAGCGTTTATTTTTTCAATATACCGGCCCTGTACTCCTGGACGGATAAAGCGTTTACGCGAAAACGGATGCCCGTTGCGACGCTGGCGCGTGTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  1:5683182/150‑1 (MQ=255)
                                                                                                                   aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  1:3511793/150‑1 (MQ=255)
                                                                                                                   aTTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTATTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAg                                    <  2:6163131/150‑1 (MQ=255)
                                                                                                                    tttGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCAGTACGCCAGGACGGATAAAGCGTTTAAGCGACAACGGATTCCCGATGCCACGCTGGCGCGTCTTATAGGGCCAACAAAAGGCCCCCGATCGTAGGTCGGGTAAgg                                   >  2:1172573/1‑150 (MQ=255)
                                                                                                                     ttgtggAGAAAGCAGCGCCGTTAATTTTTCCATATACCGGCCCTGTCCGCCTGGCCGGATAAAGCGTTTACGCGACAACGGAGGCCCGATGCGACGCGGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGc                                  <  1:4813812/150‑1 (MQ=255)
                                                                                                                             aaaGCAGAGCCGTTAATTTTTCCATATACCGGCGCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgc                          <  2:2911066/150‑1 (MQ=255)
                                                                                                                             aaaGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAAAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgc                          >  2:6717282/1‑150 (MQ=255)
                                                                                                                                 cAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgcctta                      <  2:4743264/150‑4 (MQ=39)
                                                                                                                                  agagCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAAcgccttat                     <  1:5914734/150‑5 (MQ=39)
                                                                                                                                  agagCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAACGTAGGTCGGATAAGGCGTTCAAGCCGCAt                     >  2:2773874/1‑150 (MQ=255)
                                                                                                                                   gagCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATc                    <  2:3665649/150‑1 (MQ=255)
                                                                                                                                      ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                 <  2:926004/150‑1 (MQ=255)
                                                                                                                                        gTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAcc               <  2:4045708/150‑2 (MQ=39)
                                                                                                                                        gTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa               >  1:124155/1‑150 (MQ=255)
                                                                                                                                         ttAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa              >  1:2693687/1‑150 (MQ=255)
                                                                                                                                             tttttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCa          <  1:4363630/150‑1 (MQ=255)
                                                                                                                                                ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGc       >  2:3939467/1‑150 (MQ=255)
                                                                                                                                                ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGc       >  1:4591879/1‑150 (MQ=255)
                                                                                                                                                     tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGACTCGTCGGTCGGAGACGGCGTTCACGCCGCCTCCGCCACACCATGCTggcg  >  2:4373176/1‑148 (MQ=38)
                                                                                                                                                     tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt  <  1:5257361/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAT  >  NC_000913/173304‑173602

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