Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 52,760 A→G G219G (GGT→GGC pdxA ← PdxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091352,7600AG99.6% 774.1 / ‑5.2 241G219G (GGT→GGCpdxAPdxA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (133/107);  total (134/107)
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.75e-01

Only 100 of 250 total aligned reads displayed.

CAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTG  >  NC_000913/52670‑52849
                                                                                          |                                                                                         
 aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGccc                                                                                        <  1:3463504/93‑1 (MQ=255)
 aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTGCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGccc                                                                                        <  2:2922233/93‑1 (MQ=255)
  gggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCa                                                                                       >  1:1765534/1‑93 (MQ=255)
   ggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCat                                                                                      <  1:5072580/93‑1 (MQ=255)
    gTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCata                                                                                     <  1:5486769/93‑1 (MQ=255)
    gTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCata                                                                                     >  2:6694072/1‑93 (MQ=255)
    gTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCata                                                                                     <  1:5515761/93‑1 (MQ=255)
     tATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCatat                                                                                    >  2:1195685/1‑93 (MQ=255)
     tATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCatat                                                                                    <  2:3616894/93‑1 (MQ=255)
      aTCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATg                                                                                   <  2:5808899/93‑1 (MQ=255)
        cGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGAc                                                                                 <  2:4583570/93‑1 (MQ=255)
        cGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGAc                                                                                 <  1:4439845/93‑1 (MQ=255)
        cGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGAc                                                                                 <  1:3708537/93‑1 (MQ=255)
         ggcaggcaGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGAcc                                                                                <  1:5620872/93‑1 (MQ=255)
           caggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGc                                                                              >  2:261887/1‑93 (MQ=255)
             ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                            >  1:245942/1‑93 (MQ=255)
             ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                            >  2:4003055/1‑93 (MQ=255)
             ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                            >  2:1368175/1‑93 (MQ=255)
             ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                            <  2:4168218/93‑1 (MQ=255)
             ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCt                                                                            <  1:1515391/93‑1 (MQ=255)
              gcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCtt                                                                           >  2:3547066/1‑93 (MQ=255)
                 gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGc                                                                        >  1:1994220/1‑93 (MQ=255)
                 gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGc                                                                        >  2:2411872/1‑93 (MQ=255)
                 gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGc                                                                        >  1:3347020/1‑93 (MQ=255)
                  cGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGcc                                                                       <  1:2562138/93‑1 (MQ=255)
                   ggCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGccg                                                                      <  1:6437969/93‑2 (MQ=255)
                   ggCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGccc                                                                      <  1:3797445/93‑1 (MQ=255)
                    gCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCcg                                                                     <  1:3480368/93‑1 (MQ=255)
                    gCCCGTTCATTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCcg                                                                     <  2:5189354/93‑1 (MQ=255)
                       cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGt                                                                  >  1:6113814/1‑93 (MQ=255)
                        gTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTg                                                                 >  1:6290931/1‑93 (MQ=255)
                        gTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTg                                                                 <  2:5692139/93‑1 (MQ=255)
                           gAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCgg                                                              <  2:3821220/93‑1 (MQ=255)
                            aGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGa                                                             <  2:3578441/93‑1 (MQ=255)
                             gTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAt                                                            <  1:790020/93‑1 (MQ=255)
                             gTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAt                                                            >  2:6537071/1‑93 (MQ=255)
                              tttCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAtt                                                           <  1:1059513/93‑1 (MQ=255)
                              tttCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAtt                                                           >  1:3560835/1‑93 (MQ=255)
                              tttCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGAtt                                                           <  1:3151295/93‑1 (MQ=255)
                               ttCAGCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTc                                                          <  2:5358274/93‑1 (MQ=255)
                                   tCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGcc                                                      >  1:5004335/1‑93 (MQ=255)
                                    ccccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGccc                                                     <  2:1393392/93‑1 (MQ=255)
                                    ccccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGccc                                                     >  2:5010683/1‑93 (MQ=255)
                                     cccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCg                                                    <  2:2094018/93‑1 (MQ=255)
                                     cccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCg                                                    >  1:2278044/1‑93 (MQ=255)
                                      ccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGc                                                   <  2:5563097/93‑1 (MQ=255)
                                       cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCa                                                  >  2:303769/1‑93 (MQ=255)
                                       cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCa                                                  >  1:1662564/1‑93 (MQ=255)
                                         gcgcCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGa                                                <  1:1448086/93‑1 (MQ=255)
                                            cccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCa                                             <  1:5383378/93‑1 (MQ=255)
                                            cccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCa                                             >  1:5697468/1‑93 (MQ=255)
                                             ccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAg                                            >  2:5423379/1‑93 (MQ=255)
                                             ccGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAg                                            >  2:4218780/1‑93 (MQ=255)
                                              cGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGa                                           >  2:6025890/1‑93 (MQ=255)
                                              cGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGa                                           >  1:3590527/1‑93 (MQ=255)
                                                 aGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATg                                        >  2:4071064/1‑93 (MQ=255)
                                                 aGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATg                                        <  1:6039328/93‑1 (MQ=255)
                                                  gCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgc                                       >  1:1699323/1‑93 (MQ=255)
                                                  gCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgc                                       >  1:194455/1‑93 (MQ=255)
                                                   cTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcg                                      <  1:2162267/93‑1 (MQ=255)
                                                   cTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcg                                      <  1:2970635/93‑1 (MQ=255)
                                                   cTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcg                                      <  1:4681095/93‑1 (MQ=255)
                                                    tCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcgc                                     >  2:5106461/1‑93 (MQ=255)
                                                    tCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgcgc                                     >  1:2718637/1‑93 (MQ=255)
                                                     cATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGACTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGcgcg                                    >  2:2649689/1‑93 (MQ=255)
                                                      aTTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCgg                                   <  1:2391443/93‑1 (MQ=255)
                                                       ttGAGCACCAGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGt                                  <  1:4618790/93‑1 (MQ=255)
                                                          aGCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCg                               >  2:6264232/1‑93 (MQ=255)
                                                           gCACCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCgg                              >  1:4225166/1‑93 (MQ=255)
                                                            cACCGGAATAATGTTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGc                             <  1:4493628/93‑1 (MQ=255)
                                                             aCCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCa                            <  2:3692372/93‑1 (MQ=255)
                                                             aCCGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCa                            >  2:2213256/1‑93 (MQ=255)
                                                              ccGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCaa                           >  1:5924616/1‑93 (MQ=255)
                                                              ccGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCaa                           <  2:1813238/93‑1 (MQ=255)
                                                               cGGAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAAt                          >  1:1177643/1‑93 (MQ=255)
                                                                 gAATAATGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATAc                        <  2:5304332/93‑1 (MQ=255)
                                                                   ataatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCa                      <  2:3324390/93‑1 (MQ=255)
                                                                     aatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaaa                    <  1:724837/93‑1 (MQ=255)
                                                                     aatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaaa                    >  2:1930849/1‑93 (MQ=255)
                                                                     aatGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAACAGAATGCGCGGTTCGGCAATACCaaa                    >  2:3900515/1‑93 (MQ=255)
                                                                      atGGTGTCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAAt                   >  1:1011871/1‑93 (MQ=255)
                                                                          tgtCTATCTCTTCCGTGCCCATATGCCCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGCCCAGAATGCGCGGTTCGGCAATACCAAATTTgc               <  1:4021043/93‑2 (MQ=255)
                                                                          tgtCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg               >  2:498986/1‑93 (MQ=255)
                                                                          tgtCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg               >  2:6351834/1‑93 (MQ=255)
                                                                           gtCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGt              >  2:1481095/1‑93 (MQ=255)
                                                                            tCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTc             >  1:1243449/1‑93 (MQ=255)
                                                                            tCTATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTc             >  1:2871436/1‑93 (MQ=255)
                                                                              tATCTCTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCg           >  1:2811469/1‑93 (MQ=255)
                                                                                tctcTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCa         >  2:3084236/1‑93 (MQ=255)
                                                                                tctcTTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCa         >  2:1555243/1‑93 (MQ=255)
                                                                                  tctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCaaa       <  1:6524781/93‑1 (MQ=255)
                                                                                   ctTCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAAt      >  2:2972047/1‑93 (MQ=255)
                                                                                    ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATc     <  2:731530/93‑1 (MQ=255)
                                                                                    ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATc     >  1:1378324/1‑93 (MQ=255)
                                                                                    ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATc     <  1:277838/93‑1 (MQ=255)
                                                                                    ttCCGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATc     <  2:6188549/93‑1 (MQ=255)
                                                                                      ccGTGCCCATATGACCGCTTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt   <  2:4412651/93‑1 (MQ=255)
                                                                                      ccGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt   >  2:6485999/1‑93 (MQ=255)
                                                                                      ccGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt   >  1:4866112/1‑93 (MQ=255)
                                                                                       cGTGCCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTg  <  2:394890/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTG  >  NC_000913/52670‑52849

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