Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,103 2 bp→GC coding (232‑233/1263 nt) wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,1030CG99.8% 1910.3 / NA 560R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (250/309);  minor base T (1/0);  total (251/309)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,969,1040GC99.8% 1850.9 / NA 543R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (236/306);  minor base A (1/0);  total (237/306)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 578 total aligned reads displayed.

TCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC  >  NC_000913/3969011‑3969188
                                                                                            ||                                                                                    
     atCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctgg                                                                                  >  1:414363/1‑93 (MQ=255)
      tcgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggc                                                                                 >  1:2503510/1‑93 (MQ=255)
       cgtcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggct                                                                                <  2:3067536/93‑1 (MQ=255)
       cgtcgAACCTGATGTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggct                                                                                <  1:2050490/93‑1 (MQ=255)
         tcgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                                              <  1:8761162/93‑1 (MQ=255)
         tagAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                                              <  1:3740448/91‑1 (MQ=255)
          cgAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCAAGTTGAAGCGGATGCCTGGCTGAt                                                                             <  1:973117/93‑1 (MQ=255)
           gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                            >  2:7997420/1‑93 (MQ=255)
           gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                            >  2:7587122/1‑93 (MQ=255)
           gAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                            <  2:7056995/93‑1 (MQ=255)
            aaCCTGATTTGGCTAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTg                                                                           >  1:2123156/1‑93 (MQ=255)
             aCCTGATTTGGCGATTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                                          <  2:7302592/93‑1 (MQ=255)
               cTGATTTGGCGAGTGTAGTAAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                                        <  2:466330/93‑1 (MQ=255)
                 gATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                                      >  2:6805288/1‑93 (MQ=255)
                   tttGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                                    <  2:8745731/93‑1 (MQ=255)
                     tGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCa                                                                  <  2:2737655/93‑1 (MQ=255)
                      ggCGAGTGTAGTAAAAACTGCCGTAGAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 <  1:2240739/93‑1 (MQ=255)
                       gCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAcc                                                                >  1:7244799/1‑92 (MQ=255)
                       gCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                >  1:6154592/1‑93 (MQ=255)
                        cGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGc                                                               >  2:7473068/1‑93 (MQ=255)
                           gtgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGt                                                            >  2:8499737/1‑93 (MQ=255)
                             gtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                          >  1:3353511/1‑93 (MQ=255)
                             gtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                          <  2:3864023/93‑1 (MQ=255)
                             gtagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                          >  2:2089261/1‑93 (MQ=255)
                              tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTa                                                         >  1:2713097/1‑93 (MQ=255)
                               agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                                        >  2:522488/1‑93 (MQ=255)
                                gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                                       >  1:8453028/1‑93 (MQ=255)
                                 taAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGGTGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                      >  2:1783022/1‑93 (MQ=255)
                                 taAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                      <  1:8363942/93‑1 (MQ=255)
                                 taAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                      <  1:7002946/93‑1 (MQ=255)
                                 taAAAACCGCCGTAGAAGGCGCTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGTCTGATTGCTGTACCCACGCCGTTTAAgg                                                      <  1:1083365/93‑1 (MQ=255)
                                  aaaaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                     <  1:4963246/93‑1 (MQ=255)
                                  aaaaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                                     <  1:472905/93‑1 (MQ=255)
                                    aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGGTGTACCCACGCCGTTTAAGGGCg                                                   >  2:4804078/1‑93 (MQ=255)
                                    aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                   >  2:91592/1‑93 (MQ=255)
                                    aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                                   >  1:4062986/1‑93 (MQ=255)
                                     aaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgggagt                                                  >  1:1322644/1‑90 (MQ=255)
                                     aaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGa                                                  <  2:5057771/93‑1 (MQ=255)
                                       ctcccGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATc                                                <  1:7873444/90‑1 (MQ=255)
                                        tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                               <  1:1507597/93‑1 (MQ=255)
                                        tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                               <  2:2172201/93‑1 (MQ=255)
                                         gCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCAt                                              <  1:6398365/93‑1 (MQ=255)
                                         gCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCAt                                              <  1:7460094/93‑1 (MQ=255)
                                           cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                                            <  1:387773/93‑1 (MQ=255)
                                            gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCCCGCCGTTTAAGGGCGGTCATGac                                           >  2:2869815/1‑92 (MQ=255)
                                            gTAAAAGGCGTTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                                           <  1:1867149/93‑1 (MQ=255)
                                             tAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                                          >  2:371955/1‑93 (MQ=255)
                                               gAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                                        >  2:4563168/1‑93 (MQ=255)
                                                aaGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAg                                       >  1:1018684/1‑93 (MQ=255)
                                                 aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                      >  2:1091520/1‑93 (MQ=255)
                                                 aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                      <  1:4244056/93‑1 (MQ=255)
                                                 aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                      >  2:660915/1‑93 (MQ=255)
                                                 aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                      <  2:819925/93‑1 (MQ=255)
                                                 aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGa                                      >  2:7716073/1‑93 (MQ=255)
                                                  ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGCTCATGAGCCAGat                                     >  2:4367093/1‑93 (MQ=255)
                                                  ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                     >  1:6580603/1‑93 (MQ=255)
                                                  ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                     >  1:5245923/1‑93 (MQ=255)
                                                  ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                     >  2:1263324/1‑93 (MQ=255)
                                                   gCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                                    >  1:6744072/1‑93 (MQ=255)
                                                       ttttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                <  1:7768302/93‑1 (MQ=255)
                                                       ttttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                <  2:1611269/93‑1 (MQ=255)
                                                       ttttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                                >  2:834623/1‑93 (MQ=255)
                                                        tttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                               <  2:2551505/93‑1 (MQ=255)
                                                        tttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                               <  1:8581540/93‑1 (MQ=255)
                                                          tttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTa                             <  2:2396305/93‑1 (MQ=255)
                                                            tACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACg                           <  1:3641977/93‑1 (MQ=255)
                                                             aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGt                          <  2:6250858/93‑1 (MQ=255)
                                                             aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGt                          >  2:5495811/1‑93 (MQ=255)
                                                               gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                        <  1:5834578/93‑1 (MQ=255)
                                                               gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg                        >  1:6257647/1‑93 (MQ=255)
                                                                agcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa                       >  2:585919/1‑93 (MQ=255)
                                                                  cgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt                     <  1:5460502/93‑1 (MQ=255)
                                                                    agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg                   <  2:4207824/93‑1 (MQ=255)
                                                                    agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg                   >  1:310954/1‑93 (MQ=255)
                                                                    agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg                   >  2:5375343/1‑93 (MQ=255)
                                                                     gcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCgg                  >  2:1919158/1‑93 (MQ=255)
                                                                      cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                 <  1:7992860/93‑1 (MQ=255)
                                                                      cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                 <  1:5677712/93‑1 (MQ=255)
                                                                      cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgc                 <  2:5989112/93‑1 (MQ=255)
                                                                       acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  2:3579680/93‑1 (MQ=255)
                                                                       acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  1:8255487/93‑1 (MQ=255)
                                                                       acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  2:341233/93‑1 (MQ=255)
                                                                        cgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctg               >  1:4506672/1‑93 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             >  1:5612837/1‑93 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             >  1:3241805/1‑93 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             >  2:3929445/1‑93 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             <  1:482688/93‑1 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             <  2:410444/93‑1 (MQ=255)
                                                                          acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             <  2:3215310/93‑1 (MQ=255)
                                                                            gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg           >  2:1440559/1‑93 (MQ=255)
                                                                            gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg           >  2:5756433/1‑93 (MQ=255)
                                                                             ccAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgc          <  1:1367715/93‑1 (MQ=255)
                                                                               aGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgctc        <  2:2251125/93‑1 (MQ=255)
                                                                                gTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTcc       <  1:2356748/93‑1 (MQ=255)
                                                                                gTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTcc       >  1:605725/1‑93 (MQ=255)
                                                                                 ttGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGACTCGGCTGCTCGCTCCa      >  1:1697210/1‑93 (MQ=255)
                                                                                   gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAtt    <  1:7373901/93‑1 (MQ=255)
                                                                                   gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGACAGATATGACCTACGTTGAATCGGCTGCTCGCTCCAtt    >  1:4342295/1‑93 (MQ=255)
                                                                                    aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTg   <  1:2828405/93‑1 (MQ=255)
                                                                                     aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc  <  1:7565354/93‑1 (MQ=255)
                                                                                            ||                                                                                    
TCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC  >  NC_000913/3969011‑3969188

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