Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,247 +C intergenic (+10/+29) yghB → / ← yqhC hypothetical protein/putative ARAC‑type regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,2441.C99.5% 689.7 / ‑4.0 212intergenic (+7/+32)yghB/yqhChypothetical protein/putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (127/84);  total (128/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 257 total aligned reads displayed.

TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAG  >  NC_000913/3152152‑3152333
                                                                                              |                                                                                         
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                            >  1:5270788/1‑93 (MQ=255)
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                            >  2:3884367/1‑93 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          <  2:424253/93‑1 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                        <  2:1106554/93‑1 (MQ=255)
   tgatCCTGCAAATTGCTTTTTTAACCGCTGGTTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                        <  2:360989/93‑1 (MQ=255)
    gatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcccc                                                                                       >  2:5953456/1‑92 (MQ=255)
      tCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgt                                                                                     >  1:3558631/1‑90 (MQ=255)
      tCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgt                                                                                     >  1:2553317/1‑90 (MQ=255)
      tCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgt                                                                                     >  2:2450874/1‑90 (MQ=255)
       ccTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgtt                                                                                    >  2:1527132/1‑89 (MQ=255)
       ccTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccgtt                                                                                    >  2:1824034/1‑89 (MQ=255)
          gCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtccccgttccc                                                                                  <  2:4285621/93‑1 (MQ=255)
            cAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGGTTGTGGGGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCgg                                                                                >  2:3640438/1‑93 (MQ=255)
             aaTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGt                                                                               <  2:5118022/93‑1 (MQ=255)
             aaTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGt                                                                               >  1:2127752/1‑93 (MQ=255)
             aaTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGt                                                                               <  1:28577/93‑1 (MQ=255)
              aTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGtt                                                                              >  1:441318/1‑93 (MQ=255)
                tGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGATGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGc                                                                            <  2:6262009/93‑1 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                                                           <  1:5839707/93‑1 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                                                           >  1:1971919/1‑93 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGGTCCCCGTTGCt                                                                           >  2:261540/1‑93 (MQ=255)
                  ccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCCGctgctg                                                                          >  1:2678827/1‑93 (MQ=255)
                   cTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGt                                                                         >  1:644492/1‑93 (MQ=255)
                    ttGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTa                                                                        >  2:2800641/1‑93 (MQ=255)
                    ttGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTa                                                                        <  2:4796000/93‑1 (MQ=255)
                       ttAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCg                                                                     <  1:3120504/93‑1 (MQ=255)
                       ttAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCg                                                                     >  1:6019384/1‑93 (MQ=255)
                       ttAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCg                                                                     <  1:4097685/93‑1 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    <  2:188907/93‑1 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    >  2:3497890/1‑93 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    <  2:6129426/93‑1 (MQ=255)
                         aaCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCggg                                                                   <  2:2108128/93‑1 (MQ=255)
                         aaCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCggg                                                                   >  1:4996067/1‑93 (MQ=255)
                          aCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                                                                  >  2:2926752/1‑93 (MQ=255)
                            cGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAAc                                                                >  1:6259585/1‑93 (MQ=255)
                             gCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACg                                                               >  1:2242600/1‑93 (MQ=255)
                             gCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACg                                                               >  2:942718/1‑93 (MQ=255)
                               tGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTa                                                             >  2:3149470/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAt                                                            >  1:5989043/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATAATGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTGTACCGGGAAAGTAt                                                            >  2:3957457/1‑93 (MQ=255)
                                 gCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAtt                                                           >  1:5205007/1‑93 (MQ=255)
                                    tGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTaa                                                        >  1:2448658/1‑93 (MQ=255)
                                     gTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAAt                                                       >  1:2111610/1‑93 (MQ=255)
                                     gTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGCTGCTTTACCGGGAACGTCTTTAAt                                                       >  2:2354804/1‑93 (MQ=255)
                                       tAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTc                                                     <  2:2652577/93‑1 (MQ=255)
                                         ggCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCCGTTGCTGTACCGGGGACGTATTTTATTccc                                                   >  1:1118021/1‑93 (MQ=255)
                                          gCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTcccc                                                  <  1:909874/93‑1 (MQ=255)
                                          gCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTcccc                                                  <  1:2410162/93‑1 (MQ=255)
                                           cACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCt                                                 >  1:1532449/1‑93 (MQ=255)
                                           cACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCt                                                 <  1:1424451/93‑1 (MQ=255)
                                            aCGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTg                                                >  1:5460604/1‑93 (MQ=255)
                                            aCGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGGTCCCGGGTGCTGTACCCGGAACGGATTTTATTCCCCtt                                                >  1:2302543/1‑92 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                              <  2:769354/93‑1 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                              <  2:5198969/93‑1 (MQ=255)
                                               cTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                             >  2:2013758/1‑93 (MQ=255)
                                               cTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                             >  1:323786/1‑93 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATATTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                            >  1:3066363/1‑93 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                            >  2:4083868/1‑93 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                          >  1:3536225/1‑93 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                          <  2:2148925/93‑1 (MQ=255)
                                                   ttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGcc                                         >  1:3429403/1‑93 (MQ=255)
                                                    tgtgATGAATAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGccc                                        <  1:4351373/93‑1 (MQ=255)
                                                        gtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAt                                    >  1:6177034/1‑93 (MQ=255)
                                                        gtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAt                                    <  2:5541784/93‑1 (MQ=255)
                                                         tgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCAtt                                   >  2:5367961/1‑93 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGGTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                 >  1:184700/1‑93 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCtt                                <  2:5709736/93‑1 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCtt                                <  2:3419961/93‑1 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCtt                                <  2:6096892/93‑1 (MQ=255)
                                                             tAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTg                               >  2:2365442/1‑93 (MQ=255)
                                                             tAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTg                               >  1:5678522/1‑93 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGGATTTAATTCCCCTGCCTCGTCCGCCTTCTTGc                              >  2:5766256/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                             <  1:3879609/93‑1 (MQ=255)
                                                                  aaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCa                          <  1:562413/93‑1 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                        <  1:2039448/93‑1 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                        >  1:1133952/1‑93 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                        >  1:4711058/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGGACCGGGAACGTATTTAATTCCACTGCCTCGCCCGCATTCTTGCCGCATCTTc                     >  2:657140/1‑93 (MQ=255)
                                                                        aCTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcc                    <  2:494109/93‑1 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                  >  2:2864974/1‑93 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                  <  1:5425231/93‑1 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                  <  2:123936/93‑1 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTGTACCGGGGACGGATTTAATTTCCCTGCATCGCCCCCACTCTTTCCCCCTCTTcccc                  >  2:3278609/1‑93 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                 <  2:3256319/93‑1 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                 >  2:4624822/1‑93 (MQ=255)
                                                                            tAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCg                >  2:1939537/1‑93 (MQ=255)
                                                                            tAACGCCTGACGATTTTCCCCGTTCCCCGGTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTTCCCCATCTTcccccc                >  2:3333880/1‑92 (MQ=255)
                                                                              aCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGc              >  1:4901759/1‑93 (MQ=255)
                                                                                gCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGt            >  2:4754512/1‑93 (MQ=255)
                                                                                 ccTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTc           <  1:4221655/93‑1 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGGATTTTATTTCCCTGCATCGCCCGCCTTCTTTCCGCCTCTTTCCCCGGCGTcac         >  2:5484161/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTTATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTcac         >  2:197651/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac         >  2:781805/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac         >  1:5047344/1‑93 (MQ=255)
                                                                                    gACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcaca        >  1:3991097/1‑93 (MQ=255)
                                                                                    gACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTcaca        >  1:3863569/1‑93 (MQ=255)
                                                                                     aCGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTTCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCacac       >  1:1381876/1‑93 (MQ=255)
                                                                                      cGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc      <  1:3212361/93‑1 (MQ=255)
                                                                                        aTTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa    >  1:4818903/1‑93 (MQ=255)
                                                                                         tttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAg  <  1:1743975/93‑1 (MQ=255)
                                                                                              |                                                                                         
TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAG  >  NC_000913/3152152‑3152333

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