Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,510,336 A→C 6.7% intergenic (+203/+354) 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,3360AC 6.7% 457.8 / 16.1 166intergenic (+203/+354)yjhV/fecEpseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (90/63);  new base C (3/9);  total (94/72)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 170 total aligned reads displayed.

ACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510188‑4510484
                                                                                                                                                    |                                                                                                                                                    
  tttttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACag                                                                                                                                                   <  1:1780267/150‑1 (MQ=255)
    tttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACagag                                                                                                                                                 >  1:1385833/1‑150 (MQ=255)
       aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGa                                                                                                                                              <  1:1688390/150‑1 (MQ=255)
       aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGa                                                                                                                                              >  1:1511849/1‑150 (MQ=255)
       aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGa                                                                                                                                              >  1:913278/1‑150 (MQ=255)
       aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGAGCCAGAGCgg                                                                                                                                              <  2:2099217/150‑2 (MQ=255)
       aGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTACGTGAATTTCCGCTCTGGCACAGAGCGa                                                                                                                                              >  1:1101461/1‑150 (MQ=255)
        gAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCgga                                                                                                                                             <  1:401177/150‑3 (MQ=255)
         aaaaaTCGATGGTCAGATACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTTCGATGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCTTTGGTAATGAACCAGCTCCGTGACTGTCCGCTCTGTGCCAGAGCggac                                                                                                                                            <  2:51941/150‑4 (MQ=255)
         aaaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGTGCCAGAGCggac                                                                                                                                            <  1:1998255/150‑4 (MQ=255)
         aaaaaTCGATGGTCAGAAACTGGAGTCGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGTGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGAAATGAACCAGCTCCGTGAATGTCCACTCTGTGCCAGAGCggac                                                                                                                                            <  1:1731910/150‑4 (MQ=255)
             aTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATATATACTCAAAAGCACAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACACAGCGAAAttat                                                                                                                                        >  1:280120/1‑148 (MQ=255)
              tCGATTGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAttttt                                                                                                                                       <  1:865290/150‑1 (MQ=255)
               cGATGGGCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAtttttt                                                                                                                                      <  2:922284/150‑1 (MQ=255)
                gATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTg                                                                                                                                     <  2:765385/150‑1 (MQ=255)
                gATGGTCAGAAAATGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTg                                                                                                                                     <  2:1031132/150‑1 (MQ=255)
                   ggTCAGAAACTGGATTAGCACTTCCGAGCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAATTATCCTCAAAATCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtc                                                                                                                                  <  2:1435587/150‑1 (MQ=255)
                    gTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtct                                                                                                                                 >  1:1660858/1‑150 (MQ=255)
                     tCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCGCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGAtctc                                                                                                                                >  1:249140/1‑150 (MQ=255)
                      cAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTcc                                                                                                                               >  1:670688/1‑150 (MQ=255)
                       aGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAATTTTTTTGATCTccc                                                                                                                              >  2:425580/1‑150 (MQ=255)
                           aCTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCt                                                                                                                          >  2:1976908/1‑150 (MQ=255)
                                aTTAGCATTTCCGTTCATTTGTTGCTTTTGATTTACGATGGCGTCTGATCATTGATATATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATTTCCGCTCTGTCACTGAGCGAAATTTTTTGATCTCCCACCTGAAAt                                                                                                                     <  1:1840012/150‑1 (MQ=255)
                                aTTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                                     >  2:1346854/1‑150 (MQ=255)
                                aTTAGCAATTCCGTTCCATGGGTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCGCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAt                                                                                                                     <  1:582984/150‑1 (MQ=255)
                                   aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTa                                                                                                                  >  2:990234/1‑150 (MQ=255)
                                   aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTa                                                                                                                  >  2:1360713/1‑150 (MQ=255)
                                    gCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTaa                                                                                                                 <  1:1360794/150‑1 (MQ=255)
                                       aTTCCGTTCCATGGTTGCTTTTGATTTACGTTGGAGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCCGGGCCGTGAATGTCCGCTCTGGCGCAGAGCGATATTTTTTGATGTCCCCCCTGAAAtggaaagt                                                                                                              >  2:1346287/1‑143 (MQ=255)
                                        ttCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACtt                                                                                                             <  2:2068913/150‑1 (MQ=255)
                                            tttCCATTGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTc                                                                                                         <  2:1812910/149‑1 (MQ=255)
                                            gTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTc                                                                                                         <  2:575470/150‑1 (MQ=255)
                                             ttCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCa                                                                                                        <  1:1930227/150‑1 (MQ=255)
                                             gtCCATGGTGGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCCACCTCATTGGTACTGAACCAGCTACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGCCa                                                                                                        <  2:340444/149‑1 (MQ=255)
                                               ccATGGTTGCTTTTGATTTACGATGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGAGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCGAAACTTAggcctg                                                                                                      >  2:362058/1‑145 (MQ=255)
                                                  tGGGTGCGTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCa                                                                                                   <  2:874918/150‑1 (MQ=255)
                                                   ggTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCACCCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACGTAgtcacgtcac                                                                                                  >  2:1141104/1‑150 (MQ=255)
                                                            tGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCTCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTggg                                                                                         >  1:101182/1‑150 (MQ=255)
                                                              aTTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                       >  1:561640/1‑150 (MQ=255)
                                                              aTTTACGTTGGCGTCTGAACATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCACCGTGAATGTCCGCTCTGCCACAGAGCGAACATTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGtt                                                                                       >  2:2078227/1‑150 (MQ=255)
                                                                 tACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCt                                                                                    >  1:373779/1‑150 (MQ=255)
                                                                   cGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGCACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaa                                                                                  >  2:718454/1‑150 (MQ=255)
                                                                   cGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaa                                                                                  >  2:1638644/1‑150 (MQ=255)
                                                                      tGGCTTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCGCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAt                                                                               <  1:1332820/150‑1 (MQ=255)
                                                                       ggCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATAGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTGAAAtt                                                                              >  1:1472066/1‑150 (MQ=255)
                                                                          gTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTa                                                                           <  1:1219237/150‑1 (MQ=255)
                                                                          gTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTa                                                                           <  1:934862/150‑1 (MQ=255)
                                                                            cTGATCATTGATTTATCCTCAAAATCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAAc                                                                         <  2:1762942/150‑1 (MQ=255)
                                                                              gATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACtt                                                                       >  2:1438891/1‑150 (MQ=255)
                                                                              gATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACtt                                                                       >  2:793621/1‑150 (MQ=255)
                                                                                 cATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgc                                                                    >  1:1486245/1‑150 (MQ=255)
                                                                                  aTTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                                   <  2:695955/150‑1 (MQ=255)
                                                                                   ttGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGt                                                                  <  2:296694/150‑1 (MQ=255)
                                                                                   ttGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGt                                                                  >  2:175311/1‑150 (MQ=255)
                                                                                        ttATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTc                                                             <  1:2078887/150‑1 (MQ=255)
                                                                                         tATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCCTGTCACGGTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCg                                                            >  2:1289385/1‑150 (MQ=255)
                                                                                         actCCTCAAAAGCCCAACCTCATTGGTAATGAACCATCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGACACGTTTTTGTGTTTCTAAAATTTTAACTTCGCGTTTTTCg                                                            <  2:1001791/148‑1 (MQ=255)
                                                                                          aTCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGt                                                           >  1:2067734/1‑150 (MQ=255)
                                                                                            ccTCAAAAGCCCAACCTCATTGGTAATGCACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCGGAAATCTAAACGTAGTCCTGTCACGGTTTTGGGTTTCGAAAATTTTAACTTCGCGTTTTTCgggg                                                         >  2:1984505/1‑147 (MQ=255)
                                                                                            ccTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGCGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTg                                                         >  2:147215/1‑150 (MQ=255)
                                                                                            ccTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTg                                                         >  2:1685062/1‑150 (MQ=255)
                                                                                             cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  2:1634529/1‑150 (MQ=255)
                                                                                             cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        <  1:885361/150‑1 (MQ=255)
                                                                                             cTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCGAAACTTAGTCAGGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGc                                                        >  1:1665897/1‑150 (MQ=255)
                                                                                              tCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGcc                                                       <  2:1410002/150‑1 (MQ=255)
                                                                                                aaaaGCCCAAACTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGt                                                     <  1:181645/150‑1 (MQ=255)
                                                                                                 aaaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATATCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGtt                                                    >  1:1976137/1‑149 (MQ=255)
                                                                                                   aGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAg                                                  <  1:822579/150‑1 (MQ=255)
                                                                                                    gCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAgg                                                 >  2:1216327/1‑150 (MQ=255)
                                                                                                       cAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCTTTGCCGTAAGGGtt                                              >  1:1755273/1‑150 (MQ=255)
                                                                                                       cAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                                              >  2:1521144/1‑150 (MQ=255)
                                                                                                       cAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                                              >  1:1114752/1‑150 (MQ=255)
                                                                                                        aaCCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCTCCCCTGAAATCTAATCTTAGTCATGTCACGTGTTTGGGTTTCTAAAATTTTAAATTCGCGTGTTTCGTTGCCGGAAGGGttt                                             >  1:305104/1‑149 (MQ=255)
                                                                                                         aCCTCATTGGTAATGAACCAGCTACGTGAATGTCCTCTCTGTCAAATTGTTTTTTTTTTTTATGTGACCACTGATATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTtat                                            <  2:1600517/150‑1 (MQ=255)
                                                                                                            tCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTACACTTAGTCATGTCCCGTTTTTGGGTTTCTAAAATTTTATCTTCGCGTTTTTCGTTGCCGTAAGGGTTATACa                                         >  2:1444839/1‑150 (MQ=255)
                                                                                                            tCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACa                                         >  2:504688/1‑150 (MQ=255)
                                                                                                               ttGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGaa                                      >  2:685610/1‑150 (MQ=255)
                                                                                                               ttGGTAATGAAACAGCTCCGTGAATGACCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTAACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTATTTGCCGTAAGGGTTATACAGaa                                      >  2:541631/1‑150 (MQ=255)
                                                                                                                tGGTAATGAACCAGCTCCGTGAATGTCGGCTCTGTCACAGAGCGAAATTTTTTGATCTCCCCCCTGATATCTAATCTTTGTCATGTCACGTTTTTGGGTTTATAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGaaa                                     <  1:281229/150‑1 (MQ=255)
                                                                                                                 ggTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAt                                    >  1:1895222/1‑150 (MQ=255)
                                                                                                                    aaTGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTc                                 >  2:1910295/1‑150 (MQ=255)
                                                                                                                        aaCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGtt                             >  1:1271327/1‑150 (MQ=255)
                                                                                                                        aaCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGtt                             <  1:1571228/150‑1 (MQ=255)
                                                                                                                         aCCAGCTCCGTTAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTa                            <  1:2039427/150‑1 (MQ=255)
                                                                                                                         aCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATATTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTATCTAAAATTATAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTa                            >  2:1752163/1‑150 (MQ=255)
                                                                                                                            aGCTCCGTGAATGTTCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCACCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAATTTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                         <  1:333775/150‑1 (MQ=255)
                                                                                                                            aGCTCCGTGAATGTCCGCTCTGTCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGc                         <  1:713231/150‑1 (MQ=255)
                                                                                                                              cTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                       >  2:22341/1‑150 (MQ=255)
                                                                                                                                  gTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGtt                   >  2:1827550/1‑150 (MQ=255)
                                                                                                                                   tGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCCTGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTc                  <  2:868722/150‑1 (MQ=255)
                                                                                                                                    gAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCAAGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTTTACAGAAATGTCAGTTAAGCAGATTTCa                 >  1:1258771/1‑150 (MQ=255)
                                                                                                                                        gTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAACCTTAGTCATGTCCCGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAat             >  1:1838485/1‑150 (MQ=255)
                                                                                                                                        gTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAat             <  1:1442903/150‑1 (MQ=255)
                                                                                                                                         tCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTTGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            <  1:145525/150‑1 (MQ=255)
                                                                                                                                         tCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            >  2:1251000/1‑150 (MQ=255)
                                                                                                                                         tCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCACCTGAAATCTAAACTTAGTCATGGCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCCTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAatt            <  2:82817/150‑1 (MQ=255)
                                                                                                                                                tGGCACAGAGCGAAATTTTTGATCTCCACCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGTGATTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCTTTAAGCAGAGTTCAAAATTGATTGCCg    <  1:330609/150‑1 (MQ=255)
                                                                                                                                                 ggCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCg    >  1:938297/1‑150 (MQ=255)
                                                                                                                                                  gCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGt   >  2:1788739/1‑150 (MQ=255)
                                                                                                                                                   cacaGAGCGATATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTg  >  1:1284689/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTG  >  NC_000913/4510188‑4510484

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