Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,086,965 G→T 7.0% intergenic (+116/‑37) fdhD → / → yiiG formate dehydrogenase formation protein/DUF3829 family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,086,9650GT 7.0% 1251.8 / 10.9 388intergenic (+116/‑37)fdhD/yiiGformate dehydrogenase formation protein/DUF3829 family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (142/205);  new base T (38/3);  total (180/208)
Fisher's exact test for biased strand distribution p-value = 9.44e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

GTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCA  >  NC_000913/4086817‑4087114
                                                                                                                                                    |                                                                                                                                                     
    ataCCCATCCGCAGCGATTGAGCAATTAAATAAAATAAGTCCTGCCAAATATAATAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCAc                                                                                                                                                  <  1:602616/150‑1 (MQ=255)
     taCCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCAct                                                                                                                                                 <  1:3347901/150‑1 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTACTTAAAATAAGTCCTGCGAAATATATTACATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCActctc                                                                                                                                              <  1:4801595/150‑1 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTTTTATGTTAActctc                                                                                                                                              >  1:1891690/1‑150 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCActctg                                                                                                                                              >  2:7693488/1‑149 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCActctc                                                                                                                                              <  2:3711881/150‑1 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTAActctc                                                                                                                                              >  1:1811961/1‑150 (MQ=255)
        ccATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTAttttaaatata                                                                                                                                              >  1:1634184/1‑143 (MQ=255)
         cATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATGTTAActctat                                                                                                                                             >  2:1718176/1‑148 (MQ=255)
             cGCAGCGTTTGAGCTATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         <  1:298772/150‑1 (MQ=255)
             cGCAGCGTTTGAGCACTTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         <  1:7578613/150‑1 (MQ=255)
              gCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCt                                                                                                                                        <  2:5262535/150‑1 (MQ=255)
                aGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtat                                                                                                                                      <  2:4378089/150‑1 (MQ=255)
                aGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtat                                                                                                                                      <  1:1235221/150‑1 (MQ=255)
                aGCGTTTGAGCAAATAAATAAAATAAGTCCTCCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtat                                                                                                                                      <  1:1124446/150‑1 (MQ=255)
                 gCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtata                                                                                                                                     <  1:472708/150‑1 (MQ=255)
                  cGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCtatat                                                                                                                                    <  2:8436034/150‑1 (MQ=255)
                         tcAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTtata                                                                                                                             <  1:4580424/149‑1 (MQ=255)
                          cAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATaa                                                                                                                            <  2:4006785/150‑1 (MQ=255)
                          cAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATaa                                                                                                                            <  2:484423/150‑1 (MQ=255)
                           aattaattAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGTGTGTAGTGGCGTTTTCTTTTTTAATTTATTATTTTAAATCTCTTATGTATATATTTATaat                                                                                                                           >  1:2720075/1‑149 (MQ=255)
                           aattaattAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTGTTTTTTGATTTATTATTTTCACTCTCGTATCTCTATATTTATAAg                                                                                                                           >  2:8371184/1‑150 (MQ=255)
                           aattaattAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTGTTTTTTGAGTTATTATTGTCACGCTCTTATCTAAATATTTATAAg                                                                                                                           >  1:1700987/1‑150 (MQ=255)
                           aattaattAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAg                                                                                                                           >  1:4428450/1‑150 (MQ=255)
                            attaattaAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTGTTGTTTTTTAATTTATTATGTTGAATATCTTTTGTATATATTTATAAgg                                                                                                                          >  1:3247910/1‑150 (MQ=255)
                               aattaaAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCaa                                                                                                                       <  1:7301604/150‑1 (MQ=255)
                                   aaaaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTGTTCGTTTTTCATTTATGATGTTCACTCTCATATCTATATAATTATAAGGCCAATaa                                                                                                                   >  2:7470340/1‑150 (MQ=255)
                                    aaaTAAGTCCTGCGAAATATATTAAATTGAAGGACTTATTCATTTAGTTAATTTTATTATTTTATTTATAAAAAAAATGATTGTATTTTATTTTTATTTTTTCATTTATTATTTTAAATATATTATATATATATTTATAAttaaattaaa                                                                                                                  >  2:5623394/1‑150 (MQ=255)
                                     aaTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAAt                                                                                                                 <  1:5354490/150‑1 (MQ=255)
                                      aTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTAATTTATTATTTTAACTCTCTTATCTATATATTTATAAGGCAATTAAATg                                                                                                                >  1:647335/1‑150 (MQ=255)
                                      aTAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATTTTAAATCTATTATATATATATTTATAATGCAATTAAAtt                                                                                                                >  1:1779956/1‑149 (MQ=255)
                                              tGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTATCTTATCTATATATTTATAAGGCAATTAaatgaaaatgaa                                                                                                        >  2:5118028/1‑150 (MQ=255)
                                               gCGAAATATATGAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAAt                                                                                                       <  2:6556344/150‑1 (MQ=255)
                                                cGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAAtt                                                                                                      >  2:2831692/1‑150 (MQ=255)
                                                    atatatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATttatt                                                                                                  <  1:1073628/150‑1 (MQ=255)
                                                    atatatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCAATCTCTTATCTATATATTTATAAGGCAATTAAAATAAAAGGAAattatt                                                                                                  >  2:86801/1‑150 (MQ=255)
                                                     tatatTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTGTTTTTTGATTTATTTTTTTGACTCTCTTATGTTTATATTTATAAGGGAATTAAATGTAAAGGTAttttttt                                                                                                 >  2:3226736/1‑149 (MQ=255)
                                                         ttAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCtt                                                                                             >  2:2177956/1‑150 (MQ=255)
                                                         ttAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTAATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCtt                                                                                             >  2:5532046/1‑150 (MQ=255)
                                                         ttAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTGTTGTTTTTTAATTTATTATTTTGAAAAAATTATAGTTATATTTATAAAGCAATTAAATTAAAAGGAATttattttatt                                                                                             >  1:8149685/1‑145 (MQ=255)
                                                          tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            <  2:5847744/150‑1 (MQ=255)
                                                          tAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            <  1:4206435/150‑1 (MQ=255)
                                                            aaTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTGATTTATTATTTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTTTATTCCg                                                                                          >  2:6314168/1‑150 (MQ=255)
                                                             aTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAATGGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGc                                                                                         <  1:5166158/150‑1 (MQ=255)
                                                             aTTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTTATTTATTTTTTTGACACTCTTATATATATATTTATAAGGCAATTAAATTAAAAGGAATTTATTATAttaaga                                                                                         >  1:6124813/1‑146 (MQ=255)
                                                                 cAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATTTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATTAAAAGGAATTTATTATATTACGCtatt                                                                                     >  2:3361394/1‑150 (MQ=255)
                                                                  aGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTTATTTATTATTTTTACTAAATTTTAAATATAATTATAATGCCAATAAAATAAAAAGAAATTATTATATTAAGCaatta                                                                                    >  1:3886235/1‑150 (MQ=255)
                                                                    gACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTGTTTTTTAATTTATTATTTTGAAAATCTTATATATATAATTATAATGCCAATAAATGTAAAGGTATTTATTATATTACGCTAAtata                                                                                  >  1:6857301/1‑150 (MQ=255)
                                                                    gACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATtata                                                                                  <  1:3411660/150‑1 (MQ=255)
                                                                      cTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGt                                                                                <  2:6445112/150‑1 (MQ=255)
                                                                         aTTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGc                                                                             <  2:238246/150‑1 (MQ=255)
                                                                         aTTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATTTTGACTATATTTTGTATAAATTTATAAGGCAATTAAATTAAAAAGAAATTATTATATTAAGCAATtatattagc                                                                             >  1:6275359/1‑145 (MQ=255)
                                                                          ttCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGcc                                                                            <  1:6572402/150‑1 (MQ=255)
                                                                             aTTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCcatc                                                                         <  1:4926009/150‑1 (MQ=255)
                                                                              tttCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCcatca                                                                        <  2:2351855/150‑1 (MQ=255)
                                                                              tttCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCGATTATAGTCGCcatca                                                                        >  2:574344/1‑150 (MQ=255)
                                                                               ttCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCcatcat                                                                       <  2:4267947/150‑1 (MQ=255)
                                                                                  gTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTGACTCTCTTATCTATATATTTATAAGGCAATTAAATTAAAATGAATTTATTATATTACGCTATTATAGTAGCCATCATGTc                                                                    >  2:94685/1‑150 (MQ=255)
                                                                                  gTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTc                                                                    <  1:7379756/150‑1 (MQ=255)
                                                                                   tGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            >  2:2044213/1‑125 (MQ=255)
                                                                                   tGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTc                                                                                            <  1:2044213/125‑1 (MQ=255)
                                                                                   tGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTcc                                                                   <  1:5558981/150‑1 (MQ=255)
                                                                                   tGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCAGTTATTATGTTAACTCTCTTATATATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCCTCATGTcc                                                                   >  2:3428212/1‑150 (MQ=255)
                                                                                         ttattattttattTATAAACAAAATGCGTGTAGTGGCTTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGt                                                             <  1:8327294/150‑1 (MQ=255)
                                                                                         ttattattttattTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGt                                                             <  1:8703156/150‑1 (MQ=255)
                                                                                              attttatttATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTGATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGAc                                                        <  1:3072284/150‑1 (MQ=255)
                                                                                               ttttatttATAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACg                                                       >  2:8583633/1‑150 (MQ=255)
                                                                                                     ttatAAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATGTTAACTCTCTTATCTATATATTTATAAGGCAATTAAATTAAAAGGAATTTATTATATTACGCTATTATAGTCGCCATCATGTCCCACGGTCTGACGGGTtgt                                                 >  2:114464/1‑150 (MQ=255)
                                                                                                      tataAACAAAATGCGTGTAGTGGCGTTTTATTTTTTAATTTATTATTTTCACTATCTTTTCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATATTGCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTtgtg                                                >  2:7597753/1‑150 (MQ=255)
                                                                                                       ataAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATGTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTtgtgt                                               >  2:7959903/1‑149 (MQ=255)
                                                                                                         aaaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatg                                             >  2:5397787/1‑150 (MQ=255)
                                                                                                         aaaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCGATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATGATAGTCGCCATCAAGTCCCTCGGTCTGACGGGGTGtgatg                                             >  2:1232248/1‑150 (MQ=255)
                                                                                                         aaaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCGTATCTATATATTTATAAGGCACTTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatg                                             >  2:7355732/1‑150 (MQ=255)
                                                                                                         aaaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTATCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTACGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatg                                             >  2:4690777/1‑150 (MQ=255)
                                                                                                          aaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGtgatga                                            <  1:8090402/150‑1 (MQ=255)
                                                                                                          aaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCAACATGTCCCTCGGTCTGACGGGTTGtgatga                                            >  2:8700316/1‑150 (MQ=255)
                                                                                                          aaCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCATATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTAACTTCCGCTCTTATAGTCGCCCTCATGTCCCTCGGTCTGACGGGTTGtgatga                                            >  2:7102439/1‑150 (MQ=255)
                                                                                                           aCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATc                                                                                                                                         <  1:4332557/56‑1 (MQ=255)
                                                                                                           aCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTgatgat                                           >  2:4933890/1‑150 (MQ=255)
                                                                                                           aCAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATGATAGTCGCCATCATGTCCCTCGGTCGGACGGGGTGTgatgat                                           >  2:1891963/1‑150 (MQ=255)
                                                                                                                 tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  1:5074206/150‑1 (MQ=255)
                                                                                                                 tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  2:2512535/150‑1 (MQ=255)
                                                                                                                 tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  1:7156491/150‑1 (MQ=255)
                                                                                                                 tGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATaaaaaa                                     <  2:2705380/150‑1 (MQ=255)
                                                                                                                        aGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGaaa                              <  1:4973016/150‑1 (MQ=255)
                                                                                                                             cGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGaa                         <  2:6103465/150‑1 (MQ=255)
                                                                                                                             cGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGaa                         <  2:1391669/150‑1 (MQ=255)
                                                                                                                                   cTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTc                   <  2:7005291/150‑1 (MQ=255)
                                                                                                                                   attttttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTc                   <  2:3900757/149‑1 (MQ=255)
                                                                                                                                     tttttCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTccc                 >  1:3506609/1‑150 (MQ=255)
                                                                                                                                           aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           >  1:7934813/1‑150 (MQ=255)
                                                                                                                                           aTTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGc           <  2:3305987/150‑1 (MQ=255)
                                                                                                                                            tttATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGcc          >  1:7782764/1‑150 (MQ=255)
                                                                                                                                            tttATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGcc          >  2:1554874/1‑150 (MQ=255)
                                                                                                                                            tttATTATGTTCACTCTCTTATCTATATATTTATAAGGAATTTAAATTAAAAGTAATTTATTATCTACCGATATTATAGTCGCCATCATGTCCCTCGGTCTGAGGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGcc          <  1:970935/150‑1 (MQ=255)
                                                                                                                                              tattatGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCaa        <  2:3506609/150‑1 (MQ=255)
                                                                                                                                              tattatGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCaa        >  1:6985506/1‑150 (MQ=255)
                                                                                                                                                ttatGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATa      >  1:244280/1‑150 (MQ=255)
                                                                                                                                                ttatGTTCACTCTCTTATCTATATATTTATAAGGAAATTAAATGAAAAGTAATTTATTATCTTCCGCTTTTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATa      <  1:7246738/150‑1 (MQ=255)
                                                                                                                                                    gTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCa  >  2:2894362/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGGCGTTTTCTTTTTTCATTTATTATGTTCACTCTCTTATCTATATATTTATAAGGCAATTAAATGAAAAGGAATTTATTATCTTCCGCTATTATAGTCGCCATCATGTCCCTCGGTCTGACGGGTTGTGATGATAAAAAAGCCGAAACAGAAACGCTCCCGCCTGCCAATAGTCA  >  NC_000913/4086817‑4087114

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: