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

TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTT  >  NC_000913/897829‑898119
                                                                                                                                                 |                                                                                                                                                  
tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTAGCTACATAATGCCGCG‑CATGTcgcg                                                                                                                                                >  3:868075/1‑149 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGc                                                                                                                                              <  1:807099/149‑1 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGc                                                                                                                                              >  2:225974/1‑149 (MQ=255)
  cTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGc                                                                                                                                              <  4:357562/149‑1 (MQ=255)
     ggCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAtt                                                                                                                                           <  2:421858/149‑1 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAttatgta                                                                                                                                      <  4:41344/149‑6 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGtttt                                                                                                                                      <  1:850685/149‑1 (MQ=255)
          ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGtttt                                                                                                                                      <  3:786150/149‑1 (MQ=255)
            tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                                    <  3:641065/149‑1 (MQ=255)
             cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                                   <  2:382838/149‑1 (MQ=255)
             cGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                                   >  2:906046/1‑149 (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                                                                                                                                           >  2:119183/1‑82 (MQ=255)
                                                                        gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCAtt                                                                                                                                           <  1:119183/82‑1 (MQ=255)
                                                                           aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                     >  3:360020/1‑149 (MQ=255)
                                                                           aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGc                                                                     <  2:419367/149‑1 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                                    >  2:822782/1‑84 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCa                                                                                                                                    <  1:822782/84‑1 (MQ=255)
                                                                             aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCt                                                                   >  1:267678/1‑149 (MQ=255)
                                                                                cAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                                >  1:586297/1‑149 (MQ=255)
                                                                                 aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                               <  2:776360/149‑1 (MQ=255)
                                                                                 aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcgc                                                               <  2:898396/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:797667/1‑149 (MQ=255)
                                                                                        aaGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                        >  1:904736/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                                                              <  3:30342/125‑1 (MQ=255)
                                                                                                          cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCt                                                              >  4:30342/1‑125 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:141920/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:262907/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  2:676221/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  1:703760/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  3:837445/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  4:366532/1‑149 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                     >  2:220451/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:106396/1‑149 (MQ=255)
                                                                                                                                 aTAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg               >  3:116713/1‑149 (MQ=255)
                                                                                                                                 aTAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg               >  3:235676/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        <  2:1052/149‑1 (MQ=255)
                                                                                                                                        cgcg‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgat        <  4:360020/149‑1 (MQ=255)
                                                                                                                                         gcg‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg       >  3:788294/1‑149 (MQ=255)
                                                                                                                                         gcg‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg       >  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)
                                                                                                                                                 |                                                                                                                                                  
TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTT  >  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: ‑