Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,551 A→G 49.8% intergenic (‑14/‑179) yjgM ← / → yjgN GNAT family putative N‑acetyltransferase/DUF898 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,5510AG49.8% 11.0 / 241.1 163intergenic (‑14/‑179)yjgM/yjgNGNAT family putative N‑acetyltransferase/DUF898 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/53);  new base G (47/34);  total (76/87)
Fisher's exact test for biased strand distribution p-value = 4.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 181 total aligned reads displayed.

ACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTACTGATA  >  NC_000913/4479452‑4479648
                                                                                                   |                                                                                                 
  cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaa                                                                                                 <  1:3111505/100‑1 (MQ=255)
  cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaa                                                                                                 >  2:604465/1‑100 (MQ=255)
   tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaatgaaa                                                                                                >  1:4094151/1‑100 (MQ=255)
   tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaa                                                                                                <  1:4063088/100‑1 (MQ=255)
    gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaatgaaaa                                                                                               <  1:3896630/100‑1 (MQ=255)
    gACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaataaaaa                                                                                               <  2:1691130/100‑1 (MQ=255)
      cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaatgaaaata                                                                                             >  1:4130802/1‑100 (MQ=255)
      cGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCaatgaaaata                                                                                             >  1:3323127/1‑100 (MQ=255)
          tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGt                                                                                         <  2:3682527/100‑1 (MQ=255)
           gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTa                                                                                        <  2:3136292/100‑1 (MQ=255)
           gACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTa                                                                                        <  1:1320138/100‑1 (MQ=255)
            aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa                                                                                       <  1:279294/100‑1 (MQ=255)
            aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTaa                                                                                       <  1:2041861/100‑1 (MQ=255)
             cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAt                                                                                      >  1:3870742/1‑100 (MQ=255)
              gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAAtt                                                                                     >  1:1142739/1‑100 (MQ=255)
               cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATta                                                                                    <  1:11488/100‑1 (MQ=255)
                  gctctAGCAGGATTGTTCTGCAGCGTCAGTGGACGTATTACTAGTCATTCAGACGTAATGTTATTCATAGTTTGACCCAATACACATAAAGTAATtata                                                                                  <  2:828344/95‑1 (MQ=255)
                  gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATAc                                                                                 >  1:4085632/1‑100 (MQ=255)
                      ctgCAGGATAATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACAttt                                                                             >  1:776680/3‑100 (MQ=255)
                       aGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTg                                                                            <  1:1232300/100‑1 (MQ=255)
                         cAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGtt                                                                          <  1:1713415/100‑1 (MQ=255)
                         cAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGtt                                                                          <  1:1402395/100‑1 (MQ=255)
                          aGGATTATCCTGTAGCGTCAGGCGACGCATTACAGATGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGttt                                                                         <  2:2184584/100‑1 (MQ=255)
                          aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGttt                                                                         <  1:2134487/100‑1 (MQ=255)
                          aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGttt                                                                         <  2:3113264/100‑1 (MQ=255)
                          aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGttt                                                                         <  2:4063915/100‑1 (MQ=255)
                          aGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGAGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGttt                                                                         <  2:3237790/100‑1 (MQ=255)
                           ggAATATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAAGTTTAGTCATAGTGTAACCCAATAAAAATAAAGTAAGTATACATTTGTTTa                                                                        >  2:1372540/1‑100 (MQ=255)
                             atgaTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCTCAATGTTATTCATATTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAAt                                                                      <  2:2433807/97‑1 (MQ=255)
                              ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATa                                                                     <  2:285125/100‑1 (MQ=255)
                              ttATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATa                                                                     <  1:3272106/100‑1 (MQ=255)
                              ttATCCTGCAGCATCAGTCTACGCATTACAGGGAATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAAt                                                                      <  2:3046274/99‑1 (MQ=255)
                                aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATAcc                                                                   <  1:3125493/100‑1 (MQ=255)
                                aTCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATAcc                                                                   <  1:1473571/100‑1 (MQ=255)
                                  ccTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc                                                                 <  1:3577399/100‑1 (MQ=255)
                                  ccTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCAc                                                                 >  1:532357/1‑100 (MQ=255)
                                    tGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCActc                                                               >  1:505722/1‑100 (MQ=255)
                                       aGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCActctct                                                            <  1:2129678/100‑1 (MQ=255)
                                         cGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCtcttc                                                          >  1:1372418/1‑100 (MQ=255)
                                           tCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTcttctt                                                        <  1:1867001/100‑1 (MQ=255)
                                            cAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTta                                                       <  1:2874544/100‑1 (MQ=255)
                                            cAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTta                                                       >  1:1256124/1‑100 (MQ=255)
                                            cAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTta                                                       <  1:1587361/100‑1 (MQ=255)
                                             aGGCGACGCGTTACAGGTGATTGCTGCTGATTTTTATTCATAGGTTTACCCAACAACATTAAAGTAATGATAGATTTGTTTAATACCACGCTCTTCTta                                                       <  2:1014196/99‑1 (MQ=255)
                                              ggCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTtacg                                                     >  1:140600/1‑100 (MQ=255)
                                              ggCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTtacg                                                     >  1:2601141/1‑100 (MQ=255)
                                              ggCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTATCTTCTtac                                                      >  2:3530492/1‑99 (MQ=255)
                                                cGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTtacgta                                                   >  1:2360046/1‑100 (MQ=255)
                                                  aCGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTtacttacg                                                 >  2:910180/1‑100 (MQ=255)
                                                  aCGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTtacgtacg                                                 <  2:4077510/100‑1 (MQ=255)
                                                  aCGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTcacgtacg                                                 <  1:2735861/100‑1 (MQ=255)
                                                     cATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCaa                                              >  1:1007560/1‑100 (MQ=255)
                                                      aTTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATCTGTTTAATACCACTCTCTTCTTACGTACGCAAg                                             >  2:1806804/1‑100 (MQ=255)
                                                       ttACAAGTGATTATGACGAAATTTTATTCATAAGTTAACCCAGTATAACTAAAGTAATCATACATTTGTTTAATACCACACTCTTCTTACATACGCa                                               <  2:3485303/97‑1 (MQ=255)
                                                        tACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCg                                           <  1:2329503/100‑1 (MQ=255)
                                                        tACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCg                                           <  1:3109653/100‑1 (MQ=255)
                                                        tACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCg                                           <  2:2250876/100‑1 (MQ=255)
                                                          cAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGAc                                         >  2:3184730/1‑100 (MQ=255)
                                                           aGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACt                                        <  2:143546/100‑1 (MQ=255)
                                                            ggTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACtt                                       >  2:1845397/1‑100 (MQ=255)
                                                            ggTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGAc                                         >  1:2638289/1‑98 (MQ=255)
                                                             gTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTta                                      <  1:3930436/100‑1 (MQ=255)
                                                             gTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTta                                      <  2:2690808/100‑1 (MQ=255)
                                                                aTTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATaa                                   <  2:2127841/100‑1 (MQ=255)
                                                                aTTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATaa                                   >  2:3609190/1‑100 (MQ=255)
                                                                  tGCGGCGCAATGTTATTCATAGTTTAACCCTATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAg                                 <  2:754705/100‑1 (MQ=255)
                                                                  tGCGGCGCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAg                                 >  2:1212581/1‑100 (MQ=255)
                                                                    cggcgCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGta                               >  1:1678942/1‑100 (MQ=255)
                                                                    cggcgCAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGta                               >  2:1988744/1‑100 (MQ=255)
                                                                    cggcgCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGta                               >  1:1650899/1‑100 (MQ=255)
                                                                     ggcgCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGtaa                              >  2:3326674/1‑100 (MQ=255)
                                                                     ggcgCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGtaa                              <  1:2262640/100‑1 (MQ=255)
                                                                     ggcgCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGtaa                              >  2:2169272/1‑100 (MQ=255)
                                                                       cgcAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCACTTCTTACGTACGCAAGCGACTTATAAAGtaata                            >  2:1633261/1‑100 (MQ=255)
                                                                       cgcAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGtaata                            <  1:3597145/100‑1 (MQ=255)
                                                                        gcAATGTTATTCATAGTTTAACCCAAAGAAAATCAAGTAATAATCCATATGTTTAATAACACTCTCCTCTTACGTAAGCAACCGACTTATACAGcaa                              <  2:3449340/97‑1 (MQ=255)
                                                                         cAATGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATaaa                          >  2:2843511/1‑100 (MQ=255)
                                                                         cAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATaaa                          >  1:860105/1‑100 (MQ=255)
                                                                          aaTGTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTACTTACGTACGCAAGCGACTTATAAAGTAATAAAg                         <  2:37971/100‑1 (MQ=255)
                                                                          aaTGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAg                         <  1:1640560/100‑1 (MQ=255)
                                                                             gTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGttt                      >  1:1391291/1‑100 (MQ=255)
                                                                             gTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGttt                      >  1:3841145/1‑100 (MQ=255)
                                                                             gTTATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGttt                      <  2:129839/100‑1 (MQ=255)
                                                                             gTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGttt                      <  2:514747/100‑1 (MQ=255)
                                                                              ttATTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGGACGCAAGCGACTTATAAAGTAATAAAGttt                      <  2:2979148/99‑1 (MQ=255)
                                                                                aTTCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGttttt                    >  1:1048053/1‑99 (MQ=255)
                                                                                aTTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTg                   <  2:1970337/100‑1 (MQ=255)
                                                                                aTTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATAACACTCTCTTCTTACGTACGCAACCGACTTATAAAGTAATAAAGt                        >  2:488136/1‑95 (MQ=255)
                                                                                  tCATAGTTTAACCCAATGAAAATACAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGcc                 <  2:3293394/100‑1 (MQ=255)
                                                                                  tCATAGTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGcc                 <  1:3175186/100‑1 (MQ=255)
                                                                                       gTTTAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCt            >  1:762746/1‑100 (MQ=255)
                                                                                          tAACCCAATGAAAATAAATTAATTAGAGATTTGTTTAATACCTCTCTCTTCTTACGGACGCAAGTTACATATTAATTAATAAAGTTTTTGCCAGCCTAt          <  2:1201657/99‑1 (MQ=255)
                                                                                          tAACCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTAtt         >  2:3980607/1‑100 (MQ=255)
                                                                                          tAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTAtt         >  2:646858/1‑100 (MQ=255)
                                                                                            aCCCAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTAc       >  2:3844143/1‑100 (MQ=255)
                                                                                            aCCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTAc       >  1:703299/1‑100 (MQ=255)
                                                                                             cccAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTACt      >  1:1207052/1‑100 (MQ=255)
                                                                                               cAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTACTGa    >  1:1819852/1‑100 (MQ=255)
                                                                                               cAATGAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAATACTAAAGATTTTGCCAGCCTATTACt      >  1:3442898/1‑98 (MQ=255)
                                                                                               cAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTACTGa    >  1:3231476/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
ACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCCAATAAAAATAAAGTAATTATACATTTGTTTAATACCACTCTCTTCTTACGTACGCAAGCGACTTATAAAGTAATAAAGTTTTTGCCAGCCTATTACTGATA  >  NC_000913/4479452‑4479648

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