Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,819 A→T 20.8% intergenic (+68/‑33) 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,8190AT20.8% 204.7 / 51.6 96intergenic (+68/‑33)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base A (32/44);  new base T (13/7);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 8.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAT  >  NC_000913/1689720‑1689918
                                                                                                   |                                                                                                   
gCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCa                                                                                                     >  1:1302435/1‑100 (MQ=255)
gCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCa                                                                                                     >  1:704207/1‑100 (MQ=255)
   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:933860/100‑2 (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                         <  1:590111/100‑1 (MQ=255)
                                                                            gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         >  2:304435/1‑100 (MQ=255)
                                                                            gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         <  2:1179474/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:1918579/100‑1 (MQ=255)
                                                                              ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                       >  1:1194024/1‑100 (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                >  1:533612/5‑100 (MQ=255)
                                                                                     gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                >  2:1791556/5‑100 (MQ=255)
                                                                                     gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                >  1:1939553/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:1059158/4‑100 (MQ=255)
                                                                                      aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg               >  2:62764/4‑100 (MQ=255)
                                                                                       ccgatctcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa              <  1:1728078/93‑1 (MQ=255)
                                                                                       aatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa              >  2:2046366/3‑100 (MQ=255)
                                                                                        ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc             <  2:1302442/100‑1 (MQ=255)
                                                                                        atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATTGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc             >  1:1022063/2‑100 (MQ=255)
                                                                                        atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc             >  2:1585540/2‑100 (MQ=255)
                                                                                         tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa            >  2:74199/1‑100 (MQ=255)
                                                                                         tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa            <  2:1813570/100‑1 (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)
                                                                                                   |                                                                                                   
GCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAT  >  NC_000913/1689720‑1689918

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: