Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2770CT56.4% 20.9 / 129.0 94pseudogene (100/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (28/25);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGC  >  NC_000913/566140‑566399
                                                                                                                                         |                                                                                                                           
tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                                                                                                            >  1:887337/1‑139 (MQ=255)
      gTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTcgc                                                                                                                      >  2:948758/1‑139 (MQ=255)
           ttttttCGATGAACTACGGCACGATCTCGTCAAAACTCGCCATGTAATTTTCGTCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAc                                                                                                                 >  1:274392/1‑139 (MQ=255)
             ttttCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                               <  1:849756/139‑1 (MQ=255)
                tCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACTACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcaca                                                                                                            <  2:887337/139‑1 (MQ=255)
                tCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCcaca                                                                                                            >  2:945549/1‑139 (MQ=255)
                  gATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGTCATTTTTTCGCGTCACCCACTTg                                                                                                          <  2:357222/139‑1 (MQ=255)
                   aTGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGc                                                                                                         <  1:406494/139‑1 (MQ=255)
                   aTGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGc                                                                                                         <  2:931622/139‑1 (MQ=255)
                             gCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑c                                                                                              <  1:496996/139‑1 (MQ=255)
                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑c                                                                                              <  2:361668/134‑1 (MQ=255)
                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑c                                                                                              >  1:361668/1‑134 (MQ=255)
                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTggg                                                                                         <  1:948758/139‑1 (MQ=255)
                                   tctcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcg                                                                                                                       >  2:485393/1‑109 (MQ=255)
                                   tctcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGc                                                                                        >  2:742858/1‑139 (MQ=255)
                                   tctcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCGTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcg                                                                                                                       <  1:485393/109‑1 (MQ=255)
                                       gTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                               <  1:253351/113‑1 (MQ=255)
                                       gTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                               >  2:253351/1‑113 (MQ=255)
                                        tCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                               >  2:1072806/1‑112 (MQ=255)
                                        tCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                                               <  1:1072806/112‑1 (MQ=255)
                                          aaaaCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTaa                                                                                 <  1:1043569/139‑1 (MQ=255)
                                           aaaCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAg                                                                                <  1:447358/139‑1 (MQ=255)
                                            aaCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGc                                                                               >  2:1102988/1‑139 (MQ=255)
                                             aCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGcc                                                                              <  2:25041/139‑1 (MQ=255)
                                             aCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGct                                                                              <  1:657295/139‑2 (MQ=255)
                                             aCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATAAGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGct                                                                              >  1:435473/1‑138 (MQ=255)
                                                 gCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGctgact                                                                          >  1:1032599/1‑134 (MQ=255)
                                                     tGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAt                                                                      <  2:1066239/139‑1 (MQ=255)
                                                        aCTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATAAGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtggc                                                                   <  2:435473/139‑4 (MQ=255)
                                                              cATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTc                                                                                                                        >  1:114702/1‑81 (MQ=255)
                                                              cATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTc                                                                                                                        <  2:114702/81‑1 (MQ=255)
                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                            >  2:1108714/4‑139 (MQ=255)
                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                            >  2:950162/4‑139 (MQ=255)
                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTgg                                                                                          >  2:32510/1‑107 (MQ=255)
                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTgg                                                                                          <  1:32510/107‑1 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTg                                                                                                  >  2:354303/2‑100 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTg                                                                                                  <  1:354303/99‑1 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTggg                                                                                         >  2:980191/2‑108 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTggg                                                                                         <  1:980191/107‑1 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacc                                                          >  2:38969/2‑139 (MQ=255)
                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                       <  2:19326/139‑1 (MQ=255)
                                                                      tCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAt                                                                      >  1:797016/1‑122 (MQ=255)
                                                                      tCAATCACGACATAATGCAGGCCTTCACCCTTCATACGCGGGTCATTGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGctgactctat                                                                      <  2:797016/122‑10 (MQ=255)
                                                                             cGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGaa                                              >  2:678664/1‑139 (MQ=255)
                                                                               aCATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAAtt                                            <  2:900455/139‑1 (MQ=255)
                                                                                cATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAActt                                           >  1:881887/1‑139 (MQ=255)
                                                                                     tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                             >  1:144000/1‑69 (MQ=255)
                                                                                     tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                             <  2:144000/69‑1 (MQ=255)
                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTc                                                                                                                   <  2:405457/58‑1 (MQ=255)
                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTc                                                                                                                   >  1:405457/1‑58 (MQ=255)
                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                        >  1:363169/1‑116 (MQ=255)
                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                        <  2:363169/116‑1 (MQ=255)
                                                                                          gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATcc                                 >  2:1094974/1‑139 (MQ=255)
                                                                                           ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATTTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                <  1:835964/139‑1 (MQ=255)
                                                                                           ccTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                <  1:963424/139‑1 (MQ=255)
                                                                                              tCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                             <  1:945549/139‑1 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGc                                                                   <  1:51055/100‑1 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGc                                                                   >  2:51055/1‑100 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                            >  1:712908/1‑139 (MQ=255)
                                                                                                    ttCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                       >  2:329498/1‑138 (MQ=255)
                                                                                                         acgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                  <  2:1032599/138‑7 (MQ=255)
                                                                                                            cgGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTAtt                >  1:112457/1‑139 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt             <  1:38969/139‑12 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt             >  2:1210932/1‑128 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt             >  2:423621/1‑128 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt             <  2:799207/139‑12 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgg                   <  1:1024223/133‑6 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgg                   >  2:1024223/1‑128 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                              >  2:1205831/1‑122 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                              <  1:1205831/122‑1 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc             >  1:419293/1‑139 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc             >  1:62567/1‑139 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc             >  2:818210/1‑139 (MQ=255)
                                                                                                                tCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcatttca           >  1:180466/1‑126 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCGTCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTACCTCGAAACCACCGAGCCGGAATTTCGTGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  >  1:751061/1‑139 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCTCCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtg                                                                     <  2:877931/72‑2 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑c                                                                                              >  2:334720/1‑47 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑c                                                                                              <  1:334720/47‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTaa                                                                                 <  1:952487/60‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTaa                                                                                 >  2:952487/1‑60 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTa                                                                       <  2:921310/70‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTa                                                                       >  1:921310/1‑70 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTAtg                                                                     >  1:877931/1‑71 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                       >  1:881257/1‑102 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                       <  2:881257/102‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                             >  2:106105/1‑112 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                             <  1:106105/112‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                             <  2:850844/112‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCggg                             >  1:850844/1‑112 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  >  1:1191139/1‑139 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  >  1:886009/1‑139 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTGCGCGTCACCCACATGCTGTACTGCA‑CCTGTGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAGCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  <  2:608703/139‑1 (MQ=255)
                                                                                                                              gTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                       <  2:1183496/97‑1 (MQ=255)
                                                                                                                              gTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                       >  1:1183496/1‑97 (MQ=255)
                                                                                                                                         cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCg            >  1:316408/1‑113 (MQ=255)
                                                                                                                                         cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCg            <  2:316408/113‑1 (MQ=255)
                                                                                                                                         cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  >  2:38175/1‑123 (MQ=255)
                                                                                                                                         cTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc  <  1:38175/123‑1 (MQ=255)
                                                                                                                                         |                                                                                                                           
TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCA‑CCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGC  >  NC_000913/566140‑566399

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