Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,247,901 T→C D147G (GAC→GGC)  ycgC ← putative PTS system enzyme I/IIA/HPr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,247,9010TC99.5% 648.0 / ‑4.9 200D147G (GAC→GGC) ycgCputative PTS system enzyme I/IIA/HPr
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (80/118);  total (82/118)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 205 total aligned reads displayed.

GCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATA  >  NC_000913/1247812‑1247990
                                                                                         |                                                                                         
gCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCag                                                                                        >  2:1332169/1‑93 (MQ=255)
gCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCag                                                                                        >  1:118686/1‑93 (MQ=255)
gCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCag                                                                                        >  2:5440124/1‑93 (MQ=255)
gCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCag                                                                                        >  2:4242898/1‑93 (MQ=255)
   gggAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAgaga                                                                                     >  1:1878274/1‑93 (MQ=255)
          cGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgt                                                                              <  2:2476442/93‑1 (MQ=255)
           ggACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgtg                                                                             >  2:627480/1‑93 (MQ=255)
           ggACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgtg                                                                             >  2:5812874/1‑93 (MQ=255)
            gACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgtgt                                                                            <  1:1073977/93‑1 (MQ=255)
            gACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgtgt                                                                            >  1:2216241/1‑93 (MQ=255)
            gACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAgtgt                                                                            <  2:2579770/93‑1 (MQ=255)
               gTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgg                                                                         <  2:5480314/93‑1 (MQ=255)
               ctACATGCAGGCCTTTACGGTTGTTTATGACCACCGCCAGAGAACGGGCTTCTTAATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgg                                                                         <  1:772036/92‑1 (MQ=255)
                tACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgga                                                                        >  2:5727836/1‑93 (MQ=255)
                tACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgga                                                                        >  1:2255393/1‑93 (MQ=255)
                tACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgga                                                                        >  1:5893860/1‑93 (MQ=255)
                tACATGCAGGCCGTTACGGCTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCgga                                                                        >  2:1224747/1‑93 (MQ=255)
                 tcATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggag                                                                       <  2:714017/92‑1 (MQ=255)
                 tcATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggag                                                                       <  2:749473/92‑1 (MQ=255)
                  cATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggagg                                                                      >  2:4129722/1‑93 (MQ=255)
                  cATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggagg                                                                      <  1:6687107/93‑1 (MQ=255)
                   aTGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggagga                                                                     <  2:6465983/93‑1 (MQ=255)
                   aTGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCggagga                                                                     <  1:6028646/93‑1 (MQ=255)
                    tGCAGGCCGTTACGGTCTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGATGAc                                                                    >  1:3349069/1‑93 (MQ=255)
                     gCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACg                                                                   >  2:6248505/1‑93 (MQ=255)
                      cAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACgg                                                                  <  2:106531/93‑1 (MQ=255)
                       aGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTTGGAGGACGGt                                                                 <  1:4072450/93‑1 (MQ=255)
                       aGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGt                                                                 >  1:864147/1‑93 (MQ=255)
                        ggCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTa                                                                >  1:2487678/1‑93 (MQ=255)
                        ggCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTa                                                                >  1:3726807/1‑93 (MQ=255)
                         gCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTaa                                                               <  1:1024379/93‑1 (MQ=255)
                         gCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTaa                                                               >  2:4343923/1‑93 (MQ=255)
                         gCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTaa                                                               >  1:3591606/1‑93 (MQ=255)
                            gTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAAcc                                                            <  1:318703/93‑1 (MQ=255)
                             ttACGGTTTTTTATGACCACCGCCAGAGAACTGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAAccc                                                           <  1:6021808/93‑1 (MQ=255)
                             ttACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAAccc                                                           >  1:4265630/1‑93 (MQ=255)
                              tACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCa                                                          <  1:3421743/93‑1 (MQ=255)
                              tACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCa                                                          <  2:2592773/93‑1 (MQ=255)
                                cGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgt                                                        <  1:5147269/93‑1 (MQ=255)
                                cGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgt                                                        >  1:1773095/1‑93 (MQ=255)
                                   ttttttATGACCCCCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgttgt                                                     <  1:3214484/93‑1 (MQ=255)
                                   ttttttATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgttgt                                                     <  2:2907901/93‑1 (MQ=255)
                                    tttttATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgttgtt                                                    <  2:2872538/93‑1 (MQ=255)
                                      tttATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCa                                                  <  1:6159963/93‑1 (MQ=255)
                                           gACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGttt                                             <  1:6712881/93‑1 (MQ=255)
                                             ccaccGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTgg                                           <  1:6242457/93‑1 (MQ=255)
                                               gccGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGGTGTTCACGTTTGGCt                                         >  1:1751196/2‑93 (MQ=255)
                                               cccGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCt                                         <  1:6332188/92‑1 (MQ=255)
                                               accGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCt                                         <  1:2769252/93‑1 (MQ=255)
                                               accGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCt                                         >  2:4443924/1‑93 (MQ=255)
                                               accGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCt                                         >  2:4894361/1‑93 (MQ=255)
                                               accGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCt                                         <  1:1733498/93‑1 (MQ=255)
                                                ccGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGATGTTCACGTTTGGCtt                                        >  1:2831817/1‑93 (MQ=255)
                                                 cGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGCAAACCCAGTTGTTCACGTTTGGCTTc                                       >  2:1983280/1‑93 (MQ=255)
                                                     agagAACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAgc                                   >  1:1474778/1‑93 (MQ=255)
                                                      gagaACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAgcg                                  <  2:3034410/93‑1 (MQ=255)
                                                      gagaACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAgcg                                  >  1:6578383/1‑93 (MQ=255)
                                                      gagaACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGATCACGTTTGGCTTCCAgcg                                  >  2:1714442/1‑93 (MQ=255)
                                                       agaACGGGCTTCTGCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAgcgc                                 >  1:6435938/1‑93 (MQ=255)
                                                        gaACGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgca                                <  2:1763930/93‑1 (MQ=255)
                                                         aaCGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcat                               <  1:5454614/93‑1 (MQ=255)
                                                         aaCGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcat                               <  2:3637825/93‑1 (MQ=255)
                                                          aCGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcatg                              >  2:3442786/1‑93 (MQ=255)
                                                          aCGGGCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcatg                              >  2:5360277/1‑93 (MQ=255)
                                                            gggCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcatgca                            <  1:1557100/93‑1 (MQ=255)
                                                            gggCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcatgca                            <  1:1407994/93‑1 (MQ=255)
                                                            gggCTTCTTCATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCgcatgca                            <  2:63861/93‑1 (MQ=255)
                                                               gttcttcATCGTACGCAGGACATGTGCAAGAGATGTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATgg                         <  2:4906389/92‑1 (MQ=255)
                                                                ttcttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGc                        <  1:4655782/93‑1 (MQ=255)
                                                                ttcttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGc                        <  1:5541493/93‑1 (MQ=255)
                                                                 tcttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCg                       >  1:3935517/1‑93 (MQ=255)
                                                                   ttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTc                     >  1:3821409/1‑93 (MQ=255)
                                                                   ttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTc                     <  2:3607992/93‑1 (MQ=255)
                                                                   ttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTc                     <  1:2470920/93‑1 (MQ=255)
                                                                   ttcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTc                     <  1:973814/93‑1 (MQ=255)
                                                                    tcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                    <  2:5241869/93‑1 (MQ=255)
                                                                    tcATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                    <  2:1526478/93‑1 (MQ=255)
                                                                     cATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCaa                   >  1:4945442/1‑93 (MQ=255)
                                                                     cATCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCaa                   <  2:4605626/93‑1 (MQ=255)
                                                                      aTCGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCaaa                  >  1:2719674/1‑93 (MQ=255)
                                                                        cGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGa                <  2:3130879/93‑1 (MQ=255)
                                                                        cGTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGa                <  1:3612732/93‑1 (MQ=255)
                                                                         gTACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGAt               >  2:2040624/1‑93 (MQ=255)
                                                                          tACGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATa              >  1:5688071/1‑93 (MQ=255)
                                                                           aCGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa             <  2:6087499/93‑1 (MQ=255)
                                                                            cGCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAAc            <  1:526684/93‑1 (MQ=255)
                                                                             gCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt           >  1:5075408/1‑93 (MQ=255)
                                                                             gCAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt           >  2:4836864/1‑93 (MQ=255)
                                                                              cAGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACtt          <  2:3172910/93‑1 (MQ=255)
                                                                               aGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGTATGGCGTCAAAGATAACttt         <  1:164496/93‑1 (MQ=255)
                                                                               aGGACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACttt         <  1:2157389/93‑1 (MQ=255)
                                                                                ggACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg        >  1:3055442/1‑93 (MQ=255)
                                                                                ggACATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg        >  1:4364311/1‑93 (MQ=255)
                                                                                  aCATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTc      <  1:2987755/93‑1 (MQ=255)
                                                                                   cATGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCg     >  2:5058535/1‑93 (MQ=255)
                                                                                    aTGTGCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGa    <  1:4876164/93‑1 (MQ=255)
                                                                                     tgtgCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGat   <  2:4078989/93‑1 (MQ=255)
                                                                                     tgtgCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGat   >  2:1141397/1‑93 (MQ=255)
                                                                                     tgtgCCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGat   >  1:5642169/1‑93 (MQ=255)
                                                                                      gtgcCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGata  <  1:21007/93‑1 (MQ=255)
                                                                                         |                                                                                         
GCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATA  >  NC_000913/1247812‑1247990

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