Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,510,335 C→G 7.6% 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 7.6% 383.7 / 18.1 119intergenic (+202/+355)yjhV/fecEpseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (50/59);  new base G (5/4);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

Only 100 of 122 total aligned reads displayed.

TTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510192‑4510484
                                                                                                                                               |                                                                                                                                                     
tttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACagag                                                                                                                                                 <  1:817947/150‑1 (MQ=255)
      aaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggaca                                                                                                                                           <  1:250419/150‑5 (MQ=255)
      aaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggaca                                                                                                                                           <  2:402978/150‑5 (MQ=255)
         aTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAtttt                                                                                                                                        >  1:408312/1‑150 (MQ=255)
           cGATGGTCAGAAACTGGATTAGCAATTCAGATCCATGGTTGATTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATTAACCAGCTAAGTGAATGTCCGCTATGGAACATAGCGAAAtttttt                                                                                                                                      >  1:1310367/1‑150 (MQ=255)
              tgttCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAt                                                                                                                                   <  1:414164/147‑1 (MQ=255)
              tGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGAGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAt                                                                                                                                   <  2:1268657/150‑1 (MQ=255)
               ggTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggac                                                                                                                                            >  2:89846/1‑137 (MQ=255)
               ggTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggac                                                                                                                                            <  1:89846/140‑4 (MQ=255)
                  cAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGCGCGAAATTTTTTGATCTcc                                                                                                                               >  2:489036/1‑150 (MQ=255)
                  cAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTcc                                                                                                                               >  2:176175/1‑150 (MQ=255)
                   aGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTccc                                                                                                                              <  2:1498219/150‑1 (MQ=255)
                    gAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTcccc                                                                                                                             <  2:717627/150‑1 (MQ=255)
                       aCTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCt                                                                                                                          >  2:654908/1‑150 (MQ=255)
                         tGGATTAGAATTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTATGATAATTGATTTATACTCAAAAGCCCAACCTCATTGGTAATTAACCAGCTCCGTGAATGTCCTCTATGTCAAAGAGCGAAATTTTTTGATCTCCACCCTGa                                                                                                                        >  2:1242902/1‑150 (MQ=255)
                          ggATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGaa                                                                                                                       >  2:806286/1‑150 (MQ=255)
                           gATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGaaa                                                                                                                      >  1:533496/1‑150 (MQ=255)
                             ttAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATc                                                                                                                    <  1:697040/150‑1 (MQ=255)
                               aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTa                                                                                                                  >  2:279258/1‑150 (MQ=255)
                                gCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaa                                                                                                                 >  2:1428038/1‑150 (MQ=255)
                                 cAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaaa                                                                                                                <  1:1557170/150‑1 (MQ=255)
                                 cAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaaa                                                                                                                >  1:119068/1‑150 (MQ=255)
                                  aaTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAc                                                                                                               <  1:474964/150‑1 (MQ=255)
                                   aTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACt                                                                                                              <  1:608642/150‑1 (MQ=255)
                                    ttCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACtt                                                                                                             <  1:1655693/150‑1 (MQ=255)
                                      ccGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAACAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCAAAGAGCGAAATTTTTTGATCTCCCCCCTGCAATCTAAACTGAg                                                                                                           >  1:1153919/1‑150 (MQ=255)
                                       cGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGt                                                                                                          >  2:420845/1‑150 (MQ=255)
                                       cGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGt                                                                                                          >  1:1107905/1‑150 (MQ=255)
                                         ttCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCa                                                                                                        <  1:279258/150‑1 (MQ=255)
                                         ttCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCa                                                                                                        <  1:1247324/150‑1 (MQ=255)
                                         ttCCATGGTTGCTTTTGATTTACGTTGGCGTATGATCATTTATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCa                                                                                                        <  2:1433207/150‑1 (MQ=255)
                                              ttgTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCa                                                                                                   <  1:737288/148‑1 (MQ=255)
                                               ggTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAGTTAGTCATGTccc                                                                                                  >  2:1559682/1‑148 (MQ=255)
                                                   gCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGttt                                                                                              <  2:1683653/150‑1 (MQ=255)
                                                      tttGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCACAACCTCATTGGTAATGAACCAGATCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTACCCCCTGAAATCTAAAATTAGTCATGTCACGTTTTTg                                                                                           >  1:1618602/1‑150 (MQ=255)
                                                        tGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTggg                                                                                         >  1:1145987/1‑150 (MQ=255)
                                                         gATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGt                                                                                        <  1:1138270/150‑1 (MQ=255)
                                                          aTTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                       <  1:1140138/150‑1 (MQ=255)
                                                          aTTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                       <  1:1444685/150‑1 (MQ=255)
                                                          aTTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTAACGATTTTTGGtt                                                                                       >  2:289178/1‑150 (MQ=255)
                                                           tttACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGttt                                                                                      >  1:1403672/1‑150 (MQ=255)
                                                            ttACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTc                                                                                     <  2:1638981/150‑1 (MQ=255)
                                                            ttACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTc                                                                                     <  1:631635/150‑1 (MQ=255)
                                                            ttACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTc                                                                                     <  2:1421259/150‑1 (MQ=255)
                                                             tACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCt                                                                                    >  1:642005/1‑150 (MQ=255)
                                                              aCGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTa                                                                                   >  2:1020725/1‑150 (MQ=255)
                                                                gTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaaa                                                                                 <  2:404045/150‑1 (MQ=255)
                                                                     cGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGATTGTACTCTCTGGCACAGAGCGAATTTTTTTGATCTCCCCCCTTAATTCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAtttt                                                                            >  1:333893/1‑150 (MQ=255)
                                                                        cTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAAc                                                                         >  2:1693343/1‑150 (MQ=255)
                                                                         tGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACt                                                                        >  1:873686/1‑150 (MQ=255)
                                                                         tGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAATTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACt                                                                        >  2:853838/1‑150 (MQ=255)
                                                                          gATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACtt                                                                       >  2:1142539/1‑150 (MQ=255)
                                                                            tCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcg                                                                     <  1:243284/150‑1 (MQ=255)
                                                                            tCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcg                                                                     <  2:1595299/150‑1 (MQ=255)
                                                                              aTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                                   <  1:1020725/150‑1 (MQ=255)
                                                                              aTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                                   <  2:1075316/150‑1 (MQ=255)
                                                                              aTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                                   <  1:1442006/150‑1 (MQ=255)
                                                                               ttGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGTTGTCTAAAATTTTAAATTCGCGt                                                                  >  2:1136930/1‑150 (MQ=255)
                                                                                tGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGtt                                                                 <  2:1100507/150‑1 (MQ=255)
                                                                                 gATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGttt                                                                >  2:1360492/1‑150 (MQ=255)
                                                                                    ttATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTc                                                             <  1:1360492/150‑1 (MQ=255)
                                                                                      aTCCTCACAAGCCCAAACTCATTGGTAATGAAACAGCTCCGAGAAACGACGCTCTGTCAAAGATATATATTTTTTTATCTGACCACTGATATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGt                                                           <  2:1310367/150‑1 (MQ=255)
                                                                                      aTCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGt                                                           <  1:544626/150‑1 (MQ=255)
                                                                                      aTCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGt                                                           >  2:523251/1‑150 (MQ=255)
                                                                                        ccTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTg                                                         <  2:788439/150‑1 (MQ=255)
                                                                                         cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  2:808451/1‑150 (MQ=255)
                                                                                         cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  2:1650651/1‑150 (MQ=255)
                                                                                         cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  2:117446/1‑150 (MQ=255)
                                                                                          tCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTCGTTGCCg                                                      >  1:1704996/1‑150 (MQ=255)
                                                                                           cAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCg                                                      >  2:503468/1‑150 (MQ=255)
                                                                                             aaaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTa                                                    >  1:482570/1‑150 (MQ=255)
                                                                                             aaaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTa                                                    >  1:173084/1‑150 (MQ=255)
                                                                                               aGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAg                                                  <  2:119068/150‑1 (MQ=255)
                                                                                                gCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAgg                                                 <  2:1013896/150‑1 (MQ=255)
                                                                                                 cccAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAggg                                                <  1:117446/150‑1 (MQ=255)
                                                                                                  ccAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGt                                               <  1:480196/150‑1 (MQ=255)
                                                                                                       cTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATAc                                          <  1:806286/150‑1 (MQ=255)
                                                                                                         cATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAg                                        >  2:1603266/1‑150 (MQ=255)
                                                                                                         cATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAg                                        >  1:1124590/1‑150 (MQ=255)
                                                                                                         cATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAg                                        >  1:469529/1‑150 (MQ=255)
                                                                                                                aaTGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCACCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTc                                 <  1:368518/150‑1 (MQ=255)
                                                                                                                   gAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCACTGAAATCTAAACTTAGTCATGTCAAGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTAGTTGCATTAAGGGTTATACAGAAATGtaggt                              >  1:138326/1‑146 (MQ=255)
                                                                                                                    aaCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGtt                             <  2:408312/150‑1 (MQ=255)
                                                                                                                        aGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                         <  2:897898/150‑1 (MQ=255)
                                                                                                                          cTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                       <  2:474398/150‑1 (MQ=255)
                                                                                                                            ccGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCagag                     >  2:1699653/1‑150 (MQ=255)
                                                                                                                             cGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGt                    <  1:176175/150‑1 (MQ=255)
                                                                                                                              gTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGtt                   >  1:1189803/1‑150 (MQ=255)
                                                                                                                                 aaTGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAATTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaa                <  1:853838/150‑1 (MQ=255)
                                                                                                                                   atttcGCTCTGTGCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaaaa              >  1:947334/5‑150 (MQ=255)
                                                                                                                                   atttcGCTCTGTGCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCACGTCCCGTTTTTGGGTTTCTAAAAGTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCaaaa              >  2:1251095/5‑150 (MQ=255)
                                                                                                                                     ttcGCTCTGTGCCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            >  1:792005/3‑150 (MQ=255)
                                                                                                                                     tCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            <  2:482570/150‑1 (MQ=255)
                                                                                                                                     tCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            <  2:173084/150‑1 (MQ=255)
                                                                                                                                      ccGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTACGTTAAGCAGAGTTCAAAattg           >  1:287524/1‑150 (MQ=255)
                                                                                                                                        gCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattgat         <  2:873686/150‑1 (MQ=255)
                                                                                                                                         ctctGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAattgatt        >  1:767570/1‑150 (MQ=255)
                                                                                                                                               cacaGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTTCCGTAAGGGTTATAAAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTg  >  2:1506433/1‑150 (MQ=255)
                                                                                                                                               cacaGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTATTGGGTTTCTAAAATTTTAACTTCGCGATTTTCGTTGCCGTAAGCGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGACGTg  >  1:787621/1‑150 (MQ=255)
                                                                                                                                               cacaGAGAGAAATTTTTTTCTCTCCCCCCTGAAATCTATACTTAGTCATTTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAATGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTg  <  2:256949/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510192‑4510484

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