Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,237,242 +C intergenic (+55/‑342) ygjT → / → ygjU putative transport protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,237,2401.C98.7% 475.3 / ‑1.1 156intergenic (+53/‑344)ygjT/ygjUputative transport protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (86/68);  total (86/70)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 185 total aligned reads displayed.

CTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCA  >  NC_000913/3237148‑3237325
                                                                                             |                                                                                      
cTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtct                                                                                         <  1:177973/93‑1 (MQ=255)
 tGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctc                                                                                       <  1:2723858/93‑1 (MQ=255)
  gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctcc                                                                                      >  2:4288930/1‑93 (MQ=255)
  gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctcc                                                                                      >  1:3877878/1‑93 (MQ=255)
  gggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctcc                                                                                      >  1:205089/1‑93 (MQ=255)
  gggTGAATTATCGGCATGATAAGCAGCGGGGAGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctcc                                                                                      >  2:4190614/1‑93 (MQ=255)
   ggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccc                                                                                     <  2:2466666/93‑2 (MQ=255)
   ggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccc                                                                                     <  1:2098074/93‑2 (MQ=255)
    gTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccct                                                                                    <  2:1431394/93‑3 (MQ=255)
    gTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTtctccct                                                                                    <  1:223564/93‑3 (MQ=255)
        aTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                                <  1:72246/93‑1 (MQ=255)
         ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                               >  1:539296/1‑92 (MQ=255)
           aTCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccctt                                                                             >  1:2074523/1‑90 (MQ=255)
            tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcccttc                                                                            >  1:1045885/1‑89 (MQ=255)
            tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcccttc                                                                            >  1:3615634/1‑89 (MQ=255)
              ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                           <  2:3214581/93‑1 (MQ=255)
              ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                           >  1:767988/1‑93 (MQ=255)
               gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                          <  1:668071/93‑1 (MQ=255)
               gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                          <  2:2745129/93‑1 (MQ=255)
               gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                          <  2:4411050/93‑1 (MQ=255)
               gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                          <  1:1915982/93‑1 (MQ=255)
                 aTGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                        >  2:3732009/1‑93 (MQ=255)
                  tGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                       <  1:4614689/93‑1 (MQ=255)
                  tGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                       >  2:3583448/1‑93 (MQ=255)
                  tGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                       <  1:2134429/93‑1 (MQ=255)
                  tGATAAGCAGCGGGGGGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                       >  2:2631416/1‑93 (MQ=255)
                   gATAAGCAGCGGGGTGTATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                      <  2:2555723/93‑1 (MQ=255)
                   gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                      <  1:117948/93‑1 (MQ=255)
                   gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                      <  2:3275572/93‑1 (MQ=255)
                    aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgagg                                                                     <  1:3354095/93‑2 (MQ=255)
                    aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgagg                                                                     <  1:4406982/93‑2 (MQ=255)
                    aTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgagg                                                                     <  2:3195030/93‑2 (MQ=255)
                     tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTCCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaggg                                                                    <  2:1790270/93‑3 (MQ=255)
                      aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAgggg                                                                   <  2:1337911/93‑1 (MQ=255)
                       agcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGgggg                                                                  <  2:3671634/93‑1 (MQ=255)
                        gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGt                                                                 <  2:2805972/93‑1 (MQ=255)
                           gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCgg                                                              >  1:3428434/1‑93 (MQ=255)
                           gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCgg                                                              <  1:2615888/93‑1 (MQ=255)
                            cGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCggg                                                             >  1:1996204/1‑93 (MQ=255)
                              gggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGt                                                           <  1:762792/93‑1 (MQ=255)
                               ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTg                                                          >  2:2454508/1‑93 (MQ=255)
                                gTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGa                                                         >  1:3013331/1‑93 (MQ=255)
                                 tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAg                                                        <  1:864956/93‑1 (MQ=255)
                                 tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAg                                                        <  2:971537/93‑1 (MQ=255)
                                   gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                      <  1:1073197/93‑1 (MQ=255)
                                   gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                      >  1:3709509/1‑93 (MQ=255)
                                     tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAgggga                                                    >  2:2827224/1‑91 (MQ=255)
                                     tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTa                                                    >  1:4094166/1‑93 (MQ=38)
                                     tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTa                                                    >  1:3600691/1‑93 (MQ=38)
                                     tAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGAGAGGGTa                                                    >  1:2133968/1‑93 (MQ=255)
                                        tttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAtt                                                 >  1:3337201/1‑93 (MQ=38)
                                          tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAtttt                                               >  2:568835/1‑93 (MQ=38)
                                          tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAtttt                                               <  1:4191695/93‑1 (MQ=255)
                                           ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAttttt                                              >  1:252103/1‑93 (MQ=38)
                                             aaTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCg                                            >  1:3286135/1‑93 (MQ=38)
                                              aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGc                                           >  2:477637/1‑93 (MQ=38)
                                              aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGc                                           >  1:1484855/1‑93 (MQ=38)
                                               tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCa                                          >  2:1599250/1‑93 (MQ=38)
                                               tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCa                                          <  2:342694/93‑1 (MQ=37)
                                               tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCa                                          <  2:4468614/93‑1 (MQ=37)
                                                cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCAc                                         <  1:2845765/93‑1 (MQ=38)
                                                 tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                        >  1:4313680/1‑93 (MQ=38)
                                                 tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                        >  2:1776741/1‑93 (MQ=38)
                                                  gCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                       <  1:4149294/93‑1 (MQ=38)
                                                    cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                     >  2:2897040/1‑93 (MQ=38)
                                                    cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                     <  1:3982663/93‑1 (MQ=255)
                                                     tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                    >  1:2591041/1‑93 (MQ=38)
                                                        gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTg                                 >  2:2321219/1‑93 (MQ=38)
                                                           gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                              >  1:787263/1‑93 (MQ=37)
                                                            tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt                             <  1:272059/93‑1 (MQ=38)
                                                            tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt                             >  2:2034739/1‑93 (MQ=37)
                                                             gtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTg                            <  1:2499581/93‑1 (MQ=38)
                                                              tACCCTGTCATCAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                           <  2:2487338/93‑1 (MQ=25)
                                                               aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                          >  1:3205416/1‑93 (MQ=37)
                                                               aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                          >  1:1614813/1‑93 (MQ=37)
                                                               aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                          <  2:3339304/93‑1 (MQ=38)
                                                               aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                          >  1:4461096/1‑93 (MQ=37)
                                                                cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                         <  2:2808721/93‑1 (MQ=38)
                                                                cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                         >  1:2650788/1‑93 (MQ=37)
                                                                  cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                       >  2:1095170/1‑93 (MQ=37)
                                                                   tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                      >  2:2974933/1‑93 (MQ=37)
                                                                   tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                      >  2:3085465/1‑93 (MQ=37)
                                                                   tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                      >  2:3359920/1‑93 (MQ=37)
                                                                   tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                      <  2:3524666/93‑1 (MQ=37)
                                                                    gTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt                     >  1:4522276/1‑93 (MQ=37)
                                                                        ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                 >  2:4360257/1‑93 (MQ=25)
                                                                          aaCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTTCCCCCAcccc               >  1:154818/1‑92 (MQ=25)
                                                                           aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa              >  1:4338666/1‑93 (MQ=21)
                                                                           aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa              <  1:4276925/93‑1 (MQ=21)
                                                                            cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccccccccccca             >  1:973283/1‑90 (MQ=25)
                                                                            cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTaa             <  1:3042541/93‑1 (MQ=21)
                                                                             aTGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc            >  1:2320777/1‑93 (MQ=21)
                                                                               gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc          >  1:3838206/1‑93 (MQ=21)
                                                                               gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc          >  2:1538983/1‑93 (MQ=21)
                                                                                 gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTTCCCCCCCCCCCAccccc        >  2:4155983/1‑91 (MQ=18)
                                                                                 gCACCGTTTTCTCCCTCTCCCTGCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctc        <  2:3958926/93‑1 (MQ=18)
                                                                                     cGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCcccc    >  1:178225/1‑93 (MQ=18)
                                                                                      gTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcccc   <  2:3605194/93‑1 (MQ=17)
                                                                                      gTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcccc   >  2:2387873/1‑93 (MQ=18)
                                                                                       ttttCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCa  >  1:3135238/1‑93 (MQ=18)
                                                                                             |                                                                                      
CTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCA  >  NC_000913/3237148‑3237325

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: