Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,510,335 C→G 6.2% intergenic (+202/+355) yjhV → / ← fecE pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,510,3350CG 6.2% 422.1 / 15.5 131intergenic (+202/+355)yjhV/fecEpseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (64/58);  new base G (3/5);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 135 total aligned reads displayed.

ACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510186‑4510484
                                                                                                                                                     |                                                                                                                                                     
acacTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGc                                                                                                                                                       <  1:446957/150‑1 (MQ=255)
 cacTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGca                                                                                                                                                      >  1:1522965/1‑150 (MQ=255)
   cTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCATAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGcaca                                                                                                                                                    <  1:1634820/150‑1 (MQ=255)
       ttAGAAAAATCGATTGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGc                                                                                                                                                <  2:628473/150‑1 (MQ=255)
       ttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGc                                                                                                                                                <  1:539499/150‑1 (MQ=255)
        tAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCg                                                                                                                                               >  1:1777647/1‑150 (MQ=255)
         aGAAAAATCGATCGTCAGAAACTGGATTAGCAATTCCGTTCCATGGATGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCgg                                                                                                                                              <  2:1031499/150‑2 (MQ=255)
          taaaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCgga                                                                                                                                             <  1:666322/149‑3 (MQ=255)
          gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGaa                                                                                                                                             >  2:1402615/1‑150 (MQ=255)
          gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGaa                                                                                                                                             >  2:517622/1‑150 (MQ=255)
          gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGaa                                                                                                                                             >  2:836425/1‑150 (MQ=255)
          gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGGCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGaa                                                                                                                                             <  1:105049/150‑1 (MQ=255)
            aaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggaca                                                                                                                                           <  1:1645479/150‑5 (MQ=255)
              aaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggacatt                                                                                                                                         <  2:411250/150‑7 (MQ=255)
                 cGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAtttttt                                                                                                                                      >  1:1814241/1‑150 (MQ=255)
                              cTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTCCGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTg                                                                                                                         <  2:1131861/150‑1 (MQ=255)
                               tGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGa                                                                                                                        <  1:535770/150‑1 (MQ=255)
                                ggATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGaa                                                                                                                       <  1:412972/150‑1 (MQ=255)
                                  aTTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                                     <  2:1181554/150‑1 (MQ=255)
                                  aTTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                                     <  2:421639/150‑1 (MQ=255)
                                     aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTa                                                                                                                  >  2:272929/1‑150 (MQ=255)
                                     aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTa                                                                                                                  >  2:1692662/1‑150 (MQ=255)
                                      gCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaa                                                                                                                 >  2:597327/1‑150 (MQ=255)
                                      gCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaa                                                                                                                 <  1:1272979/150‑1 (MQ=255)
                                        aaTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATTTCCGCACTGGCACAGAGCGAAATTTTTTGATATCACCCCTGAAATCTaaaa                                                                                                               >  2:1518825/1‑149 (MQ=255)
                                        aaTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGACATCTAAAc                                                                                                               >  2:1336802/1‑150 (MQ=255)
                                        aaTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAc                                                                                                               >  2:1793867/1‑150 (MQ=255)
                                           tCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTa                                                                                                            >  1:1127306/1‑150 (MQ=255)
                                            ccGTTCCTTGGTTGCTTTTGCTTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAg                                                                                                           <  2:329212/150‑1 (MQ=255)
                                            ccGTTCCATGGTTTCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAg                                                                                                           <  2:1522965/150‑1 (MQ=255)
                                            ccGTTCCATGGTGGATTTTGATTTACGTTGGCGTCTGATCCTTGATTTTTCCTCAAAATCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAg                                                                                                           <  2:287664/150‑1 (MQ=255)
                                             cGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCAAATTTTGTTATTTTTTTATGTGACCACTGATATCTAAACTTAGt                                                                                                          <  2:1388827/150‑1 (MQ=255)
                                                  cATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGt                                                                                                     >  2:1762766/1‑150 (MQ=255)
                                                   aTGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTc                                                                                                    >  1:191011/1‑150 (MQ=255)
                                                     ggTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCAc                                                                                                  >  2:1089076/1‑150 (MQ=255)
                                                     ggTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAACTTAGTCATGTCACg                                                                                                 >  2:725595/1‑150 (MQ=255)
                                                          cTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGtttt                                                                                             >  2:477144/1‑150 (MQ=255)
                                                          cTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATATTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGtttt                                                                                             >  2:800961/1‑150 (MQ=255)
                                                           ttttGATTTACGTTGGCTTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGttttt                                                                                            <  2:1510692/150‑1 (MQ=255)
                                                            tttGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCAAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACATAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTg                                                                                           <  1:416892/150‑1 (MQ=255)
                                                              tGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTggg                                                                                         <  1:517622/150‑1 (MQ=255)
                                                              tGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTggg                                                                                         <  2:40709/150‑1 (MQ=255)
                                                                ttttACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                       <  2:589253/149‑1 (MQ=255)
                                                                   tACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCt                                                                                    >  1:781767/1‑150 (MQ=255)
                                                                          gCGTCTGATCATTGATTTATCCTCAAAAGCACATCCTCATTGGTAATGAAACAGCTCCGTGAATGTCCGCTCTGGCTCAGAGCGAAATTTTTTGATCTCACCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAttt                                                                             <  1:125125/150‑1 (MQ=255)
                                                                           cGTCTGTTCATTGATTTATCCTCAATAGCCCAACCTCATTTGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAATTCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAtttt                                                                            <  1:302527/150‑1 (MQ=255)
                                                                           cGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAtttt                                                                            >  1:1251433/1‑150 (MQ=255)
                                                                            gTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTa                                                                           <  2:1811288/150‑1 (MQ=255)
                                                                            gTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGGCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTATACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTa                                                                           <  2:937211/150‑1 (MQ=255)
                                                                             tCTTATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTaa                                                                          <  2:1269487/150‑1 (MQ=255)
                                                                              cTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAAc                                                                         <  2:1739844/150‑1 (MQ=255)
                                                                                 aTCATTGATTTATCCTCAAAAGCCCCACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTCACAGAGCGAAATTTTTTGATCGCCACCCTGAAATCTAAACTTAGTCATGTCACGGTTTTGGGTTTCTTAAATTTTAACTTc                                                                      >  1:280515/1‑150 (MQ=255)
                                                                                 aTCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCATTAAATCTAAAATTAGTCATGTCACGTTTTTGGTTTTATTAAATTTTAACTTc                                                                      >  1:704239/1‑150 (MQ=255)
                                                                                  tCATTGATTTATCCTCAAAAGCCCCAACTTATTTGTAATGTACCCGCTCAGTGGATGTCCGCGATGGCCCATTTCGTAATTTTTTGCTCTACCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcg                                                                     <  1:181818/150‑1 (MQ=255)
                                                                                   cATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgc                                                                    <  1:836425/150‑1 (MQ=255)
                                                                                    aTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTATTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                                   >  1:1005578/1‑150 (MQ=255)
                                                                                       gATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGttt                                                                >  2:440896/1‑150 (MQ=255)
                                                                                       gATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGttt                                                                >  2:1158073/1‑150 (MQ=255)
                                                                                       gATTTATCCTCAAAAGCCCAACCTCATAGGTAATGAACCAGCTCCGTGAATGACCGCTCTGGCACAGAGCGAAGTTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCCTGTCCCGTATTTGGGTTTCTAAAATTTTAAATACGCtttt                                                                >  1:582087/1‑150 (MQ=255)
                                                                                           tATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCg                                                            <  2:781767/150‑1 (MQ=255)
                                                                                            aTCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGt                                                           <  1:707953/150‑1 (MQ=255)
                                                                                               cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  2:973574/1‑150 (MQ=255)
                                                                                                 cAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCg                                                      >  1:812039/1‑150 (MQ=255)
                                                                                                  aaaaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGt                                                     <  1:29580/150‑1 (MQ=255)
                                                                                                      gCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAgg                                                 >  2:1816955/1‑150 (MQ=255)
                                                                                                        ccAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGt                                               >  1:1470680/1‑150 (MQ=255)
                                                                                                         cAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                                              >  1:1700024/1‑150 (MQ=255)
                                                                                                         cAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCACCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                                              >  1:336441/1‑150 (MQ=255)
                                                                                                            ccTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTtata                                           >  1:1077217/1‑150 (MQ=255)
                                                                                                            ccTCAATGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAGTTAGTCATGTCACGTTATTGGGTTTCTAAAATTTTAACTTCGCGTTTATCGTTGCCGTAAGGGTtata                                           <  1:273365/150‑1 (MQ=255)
                                                                                                                   ggTAATGAACCAGCTCCGTGAATGTCCGCTCTGGTACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                                    >  2:1696395/1‑150 (MQ=255)
                                                                                                                   ggTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCCTGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                                    >  1:469105/1‑150 (MQ=255)
                                                                                                                   ggTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                                    >  2:661814/1‑150 (MQ=255)
                                                                                                                    gTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATg                                   <  2:871110/150‑1 (MQ=255)
                                                                                                                       aTGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTcc                                <  1:604305/150‑1 (MQ=255)
                                                                                                                         gAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGt                              >  1:506141/1‑150 (MQ=255)
                                                                                                                          aaCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                          >  1:651374/1‑137 (MQ=255)
                                                                                                                          aaCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                          <  2:651374/137‑1 (MQ=255)
                                                                                                                           aCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTa                            <  1:272929/150‑1 (MQ=255)
                                                                                                                           aCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCAGAAATCTAAACTTAGTCATGTCAAGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTa                            >  1:740781/1‑150 (MQ=255)
                                                                                                                              aGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                         <  2:1419853/150‑1 (MQ=255)
                                                                                                                               gCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCa                        >  1:1517314/1‑150 (MQ=255)
                                                                                                                               gCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCa                        <  1:1793867/150‑1 (MQ=255)
                                                                                                                                cTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                       >  2:174291/1‑150 (MQ=255)
                                                                                                                                cTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                       <  2:1701300/150‑1 (MQ=255)
                                                                                                                                  ccGTGAATGTCCTCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCagag                     <  1:1334713/150‑1 (MQ=255)
                                                                                                                                  ccGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCagag                     <  1:1762766/150‑1 (MQ=255)
                                                                                                                                  ccGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCagag                     <  1:1401661/150‑1 (MQ=255)
                                                                                                                                  ccGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGGTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCagag                     >  1:1115444/1‑150 (MQ=255)
                                                                                                                                   cGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGt                    >  2:418159/1‑150 (MQ=255)
                                                                                                                                     tGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGtta                  >  1:735322/1‑149 (MQ=255)
                                                                                                                                        aatttcGCTCTGTGCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaaa               >  1:1396884/6‑150 (MQ=255)
                                                                                                                                          gTCCGCTCTGGCGCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGGTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAat             >  2:620345/1‑150 (MQ=255)
                                                                                                                                          gTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAat             <  2:523922/150‑1 (MQ=255)
                                                                                                                                          gTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAat             >  2:1034700/1‑150 (MQ=255)
                                                                                                                                            tcGCTCTGTGCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCAACTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTAAAAattt           >  2:605047/2‑149 (MQ=255)
                                                                                                                                             cGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAACTCTAAACTTATTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCTCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattga          <  1:1041386/150‑1 (MQ=255)
                                                                                                                                              gCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGTGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattgat         <  1:1820738/150‑1 (MQ=255)
                                                                                                                                                 ctGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGc      <  1:747657/150‑1 (MQ=255)
                                                                                                                                                   ggCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCg    >  1:1650723/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510186‑4510484

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