Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2770CT46.9% 34.9 / 116.8 98pseudogene (100/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base T (21/25);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAG  >  NC_000913/566147‑566415
                                                                                                                                  |                                                                                                                                          
tCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcgcg                                                                                                                                    <  1:320618/139‑1 (MQ=255)
   aTTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCa                                                                                                                                 >  1:1138878/1‑139 (MQ=255)
      ttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              <  2:1052000/139‑1 (MQ=255)
      ttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                                              <  2:1138878/139‑1 (MQ=255)
      ttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCACGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                                              <  2:1024526/139‑1 (MQ=255)
        ttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            <  1:65462/139‑1 (MQ=255)
        ttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            <  1:844405/139‑1 (MQ=255)
             tGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCt                                                                                                                       <  1:1212651/139‑1 (MQ=255)
                       cACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             <  2:425578/139‑1 (MQ=255)
                         cgaTCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTg                                                                                                           >  2:1445232/3‑139 (MQ=255)
                           aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             <  1:1244372/135‑1 (MQ=255)
                           aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             >  2:1244372/1‑135 (MQ=255)
                            tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAGTCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTCCCAGGCATCTTTTCGCGTCAc                                                                                                                                <  2:745017/115‑1 (MQ=255)
                            tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAc                                                                                                                                >  1:745017/1‑115 (MQ=255)
                            tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         <  1:980155/122‑1 (MQ=255)
                            tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         >  2:980155/1‑122 (MQ=255)
                            tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGc                                                                                                        >  1:1526490/1‑139 (MQ=255)
                                    aaaCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAg                                                                                                <  1:1043709/139‑1 (MQ=255)
                                         cGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGac                                                                                           <  1:1445232/139‑2 (MQ=255)
                                          gCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGctgact                                                                                          >  1:146140/1‑134 (MQ=255)
                                               gTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCa                                                                                                               >  1:487625/1‑113 (MQ=255)
                                               gTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCa                                                                                                               <  2:487625/113‑1 (MQ=255)
                                                 aCTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                                                                   <  2:982543/139‑1 (MQ=255)
                                                        atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              <  2:605254/86‑1 (MQ=255)
                                                        atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              >  1:605254/4‑89 (MQ=255)
                                                          cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            >  2:111709/2‑89 (MQ=255)
                                                          cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            <  1:111709/88‑1 (MQ=255)
                                                          cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                                          >  2:411580/2‑139 (MQ=255)
                                                            cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCTCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcac                                                                                                                            >  1:934350/1‑87 (MQ=255)
                                                            cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcac                                                                                                                            <  2:934350/87‑1 (MQ=255)
                                                             gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                       >  2:272038/1‑139 (MQ=255)
                                                             gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                       >  1:629648/1‑139 (MQ=255)
                                                             gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   >  2:1225461/1‑79 (MQ=255)
                                                             gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   <  1:1225461/79‑1 (MQ=255)
                                                              cTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCg                                                                      >  2:405230/1‑139 (MQ=255)
                                                              cTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCg                                                                      >  1:1194837/1‑139 (MQ=255)
                                                                    cACGACATAATGCAGGCCATCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                <  1:1033693/139‑1 (MQ=255)
                                                                      cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAc                                                                                                              <  1:1414485/91‑1 (MQ=255)
                                                                      cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAc                                                                                                              >  2:1414485/1‑91 (MQ=255)
                                                                        aCATAATACAGTCCTTCACGCTTCATACGCGGGTCAAGGTTGGCAAACTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGCGCCATGTAAAATGACTTTATGTCATCGAAACCACAGAGCCGGaact                                                            >  2:1189199/1‑137 (MQ=255)
                                                                             aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                                              >  2:482476/1‑68 (MQ=255)
                                                                             aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                                              <  1:482476/68‑1 (MQ=255)
                                                                             aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   >  2:1536830/1‑63 (MQ=255)
                                                                             aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   <  1:1536830/63‑1 (MQ=255)
                                                                             aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                       >  1:31835/1‑139 (MQ=255)
                                                                                   gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         <  1:1222080/67‑1 (MQ=255)
                                                                                   gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         >  2:1222080/1‑67 (MQ=255)
                                                                                   gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATcc                                                 >  1:1348270/1‑139 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              <  1:1304277/62‑1 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              >  2:1304277/1‑62 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            <  2:192852/64‑1 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            >  1:192852/1‑64 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               <  1:526545/93‑1 (MQ=255)
                                                                                   gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGAACCTGGGCCATGTAAGc                                                                                               >  2:526545/1‑93 (MQ=255)
                                                                                    ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGct                                                                                              <  1:946989/93‑2 (MQ=255)
                                                                                    ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGct                                                                                              >  2:946989/1‑92 (MQ=255)
                                                                                     cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                               <  1:354034/139‑1 (MQ=255)
                                                                                      ttCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                              <  2:629648/139‑1 (MQ=255)
                                                                                       tCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                 >  2:645895/1‑103 (MQ=255)
                                                                                       tCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                 <  1:645895/103‑1 (MQ=255)
                                                                                       tCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg                                             >  2:578152/1‑139 (MQ=255)
                                                                                        cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcaca                                                                                                                           <  1:1357127/60‑1 (MQ=255)
                                                                                        cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcaca                                                                                                                           >  2:1357127/1‑60 (MQ=255)
                                                                                        cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         <  2:36273/62‑1 (MQ=255)
                                                                                        cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         >  1:36273/1‑62 (MQ=255)
                                                                                               cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTAGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGcc                                                                  <  1:943995/110‑1 (MQ=255)
                                                                                               cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCACGTCACCCAGATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGcc                                                                  >  2:943995/1‑110 (MQ=255)
                                                                                                 tacgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgg                                   <  2:146140/137‑6 (MQ=255)
                                                                                                  acgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                  <  2:480145/138‑7 (MQ=255)
                                                                                                       ggTCATAGTTGGTAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                                                       <  2:73037/81‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGTAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                                                       >  1:73037/1‑81 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   <  2:526439/37‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   >  1:526439/1‑37 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               <  2:534530/73‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               >  1:534530/1‑73 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                             >  1:440344/1‑128 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                             <  2:1215909/139‑12 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   >  1:334201/1‑53 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   >  1:274819/1‑53 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  2:334201/53‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  2:274819/53‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                       <  2:246225/113‑1 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                       >  1:246225/1‑113 (MQ=255)
                                                                                                       ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                             >  1:1409973/1‑139 (MQ=255)
                                                                                                           aTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  1:104260/49‑1 (MQ=255)
                                                                                                           aTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   >  2:104260/1‑49 (MQ=255)
                                                                                                           aTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                              >  1:1462227/1‑118 (MQ=255)
                                                                                                           aTAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                              <  2:1462227/118‑1 (MQ=255)
                                                                                                              gTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCg                      <  2:1409973/139‑1 (MQ=255)
                                                                                                                  gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGt                     >  1:1123724/1‑136 (MQ=255)
                                                                                                                  gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGt                     <  2:1123724/136‑1 (MQ=255)
                                                                                                                  gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                  <  1:272038/139‑1 (MQ=255)
                                                                                                                  gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                  <  1:1146824/139‑1 (MQ=255)
                                                                                                                    aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                  <  2:1016380/137‑1 (MQ=255)
                                                                                                                    aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                  >  1:1016380/1‑137 (MQ=255)
                                                                                                                             ggCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                <  2:371996/82‑1 (MQ=255)
                                                                                                                             ggCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                >  1:371996/1‑82 (MQ=255)
                                                                                                                              gCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGAt      <  2:1228714/139‑1 (MQ=255)
                                                                                                                                  cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGAag  >  1:1081910/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                          
TCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAG  >  NC_000913/566147‑566415

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑