Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,003,357 Δ1 bp coding (469/543 nt) ycbU → putative fimbrial protein YcbU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,3560T.98.8% 330.5 / ‑0.1 82coding (468/543 nt)ycbUputative fimbrial protein YcbU
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (39/41);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTGGGGCT  >  NC_000913/1003218‑1003490
                                                                                                                                          |                                                                                                                                       
aGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGta                                                                                                                                        >  1:845530/1‑139 (MQ=255)
aGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGta                                                                                                                                        >  1:373433/1‑139 (MQ=255)
aGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGta                                                                                                                                        <  1:1359797/140‑2 (MQ=255)
 ggCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGtat                                                                                                                                       <  2:778841/140‑3 (MQ=255)
   cGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGtataa                                                                                                                                     <  1:1483456/140‑5 (MQ=255)
      aaGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTc                                                                                                                                 >  1:924757/1‑140 (MQ=255)
       aGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGtat                                                                                                                                       <  1:342864/134‑3 (MQ=255)
       aGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGtat                                                                                                                                       >  2:342864/1‑132 (MQ=255)
          ggtggtATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCcaca                                                                                                                             <  2:1325035/140‑1 (MQ=255)
            tggtATGGCAACGGGGATTGCCGTGCAAATTCTTGTTGCGCAAAGTCACCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACaaa                                                                                                                           <  2:654993/140‑1 (MQ=255)
            tggtATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGAGAACACACCCAAATAGCAACTTGGTATAAGTCCACaaa                                                                                                                           >  1:487297/1‑140 (MQ=255)
               tATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGc                                                                                                                        >  1:1507664/1‑140 (MQ=255)
                aTGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCg                                                                                                                       >  2:312875/1‑140 (MQ=255)
                      aCGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGGTAACACACACAAATATAAACTTCGTATAAGTCCACAAAGGCGGGAGCa                                                                                                                 >  2:530986/1‑140 (MQ=255)
                      aCGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCa                                                                                                                 >  2:1390380/1‑140 (MQ=255)
                      aCGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACACAAGTATTAACTTCGTGTAAGACCACAAAGGCGGGAGCa                                                                                                                 >  2:476695/1‑140 (MQ=255)
                        ggggATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAAc                                                                                                               >  1:1162439/1‑140 (MQ=255)
                        ggggATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAAc                                                                                                               >  2:165299/1‑140 (MQ=255)
                         gggATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGTGGGAGCAACg                                                                                                              >  1:935929/1‑140 (MQ=255)
                         gggATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTACGTATAGGTCCACCAAGGCGGGAGCAACg                                                                                                              >  1:701955/1‑140 (MQ=255)
                          ggATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACgg                                                                                                             >  2:1402784/1‑140 (MQ=255)
                           gATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACggg                                                                                                            <  1:1390380/140‑1 (MQ=255)
                                  gTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAAt                                                                                                     >  2:1011299/1‑140 (MQ=255)
                                         ttCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGAcgg                                                                                              <  1:1056775/140‑1 (MQ=255)
                                          tCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGAcggc                                                                                             >  2:804352/1‑140 (MQ=255)
                                             tGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGt                                                                                          >  1:1408761/1‑140 (MQ=255)
                                              gATGCGCAAAGTCAGCAAGAAATCCCGCTCATTGATGTCGAGCCTCTTCCGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGtt                                                                                         <  2:703352/140‑1 (MQ=255)
                                                tGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtct                                                                                       <  1:46356/140‑1 (MQ=255)
                                                tGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtct                                                                                       <  1:164757/140‑1 (MQ=255)
                                                tGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtct                                                                                       <  1:1457612/140‑1 (MQ=255)
                                                tGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTtct                                                                                       <  2:1526555/140‑1 (MQ=255)
                                                    cAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTAt                                                                                   >  2:1493066/1‑140 (MQ=255)
                                                    cAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTAt                                                                                   >  2:1494516/1‑140 (MQ=255)
                                                         tCAGCAAGAAATCCCTCTCAATCAGTTCCCGCCTCTTTCGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  1:476695/140‑1 (MQ=255)
                                                         tCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGTGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  2:935929/140‑1 (MQ=255)
                                                         tCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  2:1240767/140‑1 (MQ=255)
                                                         tCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  2:259650/140‑1 (MQ=255)
                                                         tCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCCTACTCCCTTAAATGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  2:487297/140‑1 (MQ=255)
                                                         tCAGCAAGAAATCCCCCTCAATCAGATCCCGCCTCTTCCGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGa                                                                              <  1:514751/140‑1 (MQ=255)
                                                           aGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGGTc                                                                            <  2:1507664/140‑1 (MQ=255)
                                                           aGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATc                                                                            <  1:804352/140‑1 (MQ=255)
                                                               aGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGt                                                                        <  2:778153/140‑1 (MQ=255)
                                                                 aaaTCCCGCTCAATCAGATCCATCCTCTTACGCACTTAAGAGCCGGGGGTAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGttt                                                                      <  2:701955/140‑1 (MQ=255)
                                                                   aTCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTAc                                                                    <  2:924757/140‑1 (MQ=255)
                                                                   aTCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTAc                                                                    >  1:68875/1‑140 (MQ=255)
                                                                   aTCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTAc                                                                    <  1:845569/140‑1 (MQ=255)
                                                                   aTCCCGCTCAATCAGATCCAGCCTCCTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTAc                                                                    <  1:576203/140‑1 (MQ=255)
                                                                    tCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTAcc                                                                   >  1:948749/1‑140 (MQ=255)
                                                                     cccGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCa                                                                  <  2:1294285/140‑1 (MQ=255)
                                                                      ccGCTCAACCAGGTCCCTCCTCTTACGCCCATAAAAGCCAGGGATAACTCACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAg                                                                 <  2:737917/140‑1 (MQ=255)
                                                                       cGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGt                                                                >  2:1393676/1‑140 (MQ=255)
                                                                        gCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTg                                                               >  2:1524266/1‑140 (MQ=255)
                                                                         cTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACTGGCGGTAACGCGCCGGCGGTTCTCTATTTTGATCTGGGTTACCAGTGa                                                              >  2:576203/1‑140 (MQ=255)
                                                                          tCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGCTTTCTTTTTTGATCTGGTTTACCCGGGaa                                                             >  2:1476580/1‑140 (MQ=255)
                                                                          tCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCCACGGCGGGTCTTTATTTTTATCTGGTTTACCAGTGaa                                                             >  1:737917/1‑140 (MQ=255)
                                                                               cAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGt                                                        <  1:1132399/140‑1 (MQ=255)
                                                                                       gCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGa                                               <  1:74226/140‑1 (MQ=255)
                                                                                       gCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGa                                               <  1:1393676/140‑1 (MQ=255)
                                                                                            ttACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGa                                                          <  1:1062351/125‑1 (MQ=255)
                                                                                            ttACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGa                                                          >  2:1062351/1‑125 (MQ=255)
                                                                                                  ccTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTTGTCTGGTTTACCAGTGAAGGAGGGATGTTGAAACGGATTATc                                    >  2:255540/1‑140 (MQ=255)
                                                                                                   cTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCt                                   <  1:1402784/140‑1 (MQ=255)
                                                                                                   cTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCt                                   <  1:165299/140‑1 (MQ=255)
                                                                                                   cTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCt                                   <  2:373433/140‑1 (MQ=255)
                                                                                                    ttAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTg                                  >  2:259839/1‑140 (MQ=255)
                                                                                                    ttAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGATTACCAGCGAAAGAAGGATGTCGAAACGGATAATCTg                                  >  1:550813/1‑140 (MQ=255)
                                                                                                         aGCCGGGGATAACACACTCAAAAATCAACTTCGTATAAGACCACAAAGACGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTc                             <  2:550813/140‑1 (MQ=255)
                                                                                                            cGGGGATAACACACTTAAATATTAACTTCTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGt                          <  1:530986/140‑1 (MQ=255)
                                                                                                             ggggATAACACACTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGGGATGATGAAACGGAGAATCTGGATTCTGtt                         >  1:577296/1‑140 (MQ=255)
                                                                                                                  tAACACACTCAAATATCACCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTAt                    <  2:845530/140‑1 (MQ=255)
                                                                                                                        acTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              <  1:1476580/140‑1 (MQ=255)
                                                                                                                        acTCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              <  1:259839/140‑1 (MQ=255)
                                                                                                                          tCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTa            >  1:231633/1‑140 (MQ=255)
                                                                                                                          tCAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGGGAAGGAGTGATGTTGAAACGGGTAAACTGGATTCTGTTCTTATTGGGATTa            >  1:476986/1‑140 (MQ=255)
                                                                                                                           cAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGTGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              >  1:1261658/1‑137 (MQ=255)
                                                                                                                           cAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGTGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              <  2:1261658/137‑1 (MQ=255)
                                                                                                                           cAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTaa           >  2:861802/1‑140 (MQ=255)
                                                                                                                           cAAATATCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAAGAGTGATGGTGACACGAATAATCTGGATTCTGTTCTTATTGGGATTaa           >  2:458004/1‑140 (MQ=255)
                                                                                                                              atatCAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGGTCTCCATTTTGATCTGGTTTACCAGTGAAGGAATGATGTTGTAAACGGTACTCTGGATTCTGTGTTTATTTGTGTTAACGt        >  2:739649/1‑140 (MQ=255)
                                                                                                                                  cAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTgggg    <  2:68875/140‑1 (MQ=255)
                                                                                                                                  cAACTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTgggg    <  2:218587/140‑1 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGGGGTGATGCGACGTCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTGGGGCt  <  2:515215/140‑1 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGGTAACCTGGATTCTGTTTTTATTGGGATTAACGTGGGGCt  >  2:509368/1‑140 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              <  1:84101/128‑1 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTGGGGc   <  2:70484/139‑1 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTGGGGc   >  1:70484/1‑139 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGGGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGAt              >  2:84101/1‑128 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGGTTACCAGTGAAGGAGTGATGATGAAACGGATAATCTGGATTCTGTTATTATTGGGATTAACGTGGGGCt  >  2:510538/1‑140 (MQ=255)
                                                                                                                                    aCTTCGTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGCTCCCTATTTTGATCTGGTTTACCAGCGAAGGAGGGATGGTGAAACGGATAAACTGCATTCTGTTGTTATAGGGATTAACACGGGGCt  >  2:469952/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                       
AGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAAGTCCACAAAGGCGGGAGCAACGGGCGGTAATGCGACGGCGGTTCTCTATTTTGATCTGGTTTACCAGTGAAGGAGTGATGTTGAAACGGATAATCTGGATTCTGTTCTTATTGGGATTAACGTGGGGCT  >  NC_000913/1003218‑1003490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: