Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 575,915 T→A 16.9% pseudogene (911/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913575,9150TA16.9% 320.0 / 45.2 130pseudogene (911/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (52/54);  new base A (9/13);  total (61/69)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGA  >  NC_000913/575816‑576014
                                                                                                   |                                                                                                   
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATcc                                                                                                   <  1:236526/100‑1 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAAcc                                                                                                   <  1:2392896/100‑1 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAAcc                                                                                                   <  2:2971971/100‑1 (MQ=255)
    tGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCtt                                                                                                 >  2:1442349/1‑100 (MQ=255)
     gTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTa                                                                                                >  1:1860559/1‑100 (MQ=255)
      tAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTAc                                                                                               <  2:2906128/100‑1 (MQ=255)
      tAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTAc                                                                                               <  2:2072795/100‑1 (MQ=255)
       aaTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTAcc                                                                                              <  2:2893172/100‑1 (MQ=255)
        aTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTAct                                                                                              >  1:3152392/1‑98 (MQ=255)
        aTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCt                                                                                             >  1:6261/1‑100 (MQ=255)
          cAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTACACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCttt                                                                                           >  1:2080708/1‑100 (MQ=255)
           aaCGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTa                                                                                          <  2:2534773/100‑1 (MQ=255)
           aaCGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTa                                                                                          <  2:326780/100‑1 (MQ=255)
           aaCGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTa                                                                                          >  2:941041/1‑100 (MQ=255)
              gAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGAt                                                                                       >  1:1210301/1‑100 (MQ=255)
                  gTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGc                                                                                    >  2:595602/1‑99 (MQ=255)
                  gTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCa                                                                                   >  1:673197/1‑100 (MQ=255)
                   tAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAg                                                                                  <  2:1264523/100‑1 (MQ=255)
                   tAGACATATTTTTGTTGAAGTAATACGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAg                                                                                  <  1:2641336/100‑1 (MQ=255)
                      aCATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTa                                                                               >  2:2045203/1‑100 (MQ=255)
                       cATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTaa                                                                              <  2:1285454/100‑1 (MQ=255)
                        atatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAg                                                                             <  2:506688/100‑1 (MQ=255)
                         tatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGc                                                                            >  2:2158271/1‑100 (MQ=255)
                         tatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGc                                                                            >  2:2953214/1‑100 (MQ=255)
                          atTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCa                                                                           >  1:1014340/1‑100 (MQ=255)
                           tttttGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCaa                                                                          >  2:1478660/1‑100 (MQ=255)
                               tgttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacgga                                                                      >  1:970535/1‑100 (MQ=255)
                               tgttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacgga                                                                      >  1:1362128/1‑100 (MQ=255)
                                gttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCACGTAAGCAacggac                                                                     >  1:2889886/1‑100 (MQ=255)
                                gttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACCCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggac                                                                     >  1:809381/1‑100 (MQ=255)
                                gatgAAGTAATAGATTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggac                                                                     >  2:2301312/3‑100 (MQ=255)
                                  tgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggacgg                                                                   >  1:864362/1‑100 (MQ=255)
                                   gAAGTAATAGGTTGCACCTACATCACCATATTTGACTAAGACCTGATCGCCCCATACTCCAAGATACTTACATATAGATTGCAGGTAAGCAAcggacgga                                                                  >  1:2256712/1‑100 (MQ=255)
                                   gAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCGTACCTTTAGATTGCAGGTAAGCAAcggacgga                                                                  <  2:2588289/100‑1 (MQ=255)
                                    aaGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCAGACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACggacggac                                                                 <  1:1680491/100‑1 (MQ=255)
                                      gTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGc                                                               <  1:1529075/100‑1 (MQ=255)
                                       tAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCa                                                              <  1:3210499/100‑1 (MQ=255)
                                        aaTAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAg                                                             >  1:1730052/1‑100 (MQ=255)
                                         aTAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGa                                                            >  1:1835069/1‑100 (MQ=255)
                                           aGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAcc                                                          <  1:1475271/100‑1 (MQ=255)
                                            ggTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCg                                                         >  1:447463/1‑100 (MQ=255)
                                            ggTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCg                                                         >  2:2535825/1‑100 (MQ=255)
                                             gTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGa                                                        >  1:2400138/1‑100 (MQ=255)
                                                  aCCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCg                                                   <  2:441508/100‑1 (MQ=255)
                                                   ccTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGa                                                  <  1:1539542/100‑1 (MQ=255)
                                                   ccTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGa                                                  >  2:2974053/1‑100 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                <  2:3039789/100‑1 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAGCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                <  1:2760722/100‑1 (MQ=255)
                                                      aCATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACt                                               <  1:2995626/100‑1 (MQ=255)
                                                      aCATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACt                                               <  2:3389290/100‑1 (MQ=255)
                                                      aCATCAACATAGTTGACTAAGTCCTGATCGCCCCATAGGACTTGACCCTTCCCTCTAGATTGCAGCTAAGCAACGGACGGACGCAGACCGAAAGCGAACt                                               <  1:631967/100‑1 (MQ=25)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  1:1047917/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  2:1731235/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  2:2614033/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  2:681600/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              <  2:925356/100‑1 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  1:2039917/1‑100 (MQ=255)
                                                         tCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGat                                            >  1:3285542/1‑100 (MQ=255)
                                                          cAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCCGACCGAAATCGAACTGata                                           >  1:677139/1‑100 (MQ=255)
                                                            aCATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATAtt                                         >  1:1313066/1‑100 (MQ=255)
                                                             cATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtg                                        <  2:79705/100‑1 (MQ=255)
                                                               tatTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgtg                                      <  1:2929815/100‑1 (MQ=255)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     >  2:1242183/1‑100 (MQ=255)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     <  1:1304061/100‑1 (MQ=255)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     <  1:1583168/100‑1 (MQ=255)
                                                                 tttGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCa                                    <  1:1927066/100‑1 (MQ=255)
                                                                  ttGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCaa                                   >  2:2867044/1‑100 (MQ=255)
                                                                   tGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAAc                                  <  2:3326148/100‑1 (MQ=255)
                                                                     aCTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAg                                >  2:171792/1‑100 (MQ=255)
                                                                       tAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCt                              <  1:239907/100‑1 (MQ=255)
                                                                         aGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTc                            <  1:2494998/100‑1 (MQ=255)
                                                                          gTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                           >  2:711804/1‑100 (MQ=255)
                                                                          gTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                           >  1:2590610/1‑100 (MQ=255)
                                                                           tCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGa                          <  1:2740253/100‑1 (MQ=255)
                                                                            ncTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGAAAGCAACTGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGc                               >  1:2798486/2‑94 (MQ=37)
                                                                            ccTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGaa                         >  1:404179/1‑100 (MQ=255)
                                                                             cTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAg                        <  2:2244880/100‑1 (MQ=255)
                                                                               gATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtt                      >  2:3106707/1‑100 (MQ=255)
                                                                                 tCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtttt                    <  1:538928/100‑1 (MQ=255)
                                                                                  ctccccATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTg                   <  2:1903478/98‑1 (MQ=255)
                                                                                  cGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTg                   <  1:1442343/100‑1 (MQ=255)
                                                                                   gCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg                  >  1:2096803/1‑100 (MQ=255)
                                                                                    ccccATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTggg                 >  2:2131958/1‑100 (MQ=255)
                                                                                    ccccATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTggg                 <  1:3389015/100‑1 (MQ=255)
                                                                                     cccATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGc                >  2:1950884/1‑100 (MQ=255)
                                                                                      ccATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCt               <  2:173917/100‑1 (MQ=255)
                                                                                      ccATACTCCAAGATCCTTACCTCTAGATTGCAGGTAAGCAACGGACTGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg                  <  1:2047448/97‑1 (MQ=255)
                                                                                       cATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGGTATTGTGCAACAGCTTCGAAGTTTTGGGCtt              <  1:2470765/100‑1 (MQ=255)
                                                                                       cATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg                  <  1:164973/96‑1 (MQ=255)
                                                                                          aCTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttat           >  2:103141/1‑100 (MQ=255)
                                                                                          aCTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGAAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttat           <  2:2339843/100‑1 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          >  2:2848690/1‑100 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          <  1:3035318/100‑1 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          <  2:1362134/100‑1 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttat           <  2:1992574/99‑1 (MQ=255)
                                                                                            tCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatta         <  2:3044012/100‑1 (MQ=255)
                                                                                             ccAAGAACCTTACCTTTAATTTGCAGGTAAGCAACGGACGGCCGCAGATCGAAATTGAACTGATATTGTGCAACAGCTTCGGAGTTTTGGGCTttattan        <  1:1475955/100‑2 (MQ=255)
                                                                                              cAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGc       <  1:1892491/100‑1 (MQ=255)
                                                                                               aaGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCa      >  1:2119211/1‑100 (MQ=255)
                                                                                                aGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCaa     >  1:3037616/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGA  >  NC_000913/575816‑576014

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