Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 142,754 A→C intergenic (‑84/‑25) can ← / → yadG carbonic anhydrase/putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913142,7540AC99.1% 785.2 / ‑2.2 222intergenic (‑84/‑25)can/yadGcarbonic anhydrase/putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (121/99);  total (123/99)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 222 total aligned reads displayed.

TCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTAAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGTGGGAT  >  NC_000913/142657‑142849
                                                                                                 |                                                                                               
tCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAg                                                                                               >  2:4046091/1‑100 (MQ=255)
    tGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACg                                                                                           <  2:2487834/100‑1 (MQ=255)
    tGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACg                                                                                           >  2:835299/1‑100 (MQ=255)
      tCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgta                                                                                         <  2:1738960/100‑1 (MQ=255)
       cTTTGATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaa                                                                                        >  2:4648993/1‑100 (MQ=255)
       cTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaa                                                                                        <  1:3130603/100‑1 (MQ=255)
       cTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaa                                                                                        >  2:1527312/1‑100 (MQ=255)
        tttCATGGAGGTTAACGAGCCGTAACGAAATAATTACGTTTGGCTCATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaag                                                                                       <  1:5210066/100‑1 (MQ=255)
        tttCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaag                                                                                       <  2:186903/100‑1 (MQ=255)
           cATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggta                                                                                    <  2:3405154/100‑1 (MQ=255)
           cATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggta                                                                                    >  2:1422604/1‑100 (MQ=255)
            aTGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggtaa                                                                                   >  1:4818064/1‑100 (MQ=255)
            aTGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggtaa                                                                                   >  1:4878717/1‑100 (MQ=255)
             tGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggtaag                                                                                  <  2:4872465/100‑1 (MQ=255)
                aGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACgtaaggtagatcg                                                                               >  1:2156071/1‑95 (MQ=255)
                 ggTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTaaa                                                                              >  2:2365710/1‑100 (MQ=255)
                   ttAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTaaaaa                                                                            <  2:3257413/100‑1 (MQ=255)
                    tAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAAt                                                                           >  1:3265757/1‑100 (MQ=255)
                     aaCGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAGCGTAAGGTAAGTAAAAAtt                                                                          >  2:18152/1‑100 (MQ=255)
                     aaCGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAAtt                                                                          >  1:2314815/1‑100 (MQ=255)
                     aaCGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGAGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAAtt                                                                          <  2:1461963/100‑1 (MQ=255)
                      aCGACCTGTAACCAAATAATTACGTGTGGCTAAGAGAGGGCAACTCCGGGACGAGTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAAgtt                                                                         >  2:809078/1‑100 (MQ=255)
                       cGACCTGTAACCAAATAATTACGTGTGGCCAATATAGAGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTa                                                                        <  1:5564206/100‑1 (MQ=255)
                         aCCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTTCGGGACGATTTAAACCACATATAAATTGTCAGAACGTAAGGTAAGTAAAAATTTATg                                                                      >  2:115000/1‑100 (MQ=255)
                         aCCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATg                                                                      <  2:3840855/100‑1 (MQ=255)
                         aCCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATg                                                                      >  1:111541/1‑100 (MQ=255)
                          ccTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGa                                                                     >  1:1803953/1‑100 (MQ=255)
                          ccTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGa                                                                     >  1:3366244/1‑100 (MQ=255)
                           cTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGAc                                                                    <  1:1251335/100‑1 (MQ=255)
                            tGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGAcc                                                                   <  1:373515/100‑1 (MQ=255)
                            tGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGAcc                                                                   >  2:5100363/1‑100 (MQ=255)
                             gTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCa                                                                  >  2:3740514/1‑100 (MQ=255)
                             gTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCa                                                                  <  1:2064742/100‑1 (MQ=255)
                              tAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCAt                                                                 <  2:2046979/100‑1 (MQ=255)
                              tAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCAt                                                                 >  2:120848/1‑100 (MQ=255)
                                aCCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTg                                                               >  1:1273643/1‑100 (MQ=255)
                                aCCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTg                                                               >  1:1407018/1‑100 (MQ=255)
                                 ccAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGc                                                              >  2:4064410/1‑100 (MQ=255)
                                  cAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCa                                                             >  2:5243053/1‑100 (MQ=255)
                                   aaaTAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCAc                                                            <  2:2262858/100‑1 (MQ=255)
                                      taatTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTgg                                                         >  1:2620126/1‑100 (MQ=255)
                                       aatTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGa                                                        <  2:4612298/100‑1 (MQ=255)
                                        atTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGaa                                                       >  1:4609688/1‑100 (MQ=255)
                                        atTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGaa                                                       >  1:2671840/1‑100 (MQ=255)
                                         ttACGTTTGGCTAATATAGGGCAACTGCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAAc                                                      <  2:739130/100‑1 (MQ=255)
                                            cGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTc                                                   <  1:5171356/100‑1 (MQ=255)
                                            cGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTc                                                   >  2:319394/1‑100 (MQ=255)
                                            cGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTc                                                   <  1:544193/100‑1 (MQ=255)
                                             gTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCa                                                  >  1:1546974/1‑100 (MQ=255)
                                             gTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCa                                                  <  2:2328275/100‑1 (MQ=255)
                                             gTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCa                                                  <  1:2853979/100‑1 (MQ=255)
                                               ttGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAAc                                                <  2:5282890/100‑1 (MQ=255)
                                               ttGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAAc                                                <  1:4241005/100‑1 (MQ=255)
                                               ttGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAAc                                                <  1:846410/100‑1 (MQ=255)
                                                tGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACa                                               <  1:3427517/100‑1 (MQ=255)
                                                tGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACa                                               >  1:1996103/1‑100 (MQ=255)
                                                tGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACa                                               <  2:4986562/100‑1 (MQ=255)
                                                gggCTAATATAGGGCAACTCCGGGACCAGTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGCAACTTCAACa                                               >  2:5224982/2‑100 (MQ=255)
                                                 ggCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAg                                              >  2:3822695/1‑100 (MQ=255)
                                                 ggCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAg                                              <  1:3031348/100‑1 (MQ=255)
                                                 ggCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAg                                              <  1:4008981/100‑1 (MQ=255)
                                                  gCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGc                                             >  2:3531061/1‑100 (MQ=255)
                                                   cTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCt                                            >  1:370721/1‑100 (MQ=255)
                                                   cTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCt                                            <  2:4517753/100‑1 (MQ=255)
                                                   cTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCt                                            <  2:2671575/100‑1 (MQ=255)
                                                     aaTATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTa                                          >  1:4369839/1‑100 (MQ=255)
                                                       tataGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaa                                        >  1:5361617/1‑100 (MQ=255)
                                                       tataGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaa                                        >  1:5048411/1‑100 (MQ=255)
                                                        ataGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaaa                                       <  2:3354673/100‑1 (MQ=255)
                                                         taGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaaaa                                      <  1:417573/100‑1 (MQ=255)
                                                         taGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaaaa                                      >  2:4438343/1‑100 (MQ=255)
                                                           gggCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTaaaaaaa                                    >  2:10321/1‑100 (MQ=255)
                                                            ggCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAGGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAAc                                   >  2:4821313/1‑100 (MQ=255)
                                                            ggCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAAc                                   >  1:2480384/1‑100 (MQ=255)
                                                            ggCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATGTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAAc                                   >  1:1876164/1‑100 (MQ=255)
                                                             gCAACTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAAcc                                  <  2:84636/100‑1 (MQ=255)
                                                                aCTCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTAt                               >  2:1770229/1‑100 (MQ=255)
                                                                  tCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATcc                             >  1:439934/1‑100 (MQ=255)
                                                                  tCCGGGACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTAAAAAAAAAAAAAAAAACCTATcc                             <  2:2641066/100‑1 (MQ=255)
                                                                      ggACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggn                         <  1:1146198/100‑2 (MQ=255)
                                                                      ggACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggc                         <  1:2367090/100‑1 (MQ=255)
                                                                      ggACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggc                         >  1:2744108/1‑100 (MQ=255)
                                                                      ggACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggc                         >  1:577152/1‑100 (MQ=255)
                                                                      ggACGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggc                         <  2:1658798/100‑1 (MQ=255)
                                                                         cGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggcggc                      <  2:2984272/100‑1 (MQ=255)
                                                                         cGATTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAggcggc                      <  2:2776201/100‑1 (MQ=255)
                                                                           aTTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAATAAACCTATCCAGGCGGCGt                    >  2:4914778/1‑100 (MQ=255)
                                                                           aTTTAAACCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGt                    <  1:1740305/100‑1 (MQ=255)
                                                                             ttcaaCCACAGATAAAGACACAGAAGGTAAAGTAAGTAAAAATTTATGACCATTTTACTTTAACGTCAACAGCTTAAAAAAACCTATTCATGCGGCGTTc                  <  2:329754/97‑1 (MQ=255)
                                                                               aaaCCACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAg                <  2:429153/100‑1 (MQ=255)
                                                                                  ccACAGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGgcg             >  1:3395518/1‑100 (MQ=255)
                                                                                   cacaGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCCGGCGGCGTTCAGgcgc            >  2:3089110/1‑100 (MQ=255)
                                                                                   cacaGATAAAGTGTCAGAACATAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGgcgc            >  1:4342746/1‑100 (MQ=255)
                                                                                      aGATAAAGTGTCAGAACGTAAGGTAAGTAAATATTTATGTCCATTGCACTGGAACTTCAACAGCTTAAAAGAACCTTTCCGGGCGGCGTTCAGGCGCTTc         <  2:541757/100‑1 (MQ=255)
                                                                                      aGATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTc         <  2:1255351/100‑1 (MQ=255)
                                                                                       gATAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCg        <  1:4095572/100‑1 (MQ=255)
                                                                                        aTAAAGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGt       <  1:3106467/100‑1 (MQ=255)
                                                                                           aaGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGTggg    >  2:2832467/1‑100 (MQ=255)
                                                                                           aaGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGTggg    <  1:3852589/100‑1 (MQ=255)
                                                                                            aGTGTCAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGTGGGa   <  2:2594555/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGCTAATATAGGGCAACTCCGGGACGATTTAAACCACAGATAAAGTGTAAGAACGTAAGGTAAGTAAAAATTTATGACCATTGCACTGGAACTTCAACAGCTTAAAAAAACCTATCCAGGCGGCGTTCAGGCGCTTCGTGGGAT  >  NC_000913/142657‑142849

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