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.4% 545.3 / ‑4.1 165intergenic (+7/+32)yghB/yqhChypothetical protein/putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (100/64);  total (100/65)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 201 total aligned reads displayed.

TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAA  >  NC_000913/3152152‑3152332
                                                                                              |                                                                                         
tgatgaTCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGAtttt                                                                                            >  1:254649/1‑93 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          <  1:1384985/93‑1 (MQ=255)
 gatgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttc                                                                                          <  2:3306351/93‑1 (MQ=255)
  atgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcc                                                                                         >  1:2802642/1‑93 (MQ=255)
  atgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttcc                                                                                         <  2:3797473/93‑1 (MQ=255)
  atgatCCTGCCAATTGCCTTGTTAACCACTGGCTTGTTAGGAACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCGACCGATTttcc                                                                                         >  1:1953257/1‑93 (MQ=255)
   tgatCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTttccc                                                                                        <  2:693683/93‑1 (MQ=255)
      tCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTtttccagt                                                                                     >  2:3800272/1‑90 (MQ=255)
      gccTGCTTATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAATGCCTGATGATTttccccgt                                                                                     >  2:4225861/2‑90 (MQ=255)
       ccTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGAGGTTTtcccccttc                                                                                   >  1:1435873/1‑93 (MQ=255)
         tGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtccccgttcc                                                                                  >  1:1076674/1‑93 (MQ=255)
          gCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTtccccgttccc                                                                                 >  1:2491702/1‑93 (MQ=255)
            cAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCgg                                                                               <  2:1047572/93‑1 (MQ=255)
             aaTTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGt                                                                              >  2:1193311/1‑93 (MQ=255)
               ttGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTg                                                                            >  1:686132/1‑93 (MQ=255)
                 tccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                                                          <  2:2032675/92‑1 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCTCGGTTGCt                                                                          >  1:102280/1‑93 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                                                          >  2:3940688/1‑93 (MQ=255)
                 gCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCt                                                                          >  2:2874475/1‑93 (MQ=255)
                  ccTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTg                                                                         >  2:426973/1‑93 (MQ=255)
                    ttGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTa                                                                       >  2:786026/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                      >  1:3626024/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                      <  1:3693162/93‑1 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTAc                                                                      >  1:4463955/1‑93 (MQ=255)
                     tGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGCTGATTAAAAAAAAATACTGTAACGCCTGTCGGTTTTCCCGCTTCGCGGTTGCTGTAc                                                                      >  2:583565/1‑93 (MQ=255)
                       ttAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCg                                                                    <  2:2613792/93‑1 (MQ=255)
                        tAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCgg                                                                   <  2:392307/93‑1 (MQ=255)
                         aaCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCggg                                                                  <  2:4746260/93‑1 (MQ=255)
                          aCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGa                                                                 <  2:4487110/93‑1 (MQ=255)
                          aCCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGCTCCCGGTTGCTGTACCGGGa                                                                 >  1:3528095/1‑93 (MQ=255)
                           ccGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGaa                                                                <  2:2974649/93‑1 (MQ=255)
                           ccGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGaa                                                                >  1:4617930/1‑93 (MQ=255)
                           ccGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGaa                                                                <  2:3556445/93‑1 (MQ=255)
                            cGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAAc                                                               >  2:3094577/1‑93 (MQ=255)
                              cTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGt                                                             >  1:509108/1‑93 (MQ=255)
                              cTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTGTACCGGGGACGt                                                             >  1:4239882/1‑93 (MQ=255)
                               tGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTa                                                            >  2:2940130/1‑93 (MQ=255)
                               tGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTa                                                            >  1:822256/1‑93 (MQ=255)
                                ggCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTAt                                                           <  1:98300/93‑1 (MQ=255)
                                 gCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGGTGCTGTACCGGGAACGTAtt                                                          >  2:2047232/1‑93 (MQ=255)
                                   ttGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTa                                                        >  1:3259571/1‑93 (MQ=255)
                                   ttGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTa                                                        >  1:625479/1‑93 (MQ=255)
                                    tGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTaa                                                       <  1:3087191/93‑1 (MQ=255)
                                     gTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAAt                                                      <  1:2008092/93‑1 (MQ=255)
                                      gtAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAAtt                                                     <  1:1535755/92‑1 (MQ=255)
                                          gCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTcccc                                                 <  1:864688/93‑1 (MQ=255)
                                            aCGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTg                                               >  1:2315194/1‑93 (MQ=255)
                                             cGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGc                                              >  1:3044442/1‑93 (MQ=255)
                                             cGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGc                                              <  2:4094867/93‑1 (MQ=255)
                                             cGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGc                                              <  1:2495554/93‑1 (MQ=255)
                                             cGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTACCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGc                                              >  1:1764828/1‑93 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                             >  2:2139401/1‑93 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                             <  1:3352796/93‑1 (MQ=255)
                                              gCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCa                                             <  2:1505496/93‑1 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                           <  2:3126991/93‑1 (MQ=255)
                                                tGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATc                                           >  2:740333/1‑93 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                         >  1:1550644/1‑93 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                         >  1:2654296/1‑93 (MQ=255)
                                                  tttGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGc                                         >  2:162551/1‑93 (MQ=255)
                                                     gtggtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCg                                      >  1:3565481/1‑93 (MQ=255)
                                                        gtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAt                                   <  1:396752/93‑1 (MQ=255)
                                                        gtgATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCAt                                   <  1:2946685/93‑1 (MQ=255)
                                                          gATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTc                                 >  2:2794483/1‑93 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                <  2:1953830/93‑1 (MQ=255)
                                                           aTTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCt                                <  1:4373270/93‑1 (MQ=255)
                                                              aaaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGc                             >  2:1265541/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            >  1:846252/1‑93 (MQ=255)
                                                               aaaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGcc                            <  2:689787/93‑1 (MQ=255)
                                                                aaaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCg                           <  2:2968364/93‑1 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTGCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                          <  1:4299963/93‑1 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                          >  2:4581580/1‑93 (MQ=255)
                                                                 aaaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGc                          >  2:4299420/1‑93 (MQ=255)
                                                                  aaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGct                         >  2:4410410/1‑92 (MQ=255)
                                                                  aaaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCa                         <  1:2872239/93‑1 (MQ=255)
                                                                   aaaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCAt                        >  1:2041718/1‑93 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                       >  1:2603620/1‑93 (MQ=255)
                                                                    aaaTACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATc                       >  1:4378668/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTTCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                    >  2:1114718/1‑93 (MQ=255)
                                                                       tACTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTc                    >  2:2327253/1‑93 (MQ=255)
                                                                        aCTGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcc                   >  2:1341864/1‑93 (MQ=255)
                                                                          tGTTACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                 >  2:4644832/1‑93 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                 >  1:1296643/1‑93 (MQ=255)
                                                                          tGTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTcccc                 <  1:3034327/93‑1 (MQ=255)
                                                                           gTAACGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccc                >  1:4737934/1‑93 (MQ=255)
                                                                             aaCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTccccccg              >  2:83091/1‑91 (MQ=255)
                                                                             aaCGCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGGACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCgg              >  2:341554/1‑93 (MQ=255)
                                                                                gCCTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGt           >  1:3293315/1‑93 (MQ=255)
                                                                                  cTGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTca         >  1:4560796/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  1:4435582/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  2:4616382/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  2:609419/1‑93 (MQ=255)
                                                                                   tGACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTcac        >  1:752612/1‑93 (MQ=255)
                                                                                    gACGATTTTCCCCGTTCCCGGTTGCTGTACCGGGCACGTATTTAATTCACCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCCGCGTcaca       >  1:2301446/1‑93 (MQ=255)
                                                                                     aCGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCacac      >  1:1292172/1‑93 (MQ=255)
                                                                                      gattttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCg    <  2:3967808/90‑1 (MQ=255)
                                                                                      cGATTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACAcc     <  1:1181428/93‑1 (MQ=255)
                                                                                       attttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa   <  2:1832435/91‑1 (MQ=255)
                                                                                        ttttcCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGaa  <  2:3690886/92‑1 (MQ=255)
                                                                                        aTTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa   >  1:3071449/1‑93 (MQ=255)
                                                                                        aTTTTCCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGa   >  1:4347594/1‑93 (MQ=255)
                                                                                              |                                                                                         
TGATGATCCTGCCAATTGCCTTGTTAACCGCTGGCTTGTTAGGCACGCTGTTTGTGGTGATTAAAAAAAAATACTGTAACGCCTGACGATTTTCCCGTTCCCGGTTGCTGTACCGGGAACGTATTTAATTCCCCTGCATCGCCCGCATTCTTGCCGCATCTTCCCCCGGCGTCACACCGAA  >  NC_000913/3152152‑3152332

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