Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,975,807 A→G L372P (CTG→CCG)  lysA ← diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,975,8070AG99.4% 1059.0 / ‑3.7 336L372P (CTG→CCG) lysAdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (166/168);  total (168/168)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

Only 100 of 339 total aligned reads displayed.

GCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975717‑2975896
                                                                                          |                                                                                         
 ccTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTa                                                                                        >  1:7310568/1‑93 (MQ=255)
   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      >  1:2880934/1‑93 (MQ=255)
   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      >  2:4415774/1‑93 (MQ=255)
   tGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTAcc                                                                                      <  1:6829892/93‑1 (MQ=255)
    gACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCa                                                                                     <  1:4208749/93‑1 (MQ=255)
    gACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCa                                                                                     >  2:4601754/1‑93 (MQ=255)
      ccATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGa                                                                                   >  1:4024805/1‑93 (MQ=255)
       cATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGAt                                                                                  >  2:5236049/1‑93 (MQ=255)
       cATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGAt                                                                                  <  2:1286222/93‑1 (MQ=255)
       cATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGAt                                                                                  <  2:6217225/93‑1 (MQ=255)
       cATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGAt                                                                                  <  1:4354296/93‑1 (MQ=255)
        aTTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATa                                                                                 >  2:2081393/1‑93 (MQ=255)
         ttATGAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                <  1:6654041/93‑1 (MQ=255)
         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCTCTGTATCATGCGGTACCAGATaa                                                                                >  2:4307340/1‑93 (MQ=255)
         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                >  2:7290654/1‑93 (MQ=255)
         ttATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATaa                                                                                >  2:5551300/1‑93 (MQ=255)
          tATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAAt                                                                               <  1:6945683/93‑1 (MQ=255)
           aTCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATc                                                                              >  2:4195304/1‑93 (MQ=255)
           aTCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATc                                                                              <  1:6344153/93‑1 (MQ=255)
            tCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCa                                                                             >  1:179694/1‑93 (MQ=255)
            tCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCa                                                                             <  1:5000795/93‑1 (MQ=255)
             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            <  2:1684527/93‑1 (MQ=255)
             cAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAc                                                                            >  1:1445592/1‑93 (MQ=255)
              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAct                                                                           >  1:2254199/1‑92 (MQ=255)
              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           >  2:117905/1‑93 (MQ=255)
              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           <  1:2768860/93‑1 (MQ=255)
              aaaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCAcc                                                                           <  2:5992270/93‑1 (MQ=255)
               aaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCt                                                                          >  1:6454519/1‑93 (MQ=255)
               aaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCt                                                                          <  1:4133115/93‑1 (MQ=255)
               aaCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCt                                                                          <  2:3289962/93‑1 (MQ=255)
                aCAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTg                                                                         >  2:2236514/1‑93 (MQ=255)
                   gtaCTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCTCCATATGCCCCTGTATCATGCGGTACCAGAGAATCACCTGCCt                                                                      >  2:7034388/3‑93 (MQ=255)
                   gAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCt                                                                      <  1:624393/93‑1 (MQ=255)
                    aaCTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCtt                                                                     >  2:6514617/1‑93 (MQ=255)
                       ttCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCAc                                                                  <  2:3474811/93‑1 (MQ=255)
                         cTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACtt                                                                >  1:829785/1‑93 (MQ=255)
                          tGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTc                                                               >  1:3684908/1‑93 (MQ=255)
                           ggTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcc                                                              <  1:721676/93‑1 (MQ=255)
                            gTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTcgg                                                             >  1:2541944/1‑91 (MQ=255)
                             tAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCgg                                                            >  2:3113606/1‑93 (MQ=255)
                              aaCAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGc                                                           >  1:5204950/1‑93 (MQ=255)
                                cAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCaa                                                         >  1:527425/1‑93 (MQ=255)
                                 aGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAg                                                        >  2:4950877/1‑93 (MQ=255)
                                  gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       >  2:7886138/1‑93 (MQ=255)
                                  gCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAgg                                                       <  2:406931/93‑1 (MQ=255)
                                   cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      <  2:4856117/93‑1 (MQ=255)
                                   cGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgc                                                      <  1:3548703/93‑1 (MQ=255)
                                     gACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcgc                                                    <  1:5058728/93‑1 (MQ=255)
                                     gACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcgc                                                    >  2:316959/1‑93 (MQ=255)
                                     gACGGCTATTGTAGTTGGATGAAATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGgcga                                                    >  2:2197941/1‑92 (MQ=255)
                                      aCGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGcgcg                                                   <  2:7674728/93‑1 (MQ=255)
                                         gCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGt                                                >  2:4930914/1‑93 (MQ=255)
                                         gCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGt                                                >  1:1565410/1‑93 (MQ=255)
                                           tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  1:4218563/1‑93 (MQ=255)
                                           tATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGttt                                              >  1:3327504/1‑93 (MQ=255)
                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             <  1:6417410/93‑1 (MQ=255)
                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             <  2:2824726/93‑1 (MQ=255)
                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             <  2:1429939/93‑1 (MQ=255)
                                            aTTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTc                                             >  1:7639631/1‑93 (MQ=255)
                                             ttGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCa                                            >  1:173000/1‑93 (MQ=255)
                                              tGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCaa                                           >  1:866458/1‑93 (MQ=255)
                                              tGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGATTCaa                                           >  2:7618485/1‑93 (MQ=255)
                                               gTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAAc                                          <  2:4636641/93‑1 (MQ=255)
                                               gTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAAc                                          >  2:2533308/1‑93 (MQ=255)
                                                    tGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTc                                     <  1:1208746/93‑1 (MQ=255)
                                                    tGGATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTc                                     <  2:1366629/93‑1 (MQ=255)
                                                     ggATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTcc                                    >  2:6425328/1‑93 (MQ=255)
                                                     ggATGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGGTTCAACATTTcc                                    >  1:5231789/1‑93 (MQ=255)
                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  2:5818771/93‑1 (MQ=255)
                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 <  1:37411/93‑1 (MQ=255)
                                                        tGACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTccccc                                 >  1:3993447/1‑93 (MQ=255)
                                                         gACATTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCt                                >  1:2194518/1‑93 (MQ=255)
                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             <  1:7811526/93‑1 (MQ=255)
                                                            aTTGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTcc                             >  2:7136812/1‑93 (MQ=255)
                                                             ttGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCct                            <  1:5503652/93‑1 (MQ=255)
                                                             ttGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCct                            >  1:490835/1‑93 (MQ=255)
                                                             ttGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCct                            <  1:1038556/93‑1 (MQ=255)
                                                             ttGATGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCct                            >  1:944546/1‑93 (MQ=255)
                                                                aTGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgct                         <  2:672849/93‑1 (MQ=255)
                                                                 tGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgctg                        <  1:5918248/93‑1 (MQ=255)
                                                                 tGCGCCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCctgctg                        >  1:4414874/1‑93 (MQ=255)
                                                                  gcgcCATATGCCCCTGTATCATGCGGTACCAGATAATCCCCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTgg                       <  1:2640326/93‑1 (MQ=255)
                                                                  gcgcCATATGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTgg                       <  2:2980493/93‑1 (MQ=255)
                                                                       atatGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTaaa                  <  1:2687792/93‑1 (MQ=255)
                                                                         atGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGa                >  1:2417068/1‑93 (MQ=255)
                                                                         atGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGa                <  2:831876/93‑1 (MQ=255)
                                                                         atGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGa                <  2:6140301/93‑1 (MQ=255)
                                                                         atGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGa                <  1:7708814/93‑1 (MQ=255)
                                                                          tGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCTCCCTTCCTGCTGGGTCAAGAc               >  2:1740788/1‑93 (MQ=255)
                                                                          tGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGAc               >  2:5112591/1‑93 (MQ=255)
                                                                          tGCCCCTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGAc               <  1:1935772/93‑1 (MQ=255)
                                                                             cccTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATc            <  1:4705807/93‑1 (MQ=255)
                                                                               cTGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGc          >  2:7566151/1‑93 (MQ=255)
                                                                                tGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGcc         >  2:3717163/1‑93 (MQ=255)
                                                                                tGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGcc         <  1:1282026/93‑1 (MQ=255)
                                                                                tGTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGcc         >  2:3576874/1‑93 (MQ=255)
                                                                                 gTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGccc        <  1:748405/93‑1 (MQ=255)
                                                                                 gTATCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGccc        <  1:2324405/93‑1 (MQ=255)
                                                                                  tctCATTCCGTCCCAGATAATCCCCTGACTTCGCTTCCGGCAAGGCGCGGGTTTCAACTATTCCCCCTTCCTGCTGGGTAAAGACATCGCCCg       <  2:5793894/91‑1 (MQ=255)
                                                                                  aaTCATGCGGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCg       >  2:6943971/2‑93 (MQ=255)
                                                                                          |                                                                                         
GCCTGACCATTATCAAACAGAACTTCTGGTAACAGCGGACGGCTATTGTAGTTGGATGACATTGATGCGCCATATGCCCCTGTATCATGCAGTACCAGATAATCACCTGCCTTCACTTCCGGCAAGGCGCGGGTTTCAACATTTCCCCCTTCCTGCTGGGTAAAGACATCGCCCGATTCA  >  NC_000913/2975717‑2975896

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