Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,821 T→G 20.0% intergenic (+70/‑31) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,8210TG20.0% 204.0 / 52.8 99intergenic (+70/‑31)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base T (33/46);  new base G (13/7);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 8.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 102 total aligned reads displayed.

GTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGG  >  NC_000913/1689723‑1689920
                                                                                                  |                                                                                                   
gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTc                                                                                                    >  1:351585/1‑100 (MQ=255)
 tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTcc                                                                                                   <  1:1057111/100‑1 (MQ=255)
   tAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGc                                                                                                 <  2:1855118/100‑1 (MQ=255)
    aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                                <  1:1106362/100‑1 (MQ=255)
    aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                                <  2:699793/100‑1 (MQ=255)
    aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                                >  2:697230/1‑100 (MQ=255)
     gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCt                                                                                               >  1:2177479/1‑100 (MQ=255)
      cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                              <  1:475736/100‑1 (MQ=255)
       gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATTACACATTCAGGCGGCCAGTCCg                                                                                                  <  1:1097957/95‑1 (MQ=255)
       gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTg                                                                                              <  1:15148/99‑1 (MQ=255)
        cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAt                                                                                            <  2:2407981/100‑1 (MQ=255)
        cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaa                                                                                            <  2:1969903/100‑2 (MQ=255)
        cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaa                                                                                            <  2:855783/100‑2 (MQ=255)
          tgtTTACAACAAGCTGTAAGAGCTTTCTGAACAACATAACCTCTCTTGCTAATCTTATTAAAGGCATATAACACCTTCAGGCGTACTGGCCGCCTGaagg                                                                                          <  1:2483875/100‑4 (MQ=255)
          tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                          <  2:322785/100‑4 (MQ=255)
            tttACAATAAGCTGTAAGAACTTATGAAAAAAATTAACATCTCTTGCTATGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                          <  1:2322252/98‑4 (MQ=255)
             ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAt                                                                                       >  1:800952/1‑100 (MQ=255)
               acaacaAGCTGTAATAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAttt                                                                                     <  1:1763780/100‑1 (MQ=255)
               acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTTCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAttt                                                                                     >  1:745028/1‑100 (MQ=255)
                caacaaGCTTTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                    <  2:710951/100‑1 (MQ=255)
                caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                    >  2:200851/1‑100 (MQ=255)
                   caaGCTGTAAGTGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACATCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                  >  1:1034739/1‑99 (MQ=255)
                   caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                                 <  2:301098/100‑1 (MQ=255)
                   caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                                 <  2:284459/100‑1 (MQ=255)
                   caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                                 >  2:724782/1‑100 (MQ=255)
                    aaGCTGTAAGATCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATgg                                                                                >  1:1354805/1‑100 (MQ=255)
                    aaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATgg                                                                                <  1:2293279/100‑1 (MQ=255)
                      gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                              <  2:370061/100‑1 (MQ=255)
                           aaGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCa                                                                         >  2:804650/1‑100 (MQ=255)
                            agagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                        >  1:1960666/1‑100 (MQ=255)
                             gagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                       <  2:815686/100‑1 (MQ=255)
                              agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                      <  2:50238/100‑1 (MQ=255)
                              agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                      <  1:1995455/100‑1 (MQ=255)
                               gCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                     >  1:447073/1‑100 (MQ=255)
                                 ttACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGa                                                                   >  1:1449500/1‑100 (MQ=255)
                                    cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                                <  2:1486551/100‑1 (MQ=255)
                                      gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                              >  2:451862/1‑100 (MQ=255)
                                      gAAAAAATTAACATCTCTTGCTAAGCTGATTACAGGCTTATAACACCTTCATGCAGCCAGTCCGCGTGATTTAATTTTATGGATAAACATTATGAATa                                                                >  1:315175/1‑98 (MQ=255)
                                         aaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTGCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                            <  1:1514067/99‑1 (MQ=255)
                                          aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                          <  2:826298/100‑1 (MQ=255)
                                              tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGt                                                      >  2:1885395/1‑100 (MQ=255)
                                                       ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTATCGGTAgg                                             <  2:1354813/100‑1 (MQ=255)
                                                       ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                             >  2:2441890/1‑100 (MQ=255)
                                                         gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCg                                           <  2:771578/100‑1 (MQ=255)
                                                           tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTc                                         <  2:1466782/100‑1 (MQ=255)
                                                             aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                       >  2:2013240/1‑100 (MQ=255)
                                                              gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                      <  1:200849/100‑1 (MQ=255)
                                                              gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                      <  2:2389669/100‑1 (MQ=255)
                                                               cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                     <  1:65449/100‑1 (MQ=255)
                                                                ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTTAttgt                                    >  2:2320623/1‑100 (MQ=255)
                                                                ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgt                                    >  1:198202/1‑100 (MQ=255)
                                                                 tattaAAGGCTTATAACACCTTCAGGCGGCCAGTTCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCCttgtt                                   >  1:230021/1‑100 (MQ=255)
                                                                  attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                                  <  1:1757196/100‑1 (MQ=255)
                                                                  attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAGTTCATGTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                                  <  2:635267/100‑1 (MQ=255)
                                                                   ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                                 <  1:2166846/100‑1 (MQ=255)
                                                                   ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                                 <  2:1449508/100‑1 (MQ=255)
                                                                     aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                               >  2:2477976/1‑100 (MQ=255)
                                                                        ggCTGATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCGGGTAGCGGTAGGCGTCATTGTTGCGCTAg                            <  2:728849/100‑1 (MQ=255)
                                                                         gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                           >  2:304435/1‑100 (MQ=255)
                                                                         gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                           <  2:1179474/100‑1 (MQ=255)
                                                                         gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                           <  1:590111/100‑1 (MQ=255)
                                                                          cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                          >  1:1310146/1‑100 (MQ=255)
                                                                          cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                          <  2:2177493/100‑1 (MQ=255)
                                                                           ttATAACACCTTCATGCGGCCAGTCCGCCTGATTGCATTTTATCGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                         <  1:1467164/100‑1 (MQ=255)
                                                                           ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                         >  1:1194024/1‑100 (MQ=255)
                                                                           ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                         <  1:1918579/100‑1 (MQ=255)
                                                                            tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGt                        <  1:2216596/100‑1 (MQ=255)
                                                                              aaatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                      >  2:2452301/9‑100 (MQ=38)
                                                                                acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                    <  2:1861286/100‑1 (MQ=255)
                                                                                 tgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATAGATATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                  >  1:1391380/6‑100 (MQ=37)
                                                                                  gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                  >  2:204575/5‑100 (MQ=255)
                                                                                  gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                  >  2:1791556/5‑100 (MQ=255)
                                                                                  gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                  >  1:1939553/5‑100 (MQ=255)
                                                                                  gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                  >  1:533612/5‑100 (MQ=255)
                                                                                   ttccgatctggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                 <  2:274049/91‑1 (MQ=255)
                                                                                   ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTTTgg                 >  2:627611/1‑100 (MQ=255)
                                                                                   ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                 <  1:304432/100‑1 (MQ=255)
                                                                                   aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                 >  2:2440594/4‑100 (MQ=255)
                                                                                   aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                 >  2:62764/4‑100 (MQ=255)
                                                                                   aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                 >  2:1059158/4‑100 (MQ=255)
                                                                                    ccgatctcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                <  1:1728078/93‑1 (MQ=255)
                                                                                    aatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                >  2:2046366/3‑100 (MQ=255)
                                                                                     atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATTGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc               >  1:1022063/2‑100 (MQ=255)
                                                                                     atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc               >  2:1585540/2‑100 (MQ=255)
                                                                                      tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa              <  2:1813570/100‑1 (MQ=255)
                                                                                      tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa              >  2:74199/1‑100 (MQ=255)
                                                                                        aGGCGGACTGGCCGCCTGATTTCATTTTCTGGGTAATCATTATGAATAAATCGCCGGTAGCGGTAGGCGTCATTGTTGCGCTAGCCTTAGTCTGGAAAg             >  1:842223/1‑99 (MQ=255)
                                                                                         ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc           >  1:2436207/1‑100 (MQ=255)
                                                                                         ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGATAGGCGTAGACTGGACAggc           >  1:70157/1‑100 (MQ=255)
                                                                                           cggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGaagatcgg         >  2:1728087/1‑93 (MQ=255)
                                                                                            ggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGagatcggaa        >  1:274046/1‑92 (MQ=255)
                                                                                            ggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc        <  2:415818/100‑1 (MQ=255)
                                                                                             gcCAGTCCGACTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg       <  1:1228142/100‑1 (MQ=255)
                                                                                              ccAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc      >  2:2215714/1‑100 (MQ=255)
                                                                                                aGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt    <  2:1385650/100‑1 (MQ=255)
                                                                                                aGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAAAGGCGGCGCAt    >  1:697340/1‑100 (MQ=255)
                                                                                                 gTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATg   <  1:74198/100‑1 (MQ=255)
                                                                                                  tCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGTGCTAGGCGTAGTCTGGACAGgcggcg       >  1:2189091/1‑95 (MQ=255)
                                                                                                  tCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg  <  2:135751/100‑1 (MQ=255)
                                                                                                  tCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg  <  2:1249468/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGG  >  NC_000913/1689723‑1689920

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: