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

CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATT  >  NC_000913/897823‑898117
                                                                                                                                                  |                                                                                                                                                     
cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATg                                                                                                                                                   <  2:1293265/150‑1 (MQ=255)
cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATg                                                                                                                                                   >  1:1149128/1‑150 (MQ=255)
    gATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTcgc                                                                                                                                               >  2:979881/1‑150 (MQ=255)
        cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCT‑CATGTCGCGGCa                                                                                                                                           >  2:288828/1‑150 (MQ=255)
          cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAtt                                                                                                                                         <  1:984749/150‑1 (MQ=255)
           ggCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTg                                                                                                                                        <  2:1552442/150‑1 (MQ=255)
            gcgcGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACATCGGCTACGCGACATCCAGCGTGCAAGACGGGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAttat                                                                                                                                       <  1:827676/150‑3 (MQ=255)
                  tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                                 >  1:462045/1‑150 (MQ=255)
                    ggAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATgg                                                                                                                               <  2:1723543/150‑1 (MQ=255)
                         tGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGAc                                                                                                                          >  2:940900/1‑150 (MQ=255)
                            cGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                       >  2:531690/1‑150 (MQ=255)
                              gATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAAc                                                                                                                     <  2:552059/150‑1 (MQ=255)
                                  tATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                                 >  1:1157580/1‑150 (MQ=255)
                                   aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTgg                                                                                                                >  2:994670/1‑150 (MQ=255)
                                   aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTgg                                                                                                                >  1:1149853/1‑150 (MQ=255)
                                    tCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTggg                                                                                                               <  2:769415/150‑1 (MQ=255)
                                      gCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAt                                                                                                             >  1:1199901/1‑150 (MQ=255)
                                         tGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAtttt                                                                                                          <  1:987495/150‑1 (MQ=255)
                                          ggTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                         >  1:1126247/1‑150 (MQ=255)
                                          ggTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                         >  2:351640/1‑150 (MQ=255)
                                            tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTa                                                                                                       <  2:50709/150‑1 (MQ=255)
                                             aTCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                      <  2:671073/150‑1 (MQ=255)
                                             aTCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                      <  2:327783/150‑1 (MQ=255)
                                               cTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGa                                                                                                    >  1:875639/1‑150 (MQ=255)
                                               cTGAGGGCTCGCGCAGAAAACCAGCGGATACGCGACATCCATCGTGCAAGACGTGGATGAAGACACTAAAATTTGCCAACCTGGCTACATTTTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGa                                                                                                    <  1:859870/150‑1 (MQ=255)
                                                tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaa                                                                                                   >  2:21049/1‑150 (MQ=255)
                                                tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCAACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTTAGaa                                                                                                   >  2:303966/1‑150 (MQ=255)
                                                 gATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                                  >  1:1634899/1‑150 (MQ=255)
                                                      cTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATc                                                                                             >  2:1530132/1‑150 (MQ=255)
                                                         gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCtta                                                                                          >  2:1656647/1‑149 (MQ=255)
                                                         gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTc                                                                                          <  2:1149128/150‑1 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       >  2:1557858/1‑150 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       >  1:522549/1‑150 (MQ=255)
                                                            cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CAGGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                       <  2:95322/150‑1 (MQ=255)
                                                             aGAAAAACACCGGATACGAGACATCCAGCGTTCAAGACGTGGCTGAAGACTCTAAATTTTGCCACCGTGGGTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGGTGGGATTTTTTAAGAAAACATCTTCAtct                                                                                      <  2:1149525/150‑1 (MQ=255)
                                                               aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCa                                                                                    <  2:1110923/150‑1 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  2:1009229/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  2:1002723/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   >  1:1424031/1‑150 (MQ=255)
                                                                aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCAGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGAATTGTTTTCATGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                   <  1:869611/150‑1 (MQ=255)
                                                                 aaacttcGGATACTCGACTTCCATCTTGCTAGTCGTGTCTGAATACACTAAAATTTGCCATCCTGGCTACATATTGCCGCG‑CATGTCGCGGCTTTGTTTTCATGGAAGACGTAACGTTGGGATTTTTTAAGAAATCATCTTCATCTCATg                                                                                  <  2:734536/144‑1 (MQ=255)
                                                                 aaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACAATAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATg                                                                                  <  1:771735/150‑1 (MQ=255)
                                                                    cAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTAGCGGCATGGTTTTCATGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCTCATGCTc                                                                               >  2:1266619/1‑150 (MQ=255)
                                                                       cgtaTAAGCGACATCCAGCGGGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGcc                                                                            <  1:1431535/147‑1 (MQ=255)
                                                                        ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                           >  1:1072323/1‑150 (MQ=255)
                                                                        ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                           <  1:1649685/150‑1 (MQ=255)
                                                                         gATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTg                                                                          <  2:1586584/150‑1 (MQ=255)
                                                                            aCGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGACAACCTGGCTACATAATGCCGCG‑AATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAt                                                                       <  1:1102905/150‑1 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACTCTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      >  2:1387053/1‑150 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      >  1:414539/1‑150 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                      <  1:1310261/150‑1 (MQ=255)
                                                                             cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGACGCTGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGATtg                                                                      >  2:124696/1‑150 (MQ=255)
                                                                                 aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGcc                                                                  <  2:462045/150‑1 (MQ=255)
                                                                                   aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGTCATTGTTTTCATGGAAGACGAAACGTTGGGTTTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTg                                                                <  1:338494/150‑1 (MQ=255)
                                                                                   aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTg                                                                <  1:27995/150‑1 (MQ=255)
                                                                                     ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                              >  1:430545/1‑150 (MQ=255)
                                                                                      cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                             <  2:875639/150‑1 (MQ=255)
                                                                                       aGCGTGCAAGACTTGGCTGAAAACACTAAAATATGCCCACGTGGCTACATTATCCCGCG‑AATTTCTCGGCATTTTTTTCATTGAAGACGAAACGTTGGGATTTTTTAAGTAAACTTCTTCCTCTCCTGCTCGCCTGAATGTGCCTGCGCt                                                            <  2:258465/150‑1 (MQ=255)
                                                                                         cgggCAAGACGGGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                          <  2:887643/147‑1 (MQ=255)
                                                                                            gCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGc                                                       >  2:663508/1‑150 (MQ=255)
                                                                                             cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTAGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                      >  2:299136/1‑150 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    <  2:535918/150‑1 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    <  2:101582/150‑1 (MQ=255)
                                                                                               aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                    >  2:400098/1‑150 (MQ=255)
                                                                                                 aCGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTg                                                  <  1:979881/150‑1 (MQ=255)
                                                                                                      gCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                             >  1:493868/1‑150 (MQ=255)
                                                                                                        tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                           >  1:1413903/1‑150 (MQ=255)
                                                                                                        tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                           <  2:1109686/150‑1 (MQ=255)
                                                                                                         gAAGACACTAAAATTTGCCAACCAGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCATCATCTCATGCTCGCCTGAATGTGCCTGCGCTGTTGCAGGTGGCGgcgcc                                          >  2:752236/1‑149 (MQ=255)
                                                                                                          aaGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTc                                         <  2:223914/150‑1 (MQ=255)
                                                                                                               acTAAAATTTGCCACCCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTATCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCat                                    <  1:313640/150‑1 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAATAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  1:508508/1‑150 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  2:519729/1‑150 (MQ=255)
                                                                                                                cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                   >  2:1052490/1‑150 (MQ=255)
                                                                                                                  aaaaTTTGCCAACCTGGATAAATAATGCCGCG‑CATGTCGAGGCATTTTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGACTGCGCTGGTGAAGTTGTCGTCGCTCGCCaatat                                 >  2:1076070/1‑150 (MQ=255)
                                                                                                                    attttGCCAACCTGGCTAGATAATGCCGCG‑CATGTCGCGGCATTGTTTTCTTGGACGACGATACGTGGGGTTTTTTTAAGAAAAAATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  2:1050278/148‑1 (MQ=255)
                                                                                                                    aaTTTGCCAAGCTGGCTACATAATGCCGCG‑CATGTCGCGGGATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCCTCTCATGCTCGCCTGAAGGTGCCTGCGCTGGGGCAGGTGGCGGCGCTCGCCAttatta                               >  2:291545/1‑150 (MQ=255)
                                                                                                                    aaTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAATCATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  1:1125610/150‑1 (MQ=255)
                                                                                                                    aaTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                               <  1:1664905/150‑1 (MQ=255)
                                                                                                                     aTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATtattat                              <  1:400098/150‑1 (MQ=255)
                                                                                                                         gCCAACCTGGCTACATAAATCCGCG‑CATGTCGCGGCATTGTATTCATGCAAGACGAAACCATGGGATTAATTAACAAAACATCAACATCTCATGCTAGCCTGAATGTGCCTGAGCTCTTTAAGGAGGAATCGCTCTCCATTATTATGATc                          >  2:1409852/1‑150 (MQ=255)
                                                                                                                          ccAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                         >  2:1560036/1‑150 (MQ=255)
                                                                                                                          ccAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                         <  1:420904/150‑1 (MQ=255)
                                                                                                                           cAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                        >  2:1000562/1‑150 (MQ=255)
                                                                                                                           cAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                        >  2:621290/1‑150 (MQ=255)
                                                                                                                              ccTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTgg                     >  2:943293/1‑150 (MQ=255)
                                                                                                                               cTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGc                    <  1:940900/150‑1 (MQ=255)
                                                                                                                                    aaaacaaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa               >  1:423031/6‑150 (MQ=255)
                                                                                                                                     aCATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc              <  1:1000562/150‑1 (MQ=255)
                                                                                                                                     aCATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCATGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc              >  1:1504037/1‑150 (MQ=255)
                                                                                                                                        caaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg           >  1:910229/2‑150 (MQ=255)
                                                                                                                                         aaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  1:1374769/1‑150 (MQ=255)
                                                                                                                                         aaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCCGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  2:937657/1‑150 (MQ=255)
                                                                                                                                         aaTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc          >  1:1162064/1‑150 (MQ=255)
                                                                                                                                          aTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt         <  1:312876/150‑1 (MQ=255)
                                                                                                                                            gCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGGTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtca       >  1:1654548/1‑148 (MQ=255)
                                                                                                                                              cgcg‑CATGTCGCGGCATAGTTATCATGGAAGCCGAAACGTTGGGATTATTTAAGAAAACATCTTCATATCATGCTCGCCTGAATCTGAAAGCGCTGGTGCACGTGGCGGCGCTCGCCATTATTATGAACAGTGGCCTCGACGTGCtgatg     >  2:538735/1‑150 (MQ=255)
                                                                                                                                                 g‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtt  >  1:492554/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATT  >  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: ‑