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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,9681.A 9.6% 334.3 / 35.5 94intergenic (+39/‑21)potI/ybjOputrescine ABC transporter membrane subunit PotI/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (52/33);  new base A (7/2);  total (59/35)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

CAACATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGAT  >  NC_000913/897820‑898113
                                                                                                                                                     |                                                                                                                                                 
cAACATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCgcgc                                                                                                                                                  <  1:851320/150‑1 (MQ=255)
  aCATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCAt                                                                                                                                                >  1:1108912/1‑150 (MQ=255)
       gATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTcgc                                                                                                                                           >  1:695890/1‑150 (MQ=255)
       gATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTcgc                                                                                                                                           >  1:1258069/1‑150 (MQ=255)
        aTCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTcgcg                                                                                                                                          >  1:1437264/1‑150 (MQ=255)
         tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCgg                                                                                                                                         <  1:1059713/150‑1 (MQ=255)
          ccTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGc                                                                                                                                        <  2:234193/150‑1 (MQ=255)
            gcGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCAt                                                                                                                                      <  2:192617/149‑1 (MQ=255)
             cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAtt                                                                                                                                     <  2:352611/150‑1 (MQ=255)
                    gTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGAAATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                              >  2:404782/1‑150 (MQ=255)
                     tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTTGCTTAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCAt                                                                                                                             <  2:1286989/150‑1 (MQ=255)
                      cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATg                                                                                                                            <  1:60177/150‑1 (MQ=255)
                      cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATg                                                                                                                            >  2:1518280/1‑150 (MQ=255)
                       ggAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTTGGTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCTTTGTTTTGATgg                                                                                                                           <  2:1572516/150‑1 (MQ=255)
                          aTTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAg                                                                                                                        <  1:1513788/150‑1 (MQ=255)
                          aTTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAg                                                                                                                        <  1:1536974/150‑1 (MQ=255)
                             gTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACg                                                                                                                     >  2:783906/1‑150 (MQ=255)
                               cGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                   >  2:1607403/1‑150 (MQ=255)
                                ggATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGaaa                                                                                                                  >  2:287063/1‑150 (MQ=255)
                                    ttATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                              <  2:1077518/150‑1 (MQ=255)
                                    ttATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                              <  2:946786/150‑1 (MQ=255)
                                             tgTATCTGATGGCTCGCGCAGTAAAACAGGGGATACGGGAAATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCTCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                     <  2:1474453/149‑1 (MQ=255)
                                               tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTa                                                                                                   <  1:1518280/150‑1 (MQ=255)
                                                 tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAg                                                                                                 <  2:771919/150‑1 (MQ=255)
                                                   tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaa                                                                                               <  2:1437264/150‑1 (MQ=255)
                                                   tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGCCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                  >  2:1011444/1‑147 (MQ=255)
                                                   tGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGCCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                  <  1:1011444/147‑1 (MQ=255)
                                                    gATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                              >  2:200512/1‑150 (MQ=255)
                                                    gATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                              >  2:90682/1‑150 (MQ=255)
                                                     aTGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaaa                                                                                             >  1:1262818/1‑150 (MQ=255)
                                                       ggCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACa                                                                                           >  1:860985/1‑150 (MQ=255)
                                                               cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                   >  2:1269112/1‑150 (MQ=255)
                                                               cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                   <  2:1108912/150‑1 (MQ=255)
                                                                   aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                               <  1:1254224/150‑1 (MQ=255)
                                                                    aaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATg                                                                              >  1:206740/1‑150 (MQ=255)
                                                                     atcAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGc                                                                             <  1:1164896/148‑1 (MQ=255)
                                                                      aCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCt                                                                            >  2:214967/1‑150 (MQ=255)
                                                                         gCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGc                                                                         >  1:168971/1‑150 (MQ=255)
                                                                           ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGTAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCa                                                                                     <  1:878332/136‑1 (MQ=255)
                                                                           ggATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCa                                                                                     >  2:878332/1‑136 (MQ=255)
                                                                              tACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGaa                                                                    >  1:315823/1‑150 (MQ=255)
                                                                                cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                  >  1:1008185/1‑150 (MQ=255)
                                                                                cgcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtg                                                                  >  2:1273551/1‑150 (MQ=255)
                                                                                 gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtgt                                                                 >  2:479053/1‑150 (MQ=255)
                                                                                  cgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtgtg                                                                >  2:1133882/1‑150 (MQ=255)
                                                                                  cgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGCAAAAATAATCCCGCGCATGTGGCGGATTTGTTTTCATGGAAGCCGAAACGTTGGGATTTTTTAAGAAAACAGCTGCATCTCATGATCGCCTGAAtgtg                                                                >  2:14871/1‑150 (MQ=255)
                                                                                        ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                          >  2:1169704/1‑150 (MQ=255)
                                                                                        ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                          >  2:811074/1‑150 (MQ=255)
                                                                                         cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                         >  1:997864/1‑150 (MQ=255)
                                                                                         cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                         >  2:1439449/1‑150 (MQ=255)
                                                                                         cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGATGGGATTTTTTAAGAACACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                         >  2:41242/1‑150 (MQ=255)
                                                                                           gCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGATTGTGCCTGCGCTg                                                       >  2:1205053/1‑150 (MQ=255)
                                                                                            cGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                      <  2:1008185/150‑1 (MQ=255)
                                                                                            cGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                      <  2:527009/150‑1 (MQ=255)
                                                                                               gCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGc                                                   >  1:346140/1‑150 (MQ=255)
                                                                                                cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                  >  2:1037498/1‑150 (MQ=255)
                                                                                                cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                  >  1:943591/1‑150 (MQ=255)
                                                                                                   gACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGt                                               <  2:695890/150‑1 (MQ=255)
                                                                                                     cGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTgg                                             <  1:1269112/150‑1 (MQ=255)
                                                                                                      gTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggc                                            <  2:453751/150‑1 (MQ=255)
                                                                                                         gCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                         >  2:1291494/1‑150 (MQ=255)
                                                                                                         gCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                         <  1:287063/150‑1 (MQ=255)
                                                                                                          cTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                      <  1:1603040/136‑1 (MQ=255)
                                                                                                          cTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                      >  2:1603040/1‑136 (MQ=255)
                                                                                                           tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGTGCTGGTGCAGGTGGCGgcgc                                       >  1:1591548/1‑150 (MQ=255)
                                                                                                           tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                       <  2:168971/150‑1 (MQ=255)
                                                                                                            gAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGAGGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTTGTGCAGGTGGCGGCGCt                                      <  2:862014/150‑1 (MQ=255)
                                                                                                              aGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCg                                    <  1:783906/150‑1 (MQ=255)
                                                                                                              aGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCg                                    >  1:1340774/1‑150 (MQ=255)
                                                                                                                 cacTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCa                                 >  1:687374/1‑150 (MQ=255)
                                                                                                                 cacTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCa                                 >  2:418546/1‑150 (MQ=255)
                                                                                                                 cacTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCa                                 >  2:1610782/1‑150 (MQ=255)
                                                                                                                   cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                               >  2:54184/1‑150 (MQ=255)
                                                                                                                     aaaaTTTGCCAACATGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                             >  2:1393597/1‑150 (MQ=255)
                                                                                                                              cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                    >  2:235191/1‑150 (MQ=255)
                                                                                                                              cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                    >  2:1162673/1‑150 (MQ=255)
                                                                                                                                aCCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTg                  <  2:1262818/150‑1 (MQ=255)
                                                                                                                                 ccTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTgg                 >  1:472911/1‑150 (MQ=255)
                                                                                                                                  cTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGc                <  2:191879/150‑1 (MQ=255)
                                                                                                                                   tGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGcc               >  2:859571/1‑150 (MQ=255)
                                                                                                                                      cTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCg            <  1:418546/150‑1 (MQ=255)
                                                                                                                                      cTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCg            <  1:200512/150‑1 (MQ=255)
                                                                                                                                       tACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa           >  2:186726/1‑150 (MQ=255)
                                                                                                                                       aaaacaaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa           >  2:33489/6‑150 (MQ=255)
                                                                                                                                          acaaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGt        >  2:406094/3‑150 (MQ=255)
                                                                                                                                          acaaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGt        >  2:1641129/3‑150 (MQ=255)
                                                                                                                                           caaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                             >  1:1651197/2‑128 (MQ=255)
                                                                                                                                           caaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                             <  2:1651197/127‑1 (MQ=255)
                                                                                                                                            aaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc      >  2:222175/1‑150 (MQ=255)
                                                                                                                                            aaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc      >  2:1662102/1‑150 (MQ=255)
                                                                                                                                             aTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCAGTGGCCTCGACGTGCt     >  1:51804/1‑150 (MQ=255)
                                                                                                                                             aTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt     >  2:1598133/1‑150 (MQ=255)
                                                                                                                                               gCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtga   >  1:1626408/1‑150 (MQ=255)
                                                                                                                                                ccGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat  >  2:1315439/1‑150 (MQ=255)
                                                                                                                                                ccGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat  >  2:733323/1‑150 (MQ=255)
                                                                                                                                                   cgcATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg         <  2:1477547/140‑1 (MQ=255)
                                                                                                                                                   cgcATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg         >  1:1477547/1‑140 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAACATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGAT  >  NC_000913/897820‑898113

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