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.C98.8% 1371.1 / 5.1 423intergenic (+7/+32)yghB/yqhChypothetical protein/putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (244/173);  total (246/177)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 525 total aligned reads displayed.

TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAG  >  NC_000913/3152152‑3152333
                                                                                               |                                                                                         
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTTTAACGCCCGACGAtttt                                                                                            >  2:6051865/1‑93 (MQ=255)
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                            >  2:3847479/1‑93 (MQ=255)
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGCTTGTGGTGATTAAAAAAAAATACTGTAACCGCTGACGAtttt                                                                                            >  1:640520/1‑93 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          <  1:4000544/93‑1 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          >  1:3591738/1‑93 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          <  1:6159520/93‑1 (MQ=255)
  atgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcc                                                                                         >  1:8238507/1‑93 (MQ=255)
  atgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcc                                                                                         >  2:5113487/1‑93 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                        >  2:4703087/1‑93 (MQ=255)
    gatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcccc                                                                                       >  1:4897708/1‑92 (MQ=255)
    gatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcccc                                                                                       <  2:8388699/93‑2 (MQ=255)
     atCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccccg                                                                                      >  2:7930650/1‑91 (MQ=255)
         tGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtccccgttcc                                                                                   <  2:809858/93‑1 (MQ=255)
          gCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtccccgttccc                                                                                  <  1:2934106/93‑1 (MQ=255)
               ttGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTg                                                                             <  2:8846693/93‑1 (MQ=255)
               ttGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTg                                                                             >  2:3068573/1‑93 (MQ=255)
               ttGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGGTg                                                                             >  2:112292/1‑93 (MQ=255)
               ttGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCCGTTg                                                                             >  2:141863/1‑93 (MQ=255)
                  ccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTg                                                                          >  2:1605138/1‑93 (MQ=255)
                  ccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTg                                                                          >  2:7055568/1‑93 (MQ=255)
                   cTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGt                                                                         <  1:1354334/93‑1 (MQ=255)
                    ttGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTTCCGGTTGCTGTa                                                                        >  1:5542631/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       <  1:652216/93‑1 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       >  2:7332035/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       <  2:3670739/93‑1 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       <  1:5210708/93‑1 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       >  2:9591536/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                       <  1:575622/93‑1 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGCTGCTGTAc                                                                       >  2:9073983/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTCCCCGGTTGCTGTAc                                                                       >  1:1826650/1‑93 (MQ=255)
                      gTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTCTCCCCGCTCCCGGTTGCTGTAcc                                                                      >  1:4122210/1‑93 (MQ=255)
                      gTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAcc                                                                      >  1:2287923/1‑93 (MQ=255)
                      gTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAcc                                                                      >  1:7915592/1‑93 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    >  2:2665933/1‑93 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    >  2:5239201/1‑93 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                    <  2:6861874/93‑1 (MQ=255)
                         aaCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGGACCggg                                                                   >  1:2137494/1‑93 (MQ=255)
                           ccGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCCGGaa                                                                 >  2:6015302/1‑93 (MQ=255)
                           ccGCTGGCTTATTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGaa                                                                 <  1:8885466/93‑1 (MQ=255)
                            cGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAAc                                                                >  1:7654736/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTAt                                                            >  2:9771050/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAt                                                            >  1:4305813/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAt                                                            <  1:5961922/93‑1 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAt                                                            >  2:1318797/1‑93 (MQ=255)
                                   ttGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTa                                                         >  1:4722980/1‑93 (MQ=255)
                                    tGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTaa                                                        <  2:6474337/93‑1 (MQ=255)
                                    tGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTaa                                                        >  2:9644353/1‑93 (MQ=255)
                                         ggCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTtcc                                                   >  1:2582185/1‑93 (MQ=255)
                                            aCGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTg                                                >  1:4517905/1‑93 (MQ=255)
                                            aCGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTg                                                >  1:9922103/1‑93 (MQ=255)
                                             cGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGc                                               >  2:574148/1‑93 (MQ=255)
                                               cTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                             >  2:928142/1‑93 (MQ=255)
                                               cTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCAt                                             <  1:9343594/93‑1 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                          >  1:7070906/1‑93 (MQ=255)
                                                   ttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGcc                                         <  2:5488205/93‑1 (MQ=255)
                                                     gtggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCg                                       >  1:10017530/1‑93 (MQ=255)
                                                      tggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGc                                      >  1:7256127/1‑93 (MQ=255)
                                                       gcgtaTTAAAAAAAAATACTGTAACGCCTGACGATGTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCa                                     <  1:6876131/89‑1 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                 <  1:5675598/93‑1 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                 >  2:6564044/1‑93 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACAGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                 >  2:8502987/1‑93 (MQ=255)
                                                            ttAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCtt                                >  2:4829608/1‑93 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                              <  1:2871159/93‑1 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                              <  1:7738982/93‑1 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                              <  1:8681616/93‑1 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                              >  2:2909867/1‑93 (MQ=255)
                                                                aaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTCCATCGCCCGCATTCTTGccc                            >  2:7586512/1‑92 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                           >  2:4178031/1‑93 (MQ=255)
                                                                     aaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCt                       <  2:4995011/93‑1 (MQ=255)
                                                                     aaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCt                       >  2:6378976/1‑93 (MQ=255)
                                                                      aTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCtt                      >  2:1797595/1‑93 (MQ=255)
                                                                      aTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCtt                      >  1:961365/1‑93 (MQ=255)
                                                                      aTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCtt                      >  1:2443584/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                     >  1:5893637/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                     >  2:6213440/1‑93 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTCCCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                 <  2:2030055/93‑1 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                 <  1:2012613/93‑1 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                 >  1:3910619/1‑93 (MQ=255)
                                                                            tAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCg                >  2:8728846/1‑93 (MQ=255)
                                                                            tAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCg                >  1:2790732/1‑93 (MQ=255)
                                                                             aaCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCgg               >  2:6665665/1‑93 (MQ=255)
                                                                             aaCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCgg               >  2:3423401/1‑93 (MQ=255)
                                                                              aCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGc              <  1:7544451/93‑1 (MQ=255)
                                                                              aCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGc              <  1:1133566/93‑1 (MQ=255)
                                                                                gCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGt            <  2:9507453/93‑1 (MQ=255)
                                                                                 ccTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTc           >  1:8074358/1‑93 (MQ=255)
                                                                                 ccTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTc           >  1:3103641/1‑93 (MQ=255)
                                                                                  cTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTTCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTca          >  1:3798833/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTTCCCCAGCGTcac         >  2:4461936/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac         >  1:3396049/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac         >  2:4944734/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac         >  2:1722731/1‑93 (MQ=255)
                                                                                     aCGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCacac       >  1:1206076/1‑93 (MQ=255)
                                                                                     aCGATTTTCCCCGTTCCCGGGTGCTGTACCGGGGAGGGATTTAATTCCCCTGCATCGCCCCCCTTCTTGCCGCATCTTCCCCCCGCGTCccac       >  2:71716/1‑93 (MQ=255)
                                                                                      gattttcCCCGTTCCCGGTTGCTGCACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCTTCATCCCCCAGCGTCACACCg     >  1:3315234/4‑93 (MQ=255)
                                                                                      cGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc      >  1:5130922/1‑93 (MQ=255)
                                                                                      cGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc      <  2:2337308/93‑1 (MQ=255)
                                                                                       gATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTCACAccc     >  2:4038805/1‑92 (MQ=255)
                                                                                        ttttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGaa   <  1:8546206/92‑1 (MQ=255)
                                                                                         tttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAg  <  2:4286509/93‑1 (MQ=255)
                                                                                               |                                                                                         
TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAAG  >  NC_000913/3152152‑3152333

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