Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,973 Δ1 bp 6.5% intergenic (+44/‑16) 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,9730T. 6.5% 316.3 / 14.7 92intergenic (+44/‑16)potI/ybjOputrescine ABC transporter permease/DUF2593 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base . (3/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-02

TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTT  >  NC_000913/897829‑898119
                                                                                                                                                 |                                                                                                                                                  
tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTAGCTACATAATGCCGCGCATGTcgcg                                                                                                                                                >  3:868075/1‑149 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGc                                                                                                                                              <  1:807099/149‑1 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGc                                                                                                                                              >  2:225974/1‑149 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGc                                                                                                                                              <  4:357562/149‑1 (MQ=255)
     ggCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAtt                                                                                                                                           <  2:421858/149‑1 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAttatgta                                                                                                                                      <  4:41344/149‑6 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGtttt                                                                                                                                      <  1:850685/149‑1 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGtttt                                                                                                                                      <  3:786150/149‑1 (MQ=255)
            tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                                    <  3:641065/149‑1 (MQ=255)
             cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCAt                                                                                                                                   <  2:382838/149‑1 (MQ=255)
             cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCAt                                                                                                                                   >  2:906046/1‑149 (MQ=255)
              ggAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATg                                                                                                                                  >  4:915975/1‑149 (MQ=255)
               gAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATgg                                                                                                                                 <  3:440027/149‑1 (MQ=255)
                     tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACg                                                                                                                           >  1:404342/1‑149 (MQ=255)
                      cGGATTTATCGCCTGGTATCTGATGGCAAGCGCAGAAAAACAGCGGATTCGCTATGTCCAGCATGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGACGCGGCATTGTTTTCCTGGAAGACGa                                                                                                                          <  3:695113/149‑1 (MQ=255)
                       ggATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTAGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                         <  4:868075/149‑1 (MQ=255)
                            tATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                                    <  1:214659/149‑1 (MQ=255)
                             aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                                   <  2:768626/149‑1 (MQ=255)
                                gCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGa                                                                                                                >  1:511159/1‑149 (MQ=255)
                                  cTGGTATCTGATGGCTCGCGCAGAAAAACATCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                                    >  4:28308/1‑143 (MQ=255)
                                  cTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGCTACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                                    <  3:28308/143‑1 (MQ=255)
                                     gTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGACGACATTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                           <  3:671023/149‑1 (MQ=255)
                                        tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                        <  4:338116/149‑1 (MQ=255)
                                        tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                        <  4:693597/149‑1 (MQ=255)
                                         cTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAg                                                                                                       >  2:589990/1‑149 (MQ=255)
                                            aTGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaa                                                                                                    >  1:357448/1‑149 (MQ=255)
                                              ggCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAAc                                                                                                  >  1:502529/1‑149 (MQ=255)
                                                       aGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCacc                                                                                         <  3:317653/149‑3 (MQ=255)
                                                       aGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                         <  1:712939/149‑1 (MQ=255)
                                                         aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtctc                                                                                       >  1:419367/1‑149 (MQ=255)
                                                                gCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                                >  4:422558/1‑149 (MQ=255)
                                                                gCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                                <  1:452828/149‑1 (MQ=255)
                                                                acGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAATATCTTCATCTCATGCTCg                                                                                >  2:627459/2‑149 (MQ=255)
                                                                   gATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                             <  1:810358/149‑1 (MQ=255)
                                                                        gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCAtt                                                                                                                                           >  2:119183/1‑82 (MQ=255)
                                                                        gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCAtt                                                                                                                                           <  1:119183/82‑1 (MQ=255)
                                                                           aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                     >  3:360020/1‑149 (MQ=255)
                                                                           aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                     <  2:419367/149‑1 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                                    >  2:822782/1‑84 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                                    <  1:822782/84‑1 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCt                                                                   >  1:267678/1‑149 (MQ=255)
                                                                                cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                                >  1:586297/1‑149 (MQ=255)
                                                                                 aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                               <  2:776360/149‑1 (MQ=255)
                                                                                 aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                               <  2:898396/149‑1 (MQ=255)
                                                                                  gCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                              >  1:97883/1‑149 (MQ=255)
                                                                                   cGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTg                                                             <  4:743533/149‑1 (MQ=255)
                                                                                    gTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                            <  3:881047/149‑1 (MQ=255)
                                                                                     tGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGt                                                           <  3:915975/149‑1 (MQ=255)
                                                                                        aaGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                        >  1:797667/1‑149 (MQ=255)
                                                                                        aaGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                        >  1:904736/1‑149 (MQ=255)
                                                                                                   gAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                             <  3:754433/149‑1 (MQ=255)
                                                                                                   gAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGgcgc                                             >  4:656592/1‑149 (MQ=255)
                                                                                                      gACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCg                                          >  2:154812/1‑149 (MQ=255)
                                                                                                         acTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCa                                       >  2:730958/1‑149 (MQ=255)
                                                                                                          cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                              <  3:30342/125‑1 (MQ=255)
                                                                                                          cTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                              >  4:30342/1‑125 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:141920/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:262907/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  2:676221/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  1:703760/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:837445/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  4:366532/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  2:220451/1‑149 (MQ=255)
                                                                                                            aaaaTTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatta                                    <  1:169698/149‑1 (MQ=255)
                                                                                                               aTTTGCCGACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCAAGGGAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTTCAGGTGGCGGCGCTCGCCAttaata                                 <  4:491686/149‑1 (MQ=255)
                                                                                                               aTTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCAttatta                                 <  1:225974/149‑1 (MQ=255)
                                                                                                                   gCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGAt                             <  1:683542/149‑1 (MQ=255)
                                                                                                                    ccAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATc                            <  2:914869/149‑1 (MQ=255)
                                                                                                                     cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                           >  3:913522/1‑149 (MQ=255)
                                                                                                                     cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTGCGTGGAAGACGAAACGATGGGATTTTTTAAGAAAACATCTTCATCGCATGGTCGCCTGAATGTGCCTGCGCTTGTGTAGGTGGCGGCGCTCGCCATTATTATGATcc                           >  4:570201/1‑149 (MQ=255)
                                                                                                                      aaCCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCg                          >  2:269509/1‑149 (MQ=255)
                                                                                                                           ggCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGcc                     >  2:188123/1‑149 (MQ=255)
                                                                                                                               aCATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa                 >  1:215204/1‑149 (MQ=255)
                                                                                                                                aacaaTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTGAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc                >  1:801838/4‑149 (MQ=255)
                                                                                                                                 aTAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg               >  3:106396/1‑149 (MQ=255)
                                                                                                                                 aTAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg               >  3:116713/1‑149 (MQ=255)
                                                                                                                                 aTAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg               >  3:235676/1‑149 (MQ=255)
                                                                                                                                  tAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGt              >  4:357830/1‑149 (MQ=255)
                                                                                                                                   aaTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg             >  3:154197/1‑149 (MQ=255)
                                                                                                                                    aTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc            >  2:581156/1‑149 (MQ=255)
                                                                                                                                    aTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc            >  4:904198/1‑149 (MQ=255)
                                                                                                                                     tGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt           <  1:906046/149‑1 (MQ=255)
                                                                                                                                      gCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtg          >  1:347345/1‑149 (MQ=255)
                                                                                                                                       ccGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtga         <  4:776638/149‑1 (MQ=255)
                                                                                                                                        cgcgCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat        <  2:1052/149‑1 (MQ=255)
                                                                                                                                        cgcgCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat        <  4:360020/149‑1 (MQ=255)
                                                                                                                                         gcgcATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg       >  3:788294/1‑149 (MQ=255)
                                                                                                                                         gcgcATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg       >  3:490145/1‑149 (MQ=255)
                                                                                                                                             cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtt    <  3:352484/149‑1 (MQ=255)
                                                                                                                                             cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtt    <  1:269509/149‑1 (MQ=255)
                                                                                                                                             cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtt    <  2:630630/149‑1 (MQ=255)
                                                                                                                                               tGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAtttt  <  2:730357/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTT  >  NC_000913/897829‑898119

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: