Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,968:1 +A 6.2% intergenic (+39/‑21) potI → / → ybjO putrescine ABC transporter permease/DUF2593 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,9681.A 6.2% 369.1 / 18.6 96intergenic (+39/‑21)potI/ybjOputrescine ABC transporter permease/DUF2593 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (48/42);  new base A (4/2);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATT  >  NC_000913/897823‑898117
                                                                                                                                                  |                                                                                                                                                     
cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATg                                                                                                                                                   <  2:1293265/150‑1 (MQ=255)
cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATg                                                                                                                                                   >  1:1149128/1‑150 (MQ=255)
    gATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTcgc                                                                                                                                               >  2:979881/1‑150 (MQ=255)
        cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCTCATGTCGCGGCa                                                                                                                                           >  2:288828/1‑150 (MQ=255)
          cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAtt                                                                                                                                         <  1:984749/150‑1 (MQ=255)
           ggCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTg                                                                                                                                        <  2:1552442/150‑1 (MQ=255)
            gcgcGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACATCGGCTACGCGACATCCAGCGTGCAAGACGGGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAttat                                                                                                                                       <  1:827676/150‑3 (MQ=255)
                  tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCAt                                                                                                                                 >  1:462045/1‑150 (MQ=255)
                    ggAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATgg                                                                                                                               <  2:1723543/150‑1 (MQ=255)
                         tGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGAc                                                                                                                          >  2:940900/1‑150 (MQ=255)
                            cGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                       >  2:531690/1‑150 (MQ=255)
                              gATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAAc                                                                                                                     <  2:552059/150‑1 (MQ=255)
                                  tATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                                 >  1:1157580/1‑150 (MQ=255)
                                   aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTgg                                                                                                                >  2:994670/1‑150 (MQ=255)
                                   aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTgg                                                                                                                >  1:1149853/1‑150 (MQ=255)
                                    tCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTggg                                                                                                               <  2:769415/150‑1 (MQ=255)
                                      gCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAt                                                                                                             >  1:1199901/1‑150 (MQ=255)
                                         tGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAtttt                                                                                                          <  1:987495/150‑1 (MQ=255)
                                          ggTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                         >  1:1126247/1‑150 (MQ=255)
                                          ggTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                         >  2:351640/1‑150 (MQ=255)
                                            tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTa                                                                                                       <  2:50709/150‑1 (MQ=255)
                                             aTCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                      <  2:671073/150‑1 (MQ=255)
                                             aTCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                      <  2:327783/150‑1 (MQ=255)
                                               cTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGa                                                                                                    >  1:875639/1‑150 (MQ=255)
                                               cTGAGGGCTCGCGCAGAAAACCAGCGGATACGCGACATCCATCGTGCAAGACGTGGATGAAGACACTAAAATTTGCCAACCTGGCTACATTTTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGa                                                                                                    <  1:859870/150‑1 (MQ=255)
                                                tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaa                                                                                                   >  2:21049/1‑150 (MQ=255)
                                                tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCAACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTTAGaa                                                                                                   >  2:303966/1‑150 (MQ=255)
                                                 gATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                                  >  1:1634899/1‑150 (MQ=255)
                                                      cTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATc                                                                                             >  2:1530132/1‑150 (MQ=255)
                                                         gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCtta                                                                                          >  2:1656647/1‑149 (MQ=255)
                                                         gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTc                                                                                          <  2:1149128/150‑1 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       >  2:1557858/1‑150 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       >  1:522549/1‑150 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCAGGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       <  2:95322/150‑1 (MQ=255)
                                                             aGAAAAACACCGGATACGAGACATCCAGCGTTCAAGACGTGGCTGAAGACTCTAAATTTTGCCACCGTGGGTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGGTGGGATTTTTTAAGAAAACATCTTCAtct                                                                                      <  2:1149525/150‑1 (MQ=255)
                                                               aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCa                                                                                    <  2:1110923/150‑1 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  2:1009229/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  2:1002723/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  1:1424031/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCAGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGAATTGTTTTCATGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   <  1:869611/150‑1 (MQ=255)
                                                                 aaacttcGGATACTCGACTTCCATCTTGCTAGTCGTGTCTGAATACACTAAAATTTGCCATCCTGGCTACATATTGCCGCGCATGTCGCGGCTTTGTTTTCATGGAAGACGTAACGTTGGGATTTTTTAAGAAATCATCTTCATCTCATg                                                                                  <  2:734536/144‑1 (MQ=255)
                                                                 aaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACAATAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATg                                                                                  <  1:771735/150‑1 (MQ=255)
                                                                    cAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTAGCGGCATGGTTTTCATGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCTCATGCTc                                                                               >  2:1266619/1‑150 (MQ=255)
                                                                       cgtaTAAGCGACATCCAGCGGGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGcc                                                                            <  1:1431535/147‑1 (MQ=255)
                                                                        ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                           >  1:1072323/1‑150 (MQ=255)
                                                                        ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                           <  1:1649685/150‑1 (MQ=255)
                                                                         gATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTg                                                                          <  2:1586584/150‑1 (MQ=255)
                                                                            aCGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGACAACCTGGCTACATAATGCCGCGAATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAt                                                                       <  1:1102905/150‑1 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACTCTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      >  2:1387053/1‑150 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      >  1:414539/1‑150 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      <  1:1310261/150‑1 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGACGCTGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGATtg                                                                      >  2:124696/1‑150 (MQ=255)
                                                                                 aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGcc                                                                  <  2:462045/150‑1 (MQ=255)
                                                                                   aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGTCATTGTTTTCATGGAAGACGAAACGTTGGGTTTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTg                                                                <  1:338494/150‑1 (MQ=255)
                                                                                   aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTg                                                                <  1:27995/150‑1 (MQ=255)
                                                                                     ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                              >  1:430545/1‑150 (MQ=255)
                                                                                      cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                             <  2:875639/150‑1 (MQ=255)
                                                                                       aGCGTGCAAGACTTGGCTGAAAACACTAAAATATGCCCACGTGGCTACATTATCCCGCGAATTTCTCGGCATTTTTTTCATTGAAGACGAAACGTTGGGATTTTTTAAGTAAACTTCTTCCTCTCCTGCTCGCCTGAATGTGCCTGCGCt                                                            <  2:258465/150‑1 (MQ=255)
                                                                                         cgggCAAGACGGGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                          <  2:887643/147‑1 (MQ=255)
                                                                                            gCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGc                                                       >  2:663508/1‑150 (MQ=255)
                                                                                             cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTAGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                      >  2:299136/1‑150 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    <  2:535918/150‑1 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    <  2:101582/150‑1 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    >  2:400098/1‑150 (MQ=255)
                                                                                                 aCGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTg                                                  <  1:979881/150‑1 (MQ=255)
                                                                                                      gCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                             >  1:493868/1‑150 (MQ=255)
                                                                                                        tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                           >  1:1413903/1‑150 (MQ=255)
                                                                                                        tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                           <  2:1109686/150‑1 (MQ=255)
                                                                                                         gAAGACACTAAAATTTGCCAACCAGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCATCATCTCATGCTCGCCTGAATGTGCCTGCGCTGTTGCAGGTGGCGgcgcc                                          >  2:752236/1‑149 (MQ=255)
                                                                                                          aaGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTc                                         <  2:223914/150‑1 (MQ=255)
                                                                                                               acTAAAATTTGCCACCCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTATCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCat                                    <  1:313640/150‑1 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAATAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  1:508508/1‑150 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  2:519729/1‑150 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  2:1052490/1‑150 (MQ=255)
                                                                                                                  aaaaTTTGCCAACCTGGATAAATAATGCCGCGCATGTCGAGGCATTTTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGACTGCGCTGGTGAAGTTGTCGTCGCTCGCCaatat                                 >  2:1076070/1‑150 (MQ=255)
                                                                                                                    attttGCCAACCTGGCTAGATAATGCCGCGCATGTCGCGGCATTGTTTTCTTGGACGACGATACGTGGGGTTTTTTTAAGAAAAAATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  2:1050278/148‑1 (MQ=255)
                                                                                                                    aaTTTGCCAAGCTGGCTACATAATGCCGCGCATGTCGCGGGATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCCTCTCATGCTCGCCTGAAGGTGCCTGCGCTGGGGCAGGTGGCGGCGCTCGCCAttatta                               >  2:291545/1‑150 (MQ=255)
                                                                                                                    aaTTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAATCATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  1:1125610/150‑1 (MQ=255)
                                                                                                                    aaTTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  1:1664905/150‑1 (MQ=255)
                                                                                                                     aTTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATtattat                              <  1:400098/150‑1 (MQ=255)
                                                                                                                         gCCAACCTGGCTACATAAATCCGCGCATGTCGCGGCATTGTATTCATGCAAGACGAAACCATGGGATTAATTAACAAAACATCAACATCTCATGCTAGCCTGAATGTGCCTGAGCTCTTTAAGGAGGAATCGCTCTCCATTATTATGATc                          >  2:1409852/1‑150 (MQ=255)
                                                                                                                          ccAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                         >  2:1560036/1‑150 (MQ=255)
                                                                                                                          ccAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                         <  1:420904/150‑1 (MQ=255)
                                                                                                                           cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                        >  2:1000562/1‑150 (MQ=255)
                                                                                                                           cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                        >  2:621290/1‑150 (MQ=255)
                                                                                                                              ccTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTgg                     >  2:943293/1‑150 (MQ=255)
                                                                                                                               cTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGc                    <  1:940900/150‑1 (MQ=255)
                                                                                                                                    aaaacaaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa               >  1:423031/6‑150 (MQ=255)
                                                                                                                                     aCATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc              <  1:1000562/150‑1 (MQ=255)
                                                                                                                                     aCATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCATGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc              >  1:1504037/1‑150 (MQ=255)
                                                                                                                                        caaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg           >  1:910229/2‑150 (MQ=255)
                                                                                                                                         aaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  1:1374769/1‑150 (MQ=255)
                                                                                                                                         aaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCCGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  2:937657/1‑150 (MQ=255)
                                                                                                                                         aaTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  1:1162064/1‑150 (MQ=255)
                                                                                                                                          aTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt         <  1:312876/150‑1 (MQ=255)
                                                                                                                                            gCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGGTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtca       >  1:1654548/1‑148 (MQ=255)
                                                                                                                                              cgcgCATGTCGCGGCATAGTTATCATGGAAGCCGAAACGTTGGGATTATTTAAGAAAACATCTTCATATCATGCTCGCCTGAATCTGAAAGCGCTGGTGCACGTGGCGGCGCTCGCCATTATTATGAACAGTGGCCTCGACGTGCtgatg     >  2:538735/1‑150 (MQ=255)
                                                                                                                                                 gcATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtt  >  1:492554/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATT  >  NC_000913/897823‑898117

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: