Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,646,847 T→C intergenic (+110/+327) flxA → / ← intK Qin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,646,8470TC99.3% 487.6 / ‑4.8 138intergenic (+110/+327)flxA/intKQin prophage; uncharacterized protein/pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (78/59);  total (78/60)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 138 total aligned reads displayed.

TGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTT  >  NC_000913/1646751‑1646942
                                                                                                |                                                                                               
tgtgGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAAt                                                                                              <  1:2284540/100‑1 (MQ=255)
tgtgGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAAt                                                                                              <  2:2108173/100‑1 (MQ=255)
 gtgGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATg                                                                                             >  1:2334735/1‑100 (MQ=255)
 gtgGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATg                                                                                             >  1:1498913/1‑100 (MQ=255)
  tgGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATgg                                                                                            <  1:1128136/100‑1 (MQ=255)
   ggCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATTGa                                                                                           >  2:1825902/1‑100 (MQ=255)
     cccTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACg                                                                                         <  1:730012/100‑1 (MQ=255)
      ccTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGTTAGTTAACATTACAAGATATAAGCAATGGACGc                                                                                        <  2:2016267/100‑1 (MQ=255)
       cTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCa                                                                                       >  2:1591184/1‑100 (MQ=255)
       cTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCa                                                                                       <  2:1583731/100‑1 (MQ=255)
        tCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCa                                                                                       >  2:2541574/1‑99 (MQ=255)
         cATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACt                                                                                     >  2:2658029/1‑100 (MQ=255)
          aTCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTc                                                                                    >  2:481889/1‑100 (MQ=255)
           tCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTcc                                                                                   <  1:1789414/100‑1 (MQ=255)
           tCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTcc                                                                                   <  1:64184/100‑1 (MQ=255)
           tCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTcc                                                                                   >  1:544662/1‑100 (MQ=255)
           tCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTcc                                                                                   <  1:3027507/100‑1 (MQ=255)
                ccAGTTTTCGCCATCAGTCTTTTCTTTAATGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAAtt                                                                              <  2:1842104/100‑1 (MQ=255)
                ccACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAAtt                                                                              <  1:402496/100‑1 (MQ=255)
                ccACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAAtt                                                                              >  1:1168504/1‑100 (MQ=255)
                 cACGTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTGCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAAtt                                                                              <  2:2644433/99‑1 (MQ=255)
                  aCTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAg                                                                            <  2:1286304/100‑1 (MQ=255)
                  aCTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAg                                                                            <  1:47139/100‑1 (MQ=255)
                  aCTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAg                                                                            <  1:2395651/100‑1 (MQ=255)
                   cTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGt                                                                           >  2:1767415/1‑100 (MQ=255)
                    ttttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTc                                                                          <  2:912587/100‑1 (MQ=255)
                      ttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTa                                                                        >  1:2754303/1‑100 (MQ=255)
                      ttCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTa                                                                        >  2:1324606/1‑100 (MQ=255)
                       tCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAAAATTAAAAGATATAAGCAATGGACGAACTCCCAATTAGTCTAt                                                                       >  2:2151026/1‑100 (MQ=255)
                         gCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTAttt                                                                     <  2:739620/100‑1 (MQ=255)
                            aTCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTaaa                                                                  >  2:1727225/1‑100 (MQ=255)
                            aTCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTaaa                                                                  >  1:1735257/1‑100 (MQ=255)
                            aTCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGAGATAAGCAATGGACGCACTCCCAATTAGTCTATTTaaa                                                                  >  2:2060576/1‑100 (MQ=255)
                             tCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAAt                                                                 >  1:1642677/1‑100 (MQ=255)
                              cAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATc                                                                <  1:2373397/100‑1 (MQ=255)
                              cAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATc                                                                >  1:753677/1‑100 (MQ=255)
                              cAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATc                                                                >  1:966102/1‑100 (MQ=255)
                               aGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCg                                                               <  1:124832/100‑1 (MQ=255)
                                gCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGc                                                              <  2:1460195/100‑1 (MQ=255)
                                   ttttCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCAc                                                           >  1:2784662/1‑100 (MQ=255)
                                     tacTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAAGCGCCACg                                                          <  1:1602933/97‑1 (MQ=255)
                                      tCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAg                                                        <  2:1494172/100‑1 (MQ=255)
                                         ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGttt                                                     >  1:525197/1‑100 (MQ=255)
                                         ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGttt                                                     <  2:438737/100‑1 (MQ=255)
                                         ttAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGttt                                                     >  1:1206704/1‑100 (MQ=255)
                                          tAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGATTa                                                    >  2:572039/1‑100 (MQ=255)
                                           aaaGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTaa                                                   >  1:96120/1‑100 (MQ=255)
                                           aaaGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTaa                                                   >  2:1224250/1‑100 (MQ=255)
                                            aaGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAAc                                                  >  2:1903908/1‑100 (MQ=255)
                                             aGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACt                                                 <  2:332763/100‑1 (MQ=255)
                                             aGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACt                                                 >  1:252698/1‑100 (MQ=255)
                                              gACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTg                                                <  1:2132436/100‑1 (MQ=255)
                                              gACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTg                                                <  2:3047560/100‑1 (MQ=255)
                                                cATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGAc                                              >  2:1572296/1‑100 (MQ=255)
                                                 atatTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACa                                             <  2:763807/100‑1 (MQ=255)
                                                 atatTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACa                                             <  2:662821/100‑1 (MQ=255)
                                                 atatTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACa                                             >  1:2518835/1‑100 (MQ=255)
                                                  tattatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCGCTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACaa                                            >  2:2958489/1‑100 (MQ=255)
                                                  tattatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACaa                                            >  1:645936/1‑100 (MQ=255)
                                                   attatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAAc                                           >  2:1361487/1‑100 (MQ=255)
                                                    ttatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATATGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAAcc                                          <  2:2142739/100‑1 (MQ=255)
                                                    ttatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAAcc                                          <  2:438872/100‑1 (MQ=255)
                                                       tCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCAt                                       >  1:1970758/1‑100 (MQ=255)
                                                        cTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATg                                      <  1:2829465/100‑1 (MQ=255)
                                                             tATCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCaa                                 >  2:2339871/1‑100 (MQ=255)
                                                               tCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAAtt                               <  2:1935409/100‑1 (MQ=255)
                                                               tCATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAAtt                               >  1:2181824/1‑100 (MQ=255)
                                                                cATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTa                              >  2:951592/1‑100 (MQ=255)
                                                                cATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTa                              <  2:822448/100‑1 (MQ=255)
                                                                cATTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTa                              <  2:96831/100‑1 (MQ=255)
                                                                 aTTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTAt                             >  1:1476667/1‑100 (MQ=255)
                                                                 aTTTCTGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTAt                             <  2:2800665/100‑1 (MQ=255)
                                                                      tGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAAtt                        <  2:587304/100‑1 (MQ=255)
                                                                      tGATAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAAtt                        >  1:2980983/1‑100 (MQ=255)
                                                                        aTAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGc                      >  1:1157243/1‑100 (MQ=255)
                                                                        aTAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGc                      <  2:2980942/100‑1 (MQ=255)
                                                                         tAGTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCa                     <  1:928700/100‑1 (MQ=255)
                                                                           gTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAAc                   >  1:777495/1‑100 (MQ=255)
                                                                           gTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAAc                   >  2:893786/1‑100 (MQ=255)
                                                                           gTTAACATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAAc                   >  1:2192184/1‑100 (MQ=255)
                                                                            ttAACATTACAAGATATAAGCAATGGACGCACCCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACt                  <  1:242649/100‑1 (MQ=255)
                                                                              aaCATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTAt                >  2:1752673/1‑100 (MQ=255)
                                                                               aCATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTAtt               <  2:1206692/100‑1 (MQ=255)
                                                                                cATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTAttt              >  1:1212345/1‑100 (MQ=255)
                                                                                cATTACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTAttt              >  2:712113/1‑100 (MQ=255)
                                                                                  ttACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCt            >  2:2532635/1‑100 (MQ=255)
                                                                                   tACAAGATATAAGCAATGGACGCACTGCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTg           <  2:2029309/100‑1 (MQ=255)
                                                                                   tACAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTg           <  1:540964/100‑1 (MQ=255)
                                                                                    aCAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt          <  1:1825924/100‑1 (MQ=255)
                                                                                    aCAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt          >  2:1774843/1‑100 (MQ=255)
                                                                                    aCAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt          <  2:617531/100‑1 (MQ=255)
                                                                                    aCAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt          >  1:2866888/1‑100 (MQ=255)
                                                                                    aCAAGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGt          >  2:1381732/1‑100 (MQ=255)
                                                                                      aaGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAg        <  1:2670514/100‑1 (MQ=255)
                                                                                      aaGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAg        <  2:286634/100‑1 (MQ=255)
                                                                                       aGATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGt       >  2:1303434/1‑100 (MQ=255)
                                                                                        gATATAAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTc      >  1:2618443/1‑100 (MQ=255)
                                                                                          tataAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCAc    >  1:2125972/1‑100 (MQ=255)
                                                                                            taAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACtt  >  1:1959009/1‑100 (MQ=255)
                                                                                            taAGCAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACtt  >  1:1669684/1‑100 (MQ=255)
                                                                                                |                                                                                               
TGTGGCCCTCATCGGGCCACTTTTCGCCATCAGCCTTTTCTTTAAAGACATATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTT  >  NC_000913/1646751‑1646942

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