Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2450GA57.7% 19.9 / 93.2 79pseudogene (132/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (21/24);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTAC  >  NC_000913/566108‑566378
                                                                                                                                         |                                                                                                                                      
tGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgc                                                                                                                                       >  1:322762/1‑139 (MQ=255)
 gCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                                                                                                      <  2:760679/139‑1 (MQ=255)
      ccAAATACAAAACCAATTTAAGCCCGTGCCCCGCCCATTTTTTCGCTGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGt                                                                                                                                 >  1:336448/1‑139 (MQ=255)
               aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTgg                                                                                                                       >  2:63786/1‑139 (MQ=255)
               aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTgg                                                                                                                       >  2:923002/1‑139 (MQ=255)
               aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTgg                                                                                                                       >  1:52273/1‑139 (MQ=255)
                 aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                                                                     <  1:1059703/139‑1 (MQ=255)
                      tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTa                                                                                                                >  2:428705/1‑139 (MQ=255)
                                tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                      >  2:572078/1‑139 (MQ=255)
                                tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                      >  1:1251574/1‑139 (MQ=255)
                                tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                      >  1:571079/1‑139 (MQ=255)
                                          aTTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCa                                                                                            <  1:1137215/139‑1 (MQ=255)
                                            tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                      >  2:1206827/1‑96 (MQ=255)
                                             ttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                         <  2:322762/139‑1 (MQ=255)
                                                  gATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATg                                                                                    <  2:336448/139‑1 (MQ=255)
                                                   aTGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGc                                                                                   <  1:663073/139‑1 (MQ=255)
                                                    tGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCt                                                                                  <  1:614298/139‑1 (MQ=255)
                                                      aaCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCggg                                                                                                                                  >  2:503499/1‑90 (MQ=255)
                                                      aaCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCggg                                                                                                                                  <  1:503499/90‑1 (MQ=255)
                                                        cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                              <  1:1319049/139‑1 (MQ=255)
                                                        cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                              <  1:441724/139‑1 (MQ=255)
                                                             gCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAc                                                                         <  1:923002/139‑1 (MQ=255)
                                                              cACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                        <  2:1251574/139‑1 (MQ=255)
                                                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGTTCATAG‑TTGGCAAAGTACCAGGCATCTTCTCGCGTCACCCACATGCTGTACTGCACCTggg                                                                    <  1:812485/139‑1 (MQ=255)
                                                                     tcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCa                                                                 >  1:475428/1‑139 (MQ=255)
                                                                     tcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCa                                                                 <  1:428705/139‑1 (MQ=255)
                                                                       gTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATg                                                               >  2:565613/1‑139 (MQ=255)
                                                                       gTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTg                                                                      >  1:1093816/1‑132 (MQ=255)
                                                                       gTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTg                                                                      <  2:1093816/132‑1 (MQ=255)
                                                                        tCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGt                                                                                                                                 >  1:451561/1‑73 (MQ=255)
                                                                        tCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGt                                                                                                                                 <  2:451561/73‑1 (MQ=255)
                                                                                    aTGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                  <  1:918114/139‑1 (MQ=255)
                                                                                            ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGACATGTAAGCTGACTTTATGGCCTc                                           >  2:763227/1‑139 (MQ=255)
                                                                                             tCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGa                                         <  2:21547/139‑1 (MQ=255)
                                                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                         <  2:1077168/86‑1 (MQ=255)
                                                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                         >  1:1077168/4‑89 (MQ=255)
                                                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                                                                     <  2:133955/59‑1 (MQ=255)
                                                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                                                                     >  1:133955/1‑59 (MQ=255)
                                                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                                                                 >  1:611277/1‑95 (MQ=255)
                                                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                                                                 <  2:611277/95‑1 (MQ=255)
                                                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCa                                                                                                                     <  2:378810/58‑1 (MQ=255)
                                                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCa                                                                                                                     >  1:378810/2‑59 (MQ=255)
                                                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                         >  1:144817/2‑87 (MQ=255)
                                                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                         <  2:144817/86‑1 (MQ=255)
                                                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                    >  1:399541/1‑89 (MQ=255)
                                                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                    <  2:399541/89‑1 (MQ=255)
                                                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                  >  2:295398/1‑139 (MQ=255)
                                                                                                     cTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCg                                 <  1:743247/139‑1 (MQ=255)
                                                                                                     cTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TCGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCg                                 <  1:1334331/139‑1 (MQ=255)
                                                                                                      tCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAg                                                                                                                            >  2:277027/1‑48 (MQ=255)
                                                                                                      tCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAg                                                                                                                            <  1:277027/48‑1 (MQ=255)
                                                                                                          tCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCg                            <  2:442201/139‑1 (MQ=255)
                                                                                                             cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGaa                         >  1:534507/1‑139 (MQ=255)
                                                                                                                  tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                              >  2:510036/1‑65 (MQ=255)
                                                                                                                  tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                              <  1:510036/65‑1 (MQ=255)
                                                                                                                   aaTGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCAt                   <  1:302396/139‑1 (MQ=255)
                                                                                                                    aTGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                  >  1:167055/1‑139 (MQ=255)
                                                                                                                        aGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                    >  2:425626/1‑69 (MQ=255)
                                                                                                                        aGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                    <  1:425626/69‑1 (MQ=255)
                                                                                                                        aGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcaca                                                                                      >  1:299441/1‑67 (MQ=255)
                                                                                                                        aGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcaca                                                                                      <  2:299441/67‑1 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcac                                                                                       >  2:754993/1‑64 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcac                                                                                       <  1:754993/64‑1 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                              >  2:747300/1‑57 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                              <  1:747300/57‑1 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtg                                                >  2:167155/1‑102 (MQ=255)
                                                                                                                          gCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtg                                                <  1:167155/103‑2 (MQ=255)
                                                                                                                            cTTCACGCTTAATCCTCGGTTCATAG‑TTGGCAACGTCCCAGGCATTTTTCCGCGTCACCACCATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg          <  1:295398/139‑1 (MQ=255)
                                                                                                                             ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                         <  1:1330453/59‑1 (MQ=255)
                                                                                                                             ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                         >  2:1330453/1‑59 (MQ=255)
                                                                                                                             ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg         <  2:44760/139‑1 (MQ=255)
                                                                                                                             ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg         <  2:534507/139‑1 (MQ=255)
                                                                                                                             ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg         <  2:571079/139‑1 (MQ=255)
                                                                                                                             ttCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCgg         >  1:810665/1‑139 (MQ=255)
                                                                                                                               cACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAt              >  2:451530/1‑132 (MQ=255)
                                                                                                                               cACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAt              <  1:451530/132‑1 (MQ=255)
                                                                                                                               cACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa       >  1:1237506/1‑139 (MQ=255)
                                                                                                                                   cTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTa   <  1:572078/139‑1 (MQ=255)
                                                                                                                                    ttCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTAc  <  1:565613/139‑1 (MQ=255)
                                                                                                                                      cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATg                                                                                    <  1:621403/55‑1 (MQ=255)
                                                                                                                                      cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                              >  1:998489/1‑61 (MQ=255)
                                                                                                                                      cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                              <  2:998489/61‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTAC  >  NC_000913/566108‑566378

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: ‑