Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 575,915 T→A 41.7% 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,9150TA41.7% 62.9 / 101.7 93pseudogene (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 (26/27);  new base A (16/22);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGA  >  NC_000913/575816‑576014
                                                                                                   |                                                                                                   
aTTTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAt                                                                                                     >  2:86682/1‑100 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGTTCGCCCCATACTCCAAGAAcc                                                                                                   <  2:1004790/100‑1 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATcc                                                                                                   <  2:2233977/100‑1 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATcc                                                                                                   <  1:2194486/100‑1 (MQ=255)
  tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAAcc                                                                                                   <  1:1715925/100‑1 (MQ=255)
   ttGTAATCAACGAAAGTAGACATATTTTTGTTGAATTAATATGTCGCACCTACATCAACATATTTGACTAAGTCCTGATCCCCCCATACTCCAAGATCCt                                                                                                  >  2:2281500/1‑100 (MQ=38)
     gTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTa                                                                                                <  2:1441444/100‑1 (MQ=255)
     gTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTa                                                                                                <  1:2325188/100‑1 (MQ=255)
       aaTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTAcc                                                                                              >  1:540839/1‑100 (MQ=255)
         tCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAaccttacctt                                                                                            >  2:2244122/1‑100 (MQ=255)
           aaCGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTa                                                                                          >  1:631923/1‑100 (MQ=255)
           aaCGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTa                                                                                          >  2:219898/1‑100 (MQ=255)
             cGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGACCCATACTCCAAGAACCTTACCttt                                                                                           >  2:2180907/1‑97 (MQ=255)
              aaaaGTAGACATATTTTTGTTGAAGTAATAGATTGCACCTACATCAACATATTTGACTAAGTACTCATCGCCCCAAACTCTAAGAACCTTAACTTTAGac                                                                                       >  2:1907458/2‑99 (MQ=255)
                  gTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGc                                                                                    <  2:2394641/99‑1 (MQ=255)
                    aGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAgg                                                                                 <  1:1834051/100‑1 (MQ=255)
                    aGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAgg                                                                                 <  2:1156772/100‑1 (MQ=255)
                      aCATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTa                                                                               <  2:641512/100‑1 (MQ=255)
                       cATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTaa                                                                              >  2:2155328/1‑100 (MQ=255)
                       cATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTaa                                                                              >  1:1941961/1‑100 (MQ=255)
                       cATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTaa                                                                              >  1:184331/1‑100 (MQ=255)
                         tatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGc                                                                            <  2:631920/100‑1 (MQ=255)
                         tatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGc                                                                            >  1:1775753/1‑100 (MQ=255)
                         tatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGc                                                                            <  1:1530578/100‑1 (MQ=255)
                          atTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCa                                                                           <  2:2311371/100‑1 (MQ=255)
                            ttttGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGAAac                                                                         >  1:621833/1‑100 (MQ=255)
                            ttttGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAac                                                                         <  2:814902/100‑1 (MQ=255)
                             tttGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacg                                                                        <  2:602707/100‑1 (MQ=255)
                                gttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggac                                                                     >  2:1852339/1‑100 (MQ=255)
                                gttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggac                                                                     >  2:2161306/1‑100 (MQ=255)
                                 ttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggacg                                                                    >  1:2058235/1‑100 (MQ=255)
                                  tgAAGTAATAGGTTGCACCTACATCAAAATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAacggacgg                                                                   >  2:1166625/1‑100 (MQ=255)
                                    aaGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACggacggac                                                                 >  2:1151665/1‑100 (MQ=255)
                                      gTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGc                                                               <  1:1729936/100‑1 (MQ=255)
                                      gTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGc                                                               >  1:280226/1‑100 (MQ=255)
                                         aTAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGa                                                            >  2:292042/1‑100 (MQ=255)
                                         aTAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGa                                                            <  2:1047504/100‑1 (MQ=255)
                                          tAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAc                                                           >  1:999322/1‑100 (MQ=255)
                                           aGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAcc                                                          <  2:151709/100‑1 (MQ=255)
                                            ggTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAcc                                                          >  2:458666/1‑99 (MQ=255)
                                            ggTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAcc                                                          >  1:87035/1‑99 (MQ=255)
                                               tGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTTCAAGACCCTTACCTTTAGTTTGCATGTAATCAACGGACGGACGCAGa                                                            <  2:1047626/94‑1 (MQ=255)
                                                  aCCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCg                                                   <  1:1678319/100‑1 (MQ=255)
                                                  aCCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCg                                                   <  1:1337614/100‑1 (MQ=255)
                                                   ccTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGa                                                  >  2:1171121/1‑100 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                >  2:2275318/1‑100 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                <  2:1217592/100‑1 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                <  2:650017/100‑1 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                <  1:373320/100‑1 (MQ=255)
                                                     tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                                >  2:1645720/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  2:2339068/1‑100 (MQ=255)
                                                       cATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTg                                              >  2:597285/1‑100 (MQ=255)
                                                          cAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGata                                           >  2:1016697/1‑100 (MQ=255)
                                                          cAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACAGAAATCGAACTGcta                                           >  2:1486549/1‑100 (MQ=255)
                                                           aaCATATTTGACTAAGTCCTGATCGCCCCATATTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGatat                                          >  2:1781538/1‑100 (MQ=255)
                                                             cATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtg                                        <  1:1763107/100‑1 (MQ=255)
                                                              atatTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgt                                       <  2:2269974/100‑1 (MQ=255)
                                                               tatTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgtg                                      <  2:1495927/100‑1 (MQ=255)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     <  1:1938587/100‑1 (MQ=255)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGAACGCAACGGACGGCCCCAGCCCGAAATCGAACTGCTATtgtcc                                     >  2:2055055/1‑98 (MQ=37)
                                                                atTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     <  1:1064266/100‑1 (MQ=255)
                                                                 tttGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                     <  2:1634024/99‑1 (MQ=255)
                                                                 tttGACTAAGTCCTGATCGCCCCATACTCCAAGACCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCa                                    <  1:1099629/100‑1 (MQ=255)
                                                                   tGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAAc                                  <  1:1236403/100‑1 (MQ=255)
                                                                   tGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAAc                                  <  2:1204987/100‑1 (MQ=255)
                                                                   tGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAAc                                  >  2:2368406/1‑100 (MQ=255)
                                                                    gACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACa                                 <  2:1556488/100‑1 (MQ=255)
                                                                      cTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGc                               >  1:1120866/1‑100 (MQ=255)
                                                                       tAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCt                              <  2:833598/100‑1 (MQ=255)
                                                                        aaGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTACAGGTAAGCAACGGACACACGCAGACAGAAATCGAACTGATATATCACAACAGCt                              >  2:2238964/1‑99 (MQ=37)
                                                                          gTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                           <  1:305102/100‑1 (MQ=255)
                                                                          gTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                           <  2:1539146/100‑1 (MQ=255)
                                                                           tCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                           <  1:669772/99‑1 (MQ=255)
                                                                                  cGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGGCGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtttt                    >  1:640802/1‑99 (MQ=255)
                                                                                   gCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg                  <  2:1249824/100‑1 (MQ=255)
                                                                                   gCCCCATACTCAAAGAACCTTACCTTTAGATTGCAGGTAAGCAAAGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg                  <  1:107840/100‑1 (MQ=255)
                                                                                    ccccATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtttt                    >  1:799115/1‑97 (MQ=255)
                                                                                      ccATAATCCAAGAACCTTACCTTTAGATTGCAGGTAATCAACGGACGGACCCAGACCGACATCGAACTTATCTTGCGCCACCGCATCGAAgctccac                  >  2:2190889/1‑91 (MQ=19)
                                                                                       cctACTCCAAGAATCTTTACTTTAGATTTCATGTAATCAACGGACTGACGCAGACCGTAATCGAACTGGTATTGTGCTAGAGCTTCGAAGTTTTGGGc                <  2:748757/96‑1 (MQ=16)
                                                                                         tACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTtta            <  1:2244148/100‑1 (MQ=255)
                                                                                          aCTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttat           <  2:1930219/100‑1 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          <  1:51671/100‑1 (MQ=255)
                                                                                           cTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          >  1:1261812/1‑100 (MQ=255)
                                                                                           cTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatt          <  2:2109507/100‑1 (MQ=255)
                                                                                            tCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatta         <  1:219901/100‑1 (MQ=255)
                                                                                            tCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatta         <  1:2329046/100‑1 (MQ=255)
                                                                                            tCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatta         <  2:2061473/100‑1 (MQ=255)
                                                                                             ccAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAg        >  2:1404611/1‑100 (MQ=255)
                                                                                             ccAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAg        >  1:765559/1‑100 (MQ=255)
                                                                                             ccAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAg        >  1:1888441/1‑100 (MQ=255)
                                                                                              cAAGATCCTTACCTTTCGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTtta            <  2:2361025/95‑1 (MQ=255)
                                                                                              cAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTCTATTAGc       >  2:808751/1‑100 (MQ=255)
                                                                                               aaGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGATCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCa      >  1:637568/1‑100 (MQ=255)
                                                                                                 gATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAAc    <  1:1171132/100‑1 (MQ=255)
                                                                                                  aTCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACg   <  1:2344763/100‑1 (MQ=255)
                                                                                                   tCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGa  >  1:1798193/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: