Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,968,795 C→T R364C (CGC→TGC)  wecB → UDP‑N‑acetyl glucosamine ‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,968,7950CT85.8% 924.2 / 109.1 317R364C (CGC→TGC) wecBUDP‑N‑acetyl glucosamine ‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (138/134);  total (162/155)
Fisher's exact test for biased strand distribution p-value = 8.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 323 total aligned reads displayed.

ACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCA  >  NC_000913/3968705‑3968887
                                                                                          |                                                                                            
  ggTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTgcgc                                                                                          <  1:4837324/93‑1 (MQ=255)
   gTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGcc                                                                                         <  1:2826254/93‑1 (MQ=255)
    tGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc                                                                                        <  1:776692/93‑1 (MQ=255)
    tGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc                                                                                        <  1:3884335/93‑1 (MQ=255)
    ggCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGccc                                                                                        <  1:1481980/92‑1 (MQ=255)
      cGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCAt                                                                                      >  1:2060144/1‑93 (MQ=255)
       gTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATa                                                                                     <  1:2564886/93‑1 (MQ=255)
        tCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATaa                                                                                    >  2:1536948/1‑93 (MQ=255)
        tCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATaa                                                                                    >  1:954736/1‑93 (MQ=255)
         cTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAAc                                                                                   <  2:2227319/93‑1 (MQ=255)
         cTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAAc                                                                                   <  2:3873830/93‑1 (MQ=255)
          tGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAAcc                                                                                  >  1:1269613/1‑93 (MQ=255)
            gTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCg                                                                                <  1:4958686/93‑1 (MQ=255)
            gTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCg                                                                                <  1:5857325/93‑1 (MQ=255)
            gTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCg                                                                                <  2:4169138/93‑1 (MQ=255)
            gTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCg                                                                                <  1:519273/93‑1 (MQ=255)
              aGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTa                                                                              <  2:4872813/93‑1 (MQ=255)
               ggCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTAt                                                                             >  1:6268639/1‑93 (MQ=255)
               ggCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTAt                                                                             >  1:1000367/1‑93 (MQ=255)
               ggCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTAt                                                                             <  1:411028/93‑1 (MQ=255)
                gCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATg                                                                            >  2:6436292/1‑93 (MQ=255)
                 cACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATgg                                                                           <  1:425537/93‑1 (MQ=255)
                 cACGGATAAGCAGCGAAGTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATgg                                                                           <  2:5287456/93‑1 (MQ=255)
                  aCGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATggg                                                                          >  1:6393496/1‑92 (MQ=255)
                      aTAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATg                                                                      >  1:1956649/1‑93 (MQ=255)
                       tAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATgg                                                                     >  1:746331/1‑93 (MQ=255)
                         agcagcGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTc                                                                   <  2:5231814/93‑1 (MQ=255)
                          gcagcGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCa                                                                  <  2:1952667/93‑1 (MQ=255)
                           cagcGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAg                                                                 <  2:2964112/93‑1 (MQ=255)
                           cagcGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAg                                                                 <  2:3568514/93‑1 (MQ=255)
                            agcGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAgg                                                                >  2:5447341/1‑93 (MQ=255)
                              cGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCa                                                              <  1:3988689/93‑1 (MQ=255)
                                  ttGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGct                                                          >  2:4996000/1‑93 (MQ=255)
                                   tGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctc                                                         >  1:250655/1‑93 (MQ=255)
                                    gTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctct                                                        <  2:4439969/93‑1 (MQ=255)
                                    gTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctct                                                        <  2:5015077/93‑1 (MQ=255)
                                    gTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctct                                                        <  1:5664182/93‑1 (MQ=255)
                                    gTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctct                                                        <  2:1624873/93‑1 (MQ=255)
                                    gTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGctct                                                        >  2:5790414/1‑93 (MQ=255)
                                     tCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCtctc                                                       >  1:3711229/1‑93 (MQ=255)
                                     tCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCtctc                                                       >  1:2753931/1‑93 (MQ=255)
                                     tCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCtctc                                                       >  2:3991045/1‑93 (MQ=255)
                                       gAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGc                                                     <  1:2412258/93‑1 (MQ=255)
                                        aGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCa                                                    <  1:5516188/93‑1 (MQ=255)
                                          gAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCAtt                                                  <  2:3546524/93‑1 (MQ=255)
                                           aaGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTc                                                 >  1:5103120/1‑93 (MQ=255)
                                           aaGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTc                                                 <  2:2427011/93‑1 (MQ=255)
                                           aaGTGACGCGTCTTTTAAAAGACAAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTc                                                 <  2:5466468/93‑1 (MQ=255)
                                            aGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCt                                                >  1:5598787/1‑93 (MQ=255)
                                            aGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCt                                                <  1:2261791/93‑1 (MQ=255)
                                             gTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTg                                               >  2:1897675/1‑93 (MQ=255)
                                              tGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTgg                                              <  1:5430854/93‑1 (MQ=255)
                                               gACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGa                                             >  2:5651259/1‑93 (MQ=255)
                                                aCGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGaa                                            <  2:5716035/93‑1 (MQ=255)
                                                aCGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGaa                                            >  2:4580199/1‑93 (MQ=255)
                                                 cgcgTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAg                                           <  2:796187/93‑1 (MQ=255)
                                                  gcgTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGc                                          >  2:3159224/1‑93 (MQ=255)
                                                  gcgTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGc                                          <  2:1409972/93‑1 (MQ=255)
                                                      cTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTa                                      <  1:2678603/93‑1 (MQ=255)
                                                       ttttAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATACTCTCGCATTCTGGAAGCGTTaa                                     >  1:6621836/1‑93 (MQ=255)
                                                        tttAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTaaa                                    >  2:4821764/1‑93 (MQ=255)
                                                        tttAAAAGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTaaa                                    >  2:1103277/1‑93 (MQ=255)
                                                         ttAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTaaaa                                   <  2:702927/93‑1 (MQ=255)
                                                             aaGACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAaata                               >  1:2894703/1‑93 (MQ=255)
                                                              aGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAaataa                              >  1:6509710/1‑93 (MQ=255)
                                                               gACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAaataat                             <  2:1737226/93‑1 (MQ=255)
                                                               gACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAaataat                             >  2:3682658/1‑93 (MQ=255)
                                                               gACGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAaataat                             >  2:5177142/1‑93 (MQ=255)
                                                                aCGAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATc                            <  2:6054668/93‑1 (MQ=255)
                                                                  gAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCgg                          <  2:1248974/93‑1 (MQ=255)
                                                                  gAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCgg                          <  1:865011/93‑1 (MQ=255)
                                                                  gAAAACGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCgg                          <  2:4037398/93‑1 (MQ=255)
                                                                  gAAAACGAATATCAAGATATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCgg                          >  1:1661777/1‑93 (MQ=255)
                                                                   aaaaCGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGa                         >  2:4379579/1‑93 (MQ=255)
                                                                    aaaCGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGat                        <  1:1092839/93‑1 (MQ=255)
                                                                    aaaCGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGat                        >  1:5318272/1‑93 (MQ=255)
                                                                     aaCGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGata                       >  1:187558/1‑93 (MQ=255)
                                                                      aCGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGatat                      <  2:287883/93‑1 (MQ=255)
                                                                       cGAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATc                     <  1:3533874/93‑1 (MQ=255)
                                                                        gAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCa                    <  2:263624/93‑1 (MQ=255)
                                                                        gAATATCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCa                    <  1:59507/93‑1 (MQ=255)
                                                                        gAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCa                    >  2:6229072/1‑93 (MQ=255)
                                                                         aaTATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCAc                   >  2:2052653/1‑93 (MQ=255)
                                                                           tatCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTa                 >  2:4495589/1‑93 (MQ=255)
                                                                            atCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTAt                <  2:3414072/93‑1 (MQ=255)
                                                                             tCAAGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATg               >  1:114475/1‑93 (MQ=255)
                                                                               aaGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAg             >  1:3088628/1‑93 (MQ=255)
                                                                               aaGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAg             >  2:305951/1‑93 (MQ=255)
                                                                               aaGCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAg             >  2:4300107/1‑93 (MQ=255)
                                                                                 gCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtt           <  1:3477118/93‑1 (MQ=255)
                                                                                 gCTATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtt           >  2:473409/1‑93 (MQ=255)
                                                                                 gCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtt           >  1:4539208/1‑93 (MQ=255)
                                                                                   tATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtttt         >  2:5990150/1‑93 (MQ=255)
                                                                                   tATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtttt         >  2:6136035/1‑93 (MQ=255)
                                                                                   tATGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGtttt         >  2:6569450/1‑93 (MQ=255)
                                                                                     tGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGc       <  1:3613656/93‑1 (MQ=255)
                                                                                     tGAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGc       >  2:1809226/1‑93 (MQ=255)
                                                                                      gAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGTGTTTTGCg      <  1:3248927/93‑1 (MQ=255)
                                                                                      gAGCTGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCg      >  1:4776186/1‑93 (MQ=255)
                                                                                          cgcgCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCa  <  2:5073936/93‑1 (MQ=255)
                                                                                          |                                                                                            
ACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACCCGTATGGTGATGGTCAGGCATGCTCTCGCATTCTGGAAGCGTTAAAAAATAATCGGATATCACTATGAGTTTTGCGACCA  >  NC_000913/3968705‑3968887

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