Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2050TC50.7% 5.8 / 93.0 73pseudogene (172/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base C (20/17);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGA  >  NC_000913/566072‑566339
                                                                                                                                     |                                                                                                                                       
tACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                                                    <  2:979837/139‑1 (MQ=255)
tACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCt                                                                                                                                    >  2:663073/1‑139 (MQ=255)
     ggATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATC‑CGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              <  1:955503/139‑1 (MQ=255)
      gATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              <  1:332788/139‑1 (MQ=255)
           tatCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         <  1:408604/139‑1 (MQ=255)
           tatCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         >  1:21547/1‑139 (MQ=255)
                tGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATg                                                                                                                    >  1:442201/1‑139 (MQ=255)
                tGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATg                                                                                                                    >  1:760679/1‑139 (MQ=255)
                 gAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGc                                                                                                                   >  2:585341/1‑139 (MQ=255)
                  aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTc                                                                                                                                   >  1:946126/1‑122 (MQ=255)
                  aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTc                                                                                                                                   <  2:946126/122‑1 (MQ=255)
                  aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCa                                                                                                                  >  2:768458/1‑139 (MQ=255)
                  aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                                                                     >  2:1163962/1‑120 (MQ=255)
                  aGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGc                                                                                                                                     <  1:1163962/120‑1 (MQ=255)
                    gATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAgg                                                                                                                >  1:924006/1‑139 (MQ=255)
                    gATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAgg                                                                                                                >  2:1334331/1‑139 (MQ=255)
                      tGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         >  2:770546/1‑128 (MQ=255)
                      tGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         <  1:770546/128‑1 (MQ=255)
                       gCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCt                                                                                                             >  2:918114/1‑139 (MQ=255)
                         gTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGGTGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGAAGGCCTTc                                                                                                           >  2:812485/1‑139 (MQ=255)
                         gTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTc                                                                                                           >  2:1137215/1‑139 (MQ=255)
                              cGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCt                                                                                                      >  2:302396/1‑139 (MQ=255)
                                    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                                                                                >  1:52273/1‑139 (MQ=255)
                                                   aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTgg                                                                                >  2:923002/1‑139 (MQ=255)
                                                     aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                              <  1:1059703/139‑1 (MQ=255)
                                                         aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              >  1:1055194/1‑88 (MQ=255)
                                                         aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              <  2:1055194/88‑1 (MQ=255)
                                                          tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                                                    >  2:964853/1‑81 (MQ=255)
                                                          tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCt                                                                                                                                    <  1:964853/81‑1 (MQ=255)
                                                          tttCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTa                                                                         >  2:428705/1‑139 (MQ=255)
                                                               gCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              <  2:1322154/82‑1 (MQ=255)
                                                               gCCAGTCCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCa                                                                                                                              >  1:1322154/1‑82 (MQ=255)
                                                                 cAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         >  1:1187352/1‑85 (MQ=255)
                                                                 cAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACa                                                                                                                         <  2:1187352/85‑1 (MQ=255)
                                                                    tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                               >  1:1251574/1‑139 (MQ=255)
                                                                    tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                               >  1:571079/1‑139 (MQ=255)
                                                                    tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                               >  2:572078/1‑139 (MQ=255)
                                                                              aTTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCa                                                     <  1:1137215/139‑1 (MQ=255)
                                                                                tttttCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCAtacg                                                                                               <  1:1206827/96‑4 (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                                                                                           <  1:503499/90‑1 (MQ=255)
                                                                                          aaCTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCggg                                                                                           >  2:503499/1‑90 (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)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAAt                                                                                                                                <  2:416314/40‑1 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAAt                                                                                                                                >  1:416314/1‑40 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACa                                                                                                                         <  2:1240883/47‑1 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACa                                                                                                                         >  1:1240883/1‑47 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCa                                                                                                                  <  1:1127822/54‑1 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCa                                                                                                                  >  2:1127822/1‑54 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGcc                                                                                                              <  1:903716/58‑1 (MQ=255)
                                                                                                       tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGcc                                                                                                              >  2:903716/1‑58 (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)
                                                                                                                                     tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                              >  1:133955/1‑59 (MQ=255)
                                                                                                                                     tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCa                                                                              <  2:133955/59‑1 (MQ=255)
                                                                                                                                     tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                          <  2:611277/95‑1 (MQ=255)
                                                                                                                                     tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                          >  1:611277/1‑95 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGA  >  NC_000913/566072‑566339

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