Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,968:1 +A 7.0% 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 7.0% 349.1 / 22.5 86intergenic (+39/‑21)potI/ybjOputrescine ABC transporter permease/DUF2593 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (45/35);  new base A (3/3);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

ATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATG  >  NC_000913/897824‑898114
                                                                                                                                                 |                                                                                                                                                  
aTTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATg                                                                                                                                                >  3:776638/1‑149 (MQ=255)
aTTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCACAAAAACAGCGGATACGCGACATCCAGCGTGCACGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATg                                                                                                                                                >  4:881047/1‑149 (MQ=255)
     tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTAGCTACATAATGCCGCG‑CATGTcgcg                                                                                                                                           >  3:868075/1‑149 (MQ=255)
       cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGc                                                                                                                                         <  1:807099/149‑1 (MQ=255)
       cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGc                                                                                                                                         <  4:357562/149‑1 (MQ=255)
       cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGc                                                                                                                                         >  2:225974/1‑149 (MQ=255)
          ggCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAtt                                                                                                                                      <  2:421858/149‑1 (MQ=255)
               ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAttatgta                                                                                                                                 <  4:41344/149‑6 (MQ=255)
               ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGtttt                                                                                                                                 <  3:786150/149‑1 (MQ=255)
               ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGtttt                                                                                                                                 <  1:850685/149‑1 (MQ=255)
                 tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                               <  3:641065/149‑1 (MQ=255)
                  cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                              >  2:906046/1‑149 (MQ=255)
                  cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                              <  2:382838/149‑1 (MQ=255)
                   ggAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATg                                                                                                                             >  4:915975/1‑149 (MQ=255)
                    gAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATgg                                                                                                                            <  3:440027/149‑1 (MQ=255)
                          tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACg                                                                                                                      >  1:404342/1‑149 (MQ=255)
                           cGGATTTATCGCCTGGTATCTGATGGCAAGCGCAGAAAAACAGCGGATTCGCTATGTCCAGCATGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGACGCGGCATTGTTTTCCTGGAAGACGa                                                                                                                     <  3:695113/149‑1 (MQ=255)
                            ggATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTAGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                    <  4:868075/149‑1 (MQ=255)
                                 tATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                               <  1:214659/149‑1 (MQ=255)
                                  aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                              <  2:768626/149‑1 (MQ=255)
                                     gCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGa                                                                                                           >  1:511159/1‑149 (MQ=255)
                                       cTGGTATCTGATGGCTCGCGCAGAAAAACATCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                               >  4:28308/1‑143 (MQ=255)
                                       cTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGCTACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                               <  3:28308/143‑1 (MQ=255)
                                          gTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGACGACATTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                      <  3:671023/149‑1 (MQ=255)
                                             tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                   <  4:338116/149‑1 (MQ=255)
                                             tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                   <  4:693597/149‑1 (MQ=255)
                                              cTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAg                                                                                                  >  2:589990/1‑149 (MQ=255)
                                                 aTGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                               >  1:357448/1‑149 (MQ=255)
                                                   ggCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAAc                                                                                             >  1:502529/1‑149 (MQ=255)
                                                            aGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCacc                                                                                    <  3:317653/149‑3 (MQ=255)
                                                            aGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                    <  1:712939/149‑1 (MQ=255)
                                                              aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtctc                                                                                  >  1:419367/1‑149 (MQ=255)
                                                                     gCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           >  4:422558/1‑149 (MQ=255)
                                                                     gCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           <  1:452828/149‑1 (MQ=255)
                                                                     acGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAATATCTTCATCTCATGCTCg                                                                           >  2:627459/2‑149 (MQ=255)
                                                                        gATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                        <  1:810358/149‑1 (MQ=255)
                                                                             gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCAtt                                                                                                                                      <  1:119183/82‑1 (MQ=255)
                                                                             gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCAtt                                                                                                                                      >  2:119183/1‑82 (MQ=255)
                                                                                aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                >  3:360020/1‑149 (MQ=255)
                                                                                aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                <  2:419367/149‑1 (MQ=255)
                                                                                  aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                               <  1:822782/84‑1 (MQ=255)
                                                                                  aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                               >  2:822782/1‑84 (MQ=255)
                                                                                  aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCt                                                              >  1:267678/1‑149 (MQ=255)
                                                                                     cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                           >  1:586297/1‑149 (MQ=255)
                                                                                      aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                          <  2:898396/149‑1 (MQ=255)
                                                                                      aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                          <  2:776360/149‑1 (MQ=255)
                                                                                       gCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                         >  1:97883/1‑149 (MQ=255)
                                                                                        cGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTg                                                        <  4:743533/149‑1 (MQ=255)
                                                                                         gTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                       <  3:881047/149‑1 (MQ=255)
                                                                                          tGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGt                                                      <  3:915975/149‑1 (MQ=255)
                                                                                             aaGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                   >  1:904736/1‑149 (MQ=255)
                                                                                             aaGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                   >  1:797667/1‑149 (MQ=255)
                                                                                                        gAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                        <  3:754433/149‑1 (MQ=255)
                                                                                                        gAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                        >  4:656592/1‑149 (MQ=255)
                                                                                                           gACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCg                                     >  2:154812/1‑149 (MQ=255)
                                                                                                              acTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCa                                  >  2:730958/1‑149 (MQ=255)
                                                                                                               cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                         >  4:30342/1‑125 (MQ=255)
                                                                                                               cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                         <  3:30342/125‑1 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  3:837445/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  3:262907/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  4:366532/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  1:703760/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  2:220451/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  2:676221/1‑149 (MQ=255)
                                                                                                                tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  3:141920/1‑149 (MQ=255)
                                                                                                                 aaaaTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatta                               <  1:169698/149‑1 (MQ=255)
                                                                                                                    aTTTGCCGACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAAGGGAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTTCAGGTGGCGGCGCTCGCCAttaata                            <  4:491686/149‑1 (MQ=255)
                                                                                                                    aTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                            <  1:225974/149‑1 (MQ=255)
                                                                                                                        gCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGAt                        <  1:683542/149‑1 (MQ=255)
                                                                                                                         ccAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATc                       <  2:914869/149‑1 (MQ=255)
                                                                                                                          cAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                      >  3:913522/1‑149 (MQ=255)
                                                                                                                          cAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTGCGTGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCGCATGGTCGCCTGAATGTGCCTGCGCTTGTGTAGGTGGCGGCGCTCGCCATTATTATGATcc                      >  4:570201/1‑149 (MQ=255)
                                                                                                                           aaCCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                     >  2:269509/1‑149 (MQ=255)
                                                                                                                                ggCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGcc                >  2:188123/1‑149 (MQ=255)
                                                                                                                                    aCATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa            >  1:215204/1‑149 (MQ=255)
                                                                                                                                     aacaaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTGAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc           >  1:801838/4‑149 (MQ=255)
                                                                                                                                      aTAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg          >  3:116713/1‑149 (MQ=255)
                                                                                                                                      aTAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg          >  3:235676/1‑149 (MQ=255)
                                                                                                                                      aTAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg          >  3:106396/1‑149 (MQ=255)
                                                                                                                                       tAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGt         >  4:357830/1‑149 (MQ=255)
                                                                                                                                        aaTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg        >  3:154197/1‑149 (MQ=255)
                                                                                                                                         aTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc       >  2:581156/1‑149 (MQ=255)
                                                                                                                                         aTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc       >  4:904198/1‑149 (MQ=255)
                                                                                                                                          tGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt      <  1:906046/149‑1 (MQ=255)
                                                                                                                                           gCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtg     >  1:347345/1‑149 (MQ=255)
                                                                                                                                            ccGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtga    <  4:776638/149‑1 (MQ=255)
                                                                                                                                             cgcg‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat   <  4:360020/149‑1 (MQ=255)
                                                                                                                                             cgcg‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat   <  2:1052/149‑1 (MQ=255)
                                                                                                                                              gcg‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg  >  3:490145/1‑149 (MQ=255)
                                                                                                                                              gcg‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg  >  3:788294/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATG  >  NC_000913/897824‑898114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑