Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,555,020 A→G M280T (ATG→ACG)  yjiE ← putative transcriptional regulator LYSR‑type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,555,0200AG99.8% 1355.6 / ‑5.1 415M280T (ATG→ACG) yjiEputative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (192/222);  total (192/223)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 425 total aligned reads displayed.

AGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCG  >  NC_000913/4554928‑4555109
                                                                                            |                                                                                         
actccGAAGCTGGAGAATGTCAGCTAAGCACAATCTGCAGCTCGCGGAGGTCACGCCAGCAGTATTAGAAAACGGGATTCATGAGGGTATTCa                                                                                           >  2:1347309/5‑93 (MQ=255)
  tGCGACACTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTc                                                                                         <  1:1148474/93‑1 (MQ=255)
   gCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTcc                                                                                        <  1:4701478/93‑1 (MQ=255)
    cGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCg                                                                                       <  2:4443721/93‑1 (MQ=255)
    cGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCg                                                                                       <  2:6695130/93‑1 (MQ=255)
     gACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCgg                                                                                      >  1:7623252/1‑93 (MQ=255)
      aCGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGCTTCATGCGGGTATTCGTACGGt                                                                                     >  2:5220972/1‑93 (MQ=255)
       cGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTa                                                                                    <  1:6457923/93‑1 (MQ=255)
       cGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTa                                                                                    <  2:6474464/93‑1 (MQ=255)
       cGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTa                                                                                    <  1:1685819/93‑1 (MQ=255)
          tGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGc                                                                                 >  1:4655179/1‑93 (MQ=255)
            gtgaAGGTCAGCGAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCAt                                                                               <  2:8425777/91‑1 (MQ=255)
             agaATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATa                                                                              >  2:2680688/1‑93 (MQ=255)
             agaATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATa                                                                              >  1:6064385/1‑93 (MQ=255)
              gaATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAg                                                                             <  2:5581820/93‑1 (MQ=255)
               aaTGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAgg                                                                            >  1:5260939/1‑93 (MQ=255)
                aTGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGc                                                                           >  2:3535668/1‑93 (MQ=255)
                 tGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGcc                                                                          >  1:6834094/1‑93 (MQ=255)
                 tGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGcc                                                                          >  2:5384109/1‑93 (MQ=255)
                 tGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGcc                                                                          <  2:8020837/93‑1 (MQ=255)
                   tCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTg                                                                        >  1:2750301/1‑93 (MQ=255)
                     aGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGaa                                                                      <  2:8594078/93‑1 (MQ=255)
                      gCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAAt                                                                     <  2:5824167/93‑1 (MQ=255)
                       cTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATc                                                                    <  1:1559055/93‑1 (MQ=255)
                       cTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATc                                                                    <  2:9235897/93‑1 (MQ=255)
                        tAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCg                                                                   <  2:6314571/93‑1 (MQ=255)
                        tAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCg                                                                   <  1:6296285/93‑1 (MQ=255)
                            cacaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGAt                                                               <  1:7956138/93‑1 (MQ=255)
                            cacaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGAt                                                               >  1:5487778/1‑93 (MQ=255)
                             acaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATc                                                              >  2:8873722/1‑93 (MQ=255)
                             acaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATc                                                              >  2:6609421/1‑93 (MQ=255)
                             acaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATc                                                              <  1:309310/93‑1 (MQ=255)
                             acaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATc                                                              <  2:1284604/93‑1 (MQ=255)
                              caATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCa                                                             >  1:3537970/1‑93 (MQ=255)
                                aTCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCAcc                                                           >  1:4325632/1‑93 (MQ=255)
                                aTCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCAcc                                                           <  2:1931462/93‑1 (MQ=255)
                                aTCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCAcc                                                           <  2:3676926/93‑1 (MQ=255)
                                  ctcCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAg                                                         <  1:885448/93‑1 (MQ=255)
                                   tcCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGt                                                        <  1:3945570/93‑1 (MQ=255)
                                    ccAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGtt                                                       <  1:6755883/93‑1 (MQ=255)
                                    ccAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGGATTCGTCCGGTATGCATAGGCCTGAATCGGGATGACCAGtt                                                       >  2:9034260/1‑93 (MQ=255)
                                      aGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGTCTGAATCGGGATCACCAGTTCg                                                     <  1:5080283/93‑1 (MQ=255)
                                      aGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCg                                                     >  2:231865/1‑93 (MQ=255)
                                       gCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGt                                                    <  1:1780643/93‑1 (MQ=255)
                                          cgcgCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTccc                                                 <  1:9319396/93‑1 (MQ=255)
                                             gcAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGAt                                              <  2:9660076/93‑1 (MQ=255)
                                             gcAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGAt                                              >  1:9206286/1‑93 (MQ=255)
                                             acAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGAt                                              <  1:2230428/92‑1 (MQ=255)
                                                gTTCACGCCAGAAGCGTTCGGCAACGTGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCa                                           <  2:813826/93‑1 (MQ=255)
                                                gTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCa                                           <  2:1088739/93‑1 (MQ=255)
                                                 ttCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAg                                          >  2:6875133/1‑93 (MQ=255)
                                                 ttCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAg                                          <  1:702069/93‑1 (MQ=255)
                                                 ttCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAg                                          >  2:3778032/1‑93 (MQ=255)
                                                 ttCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAg                                          <  2:4980604/93‑1 (MQ=255)
                                                  tCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGt                                         >  1:3345804/1‑93 (MQ=255)
                                                    aCGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTAc                                       <  2:6114948/93‑1 (MQ=255)
                                                     cGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACa                                      >  1:10012068/1‑93 (MQ=255)
                                                      gCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACaa                                     >  1:5747553/1‑93 (MQ=255)
                                                       ccAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAAc                                    >  2:7189974/1‑93 (MQ=255)
                                                       ccAGAAGCGTACGGCAACGGGATTCATGCGGTTATTCGTCCGGTATGCATAGGCCTGAATCGGTATCACCAGTTCGTCCCGATTCAGTACAAc                                    <  1:2634996/93‑1 (MQ=255)
                                                        cAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACg                                   >  2:7041819/1‑93 (MQ=255)
                                                         aGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGa                                  >  1:9872330/1‑93 (MQ=255)
                                                          gAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCCGTTCGTCCCGATTCAGTACAACGAg                                 >  2:9793858/1‑93 (MQ=255)
                                                          gAAGCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTAAAACGAg                                 >  2:947984/1‑93 (MQ=255)
                                                             gCGTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCCCCAGTTCGTCCCGATTCAGTACAACGAGCtt                              <  2:2182971/93‑1 (MQ=255)
                                                              ggTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTc                             <  1:3245009/92‑1 (MQ=255)
                                                               gTTCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTcc                            <  2:2907305/93‑1 (MQ=255)
                                                                ttCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTccc                           >  1:1903975/1‑93 (MQ=255)
                                                                ttCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTccc                           >  1:2895451/1‑93 (MQ=255)
                                                                ttCGGCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTccc                           >  1:8586796/1‑93 (MQ=255)
                                                                    gCAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTg                       >  2:1601549/1‑93 (MQ=255)
                                                                     cAACGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGc                      <  1:7770184/93‑1 (MQ=255)
                                                                       aCGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGa                    >  1:4259666/1‑93 (MQ=255)
                                                                       aCGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGa                    >  2:7900306/1‑93 (MQ=255)
                                                                       aCGGGATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGa                    >  2:7408145/1‑93 (MQ=255)
                                                                          ggATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAAtt                 <  1:4179476/93‑1 (MQ=255)
                                                                           gATTCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAAttt                >  1:205830/1‑93 (MQ=255)
                                                                             ttCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCt              <  1:8141755/93‑1 (MQ=255)
                                                                             ttCATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCt              >  1:5143350/1‑93 (MQ=255)
                                                                               cATGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttg            >  2:4263754/1‑93 (MQ=255)
                                                                                aTGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgt           <  2:3065277/93‑1 (MQ=255)
                                                                                aTGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgt           <  2:399436/93‑1 (MQ=255)
                                                                                aTGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgt           >  2:7705062/1‑93 (MQ=255)
                                                                                aTGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgt           <  1:899884/93‑1 (MQ=255)
                                                                                 tGCGGGTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgtt          <  2:2586923/93‑1 (MQ=255)
                                                                                    gggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGa       >  1:7253215/1‑93 (MQ=255)
                                                                                    gggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGa       <  2:9780979/93‑1 (MQ=255)
                                                                                    gggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGa       <  1:9431855/93‑1 (MQ=255)
                                                                                     ggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGAt      >  2:1402404/1‑93 (MQ=255)
                                                                                     ggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGAt      <  2:4770767/93‑1 (MQ=255)
                                                                                     ggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGAt      <  1:8295692/93‑1 (MQ=255)
                                                                                     ggTATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGATCGTCACGATTCAGTACAACGAGCTTTACGCTGCGCATTTCTTGTTggtt      >  1:4036544/1‑91 (MQ=255)
                                                                                       tATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATgg    >  1:4484210/1‑93 (MQ=255)
                                                                                       tATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATgg    >  2:6212819/1‑93 (MQ=255)
                                                                                       tATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATgg    <  1:974389/93‑1 (MQ=255)
                                                                                       tATTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATgg    >  2:3379635/1‑93 (MQ=255)
                                                                                        aTTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGc   <  1:8498091/93‑1 (MQ=255)
                                                                                        aTTCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGc   >  2:1787633/1‑93 (MQ=255)
                                                                                         ttCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCg  >  1:4906382/1‑93 (MQ=255)
                                                                                         ttCGTCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCg  >  1:7529265/1‑93 (MQ=255)
                                                                                            |                                                                                         
AGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCG  >  NC_000913/4554928‑4555109

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