Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,277	0	C	T	53.4%	‑2.6 / 120.6	88	pseudogene (100/279 nt)	exoD	pseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base T (21/26);  total (39/49)
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.67e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGA  >  NC_000913/566140‑566416                                                                                                                                           |                                                                                                                                            tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                                                            >  2:572078/1‑139 (MQ=255) tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                                                            >  1:571079/1‑139 (MQ=255) tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCt                                                                                                                                            >  1:1251574/1‑139 (MQ=255)           aTTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCa                                                                                                                                  <  1:1137215/139‑1 (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)                         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)                                                     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‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                                                                                                       <  2:611277/95‑1 (MQ=255)                                                                  tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTg                                                                                                                       >  1:611277/1‑95 (MQ=255)                                                                  cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                               <  2:144817/86‑1 (MQ=255)                                                                  cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                               >  1:144817/2‑87 (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)                                                                           tCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCg                                                                  <  2:442201/139‑1 (MQ=255)                                                                              cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGaa                                                               >  1:534507/1‑139 (MQ=255)                                                                                   tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                    <  1:510036/65‑1 (MQ=255)                                                                                   tAATGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                    >  2:510036/1‑65 (MQ=255)                                                                                    aaTGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCAt                                                         <  1:302396/139‑1 (MQ=255)                                                                                     aTGCAGGCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                        >  1:167055/1‑139 (MQ=255)                                                                                         aGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                          <  1:425626/69‑1 (MQ=255)                                                                                         aGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                          >  2:425626/1‑69 (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                                                                                                                                    <  1:747300/57‑1 (MQ=255)                                                                                           gCCTTCACGCTTCATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                    >  2:747300/1‑57 (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                                                                                                                               >  2:1330453/1‑59 (MQ=255)                                                                                              ttCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                                               <  1:1330453/59‑1 (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                                                                                                                          >  2:621403/5‑55 (MQ=255)                                                                                                       cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATg                                                                                                                          <  1:621403/55‑1 (MQ=255)                                                                                                       cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                    <  2:998489/61‑1 (MQ=255)                                                                                                       cATACGCGGGTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                    >  1:998489/1‑61 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                               <  1:473571/42‑1 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                               >  2:473571/1‑42 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTg                                                                                                                       >  2:784643/1‑50 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTg                                                                                                                       <  1:784643/50‑1 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCa                                                                                                       >  1:1140209/1‑66 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCa                                                                                                       <  2:1140209/66‑1 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                              >  1:434455/1‑128 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                              >  1:24732/1‑128 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                              <  2:434455/139‑12 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                              <  1:1213754/59‑1 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                              >  2:1213754/1‑59 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTAt                                  >  1:1360101/1‑135 (MQ=255)                                                                                                               ggTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTAt                                  <  2:1360101/135‑1 (MQ=255)                                                                                                                 tCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                          <  2:569254/125‑1 (MQ=255)                                                                                                                 tCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                          >  1:569254/1‑125 (MQ=255)                                                                                                                 tCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCgg                            >  2:926919/1‑139 (MQ=255)                                                                                                                    tAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcaca                                                                                                                            >  2:1229199/1‑40 (MQ=255)                                                                                                                    tAG‑TTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcaca                                                                                                                            <  1:1229199/40‑1 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                              <  2:442331/112‑1 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                                              >  1:442331/1‑112 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                             <  2:1101327/113‑1 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                             >  1:1101327/1‑113 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                   <  1:61775/139‑1 (MQ=255)                                                                                                                           gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                   >  2:61775/1‑139 (MQ=255)                                                                                                                                           cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCag   >  2:587043/1‑139 (MQ=255)                                                                                                                                           cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGA‑TTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCaga  >  2:1031154/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                            TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAG‑TTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGA  >  NC_000913/566140‑566416

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