Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2770CT55.2% 21.2 / 189.6 143pseudogene (100/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (39/40);  total (69/74)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 144 total aligned reads displayed.

TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAG  >  NC_000913/566140‑566415
                                                                                                                                         |                                                                                                                                          
tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                                                                                                                           >  1:1439889/1‑139 (MQ=255)
tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCt                                                                                                                                           >  2:1283560/1‑139 (MQ=255)
      gTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTcgc                                                                                                                                     >  1:982961/1‑139 (MQ=255)
      gTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTcgc                                                                                                                                     <  1:762829/139‑1 (MQ=255)
       tCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcgcg                                                                                                                                    >  2:1224702/1‑139 (MQ=255)
        ccATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATACTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGt                                                                                                                                   <  2:846886/139‑1 (MQ=255)
                  gATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                                                                         <  2:1609310/139‑1 (MQ=255)
                        cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  2:1521619/139‑1 (MQ=255)
                         tCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACt                                                                                                                  <  1:1674721/139‑1 (MQ=255)
                         tCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACt                                                                                                                  <  2:1593016/139‑1 (MQ=255)
                             gCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCAc                                                                                                              <  2:1079270/139‑1 (MQ=255)
                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          <  1:1456451/138‑1 (MQ=255)
                                  aTCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          >  2:1456451/1‑138 (MQ=255)
                                   tctcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGc                                                                                                        >  1:1059340/1‑139 (MQ=255)
                                   tctcGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGc                                                                                                        >  1:1509708/1‑139 (MQ=255)
                                     tcGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCa                                                                                                      >  1:225699/1‑139 (MQ=255)
                                            aaCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               <  2:1439889/139‑1 (MQ=255)
                                            aaCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               <  2:841350/139‑1 (MQ=255)
                                            aaCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                                                               <  2:1509708/139‑1 (MQ=255)
                                              cTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGctg                                                                                             <  1:1246942/139‑3 (MQ=255)
                                                 gCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             >  2:119090/1‑120 (MQ=255)
                                                 gCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             <  1:119090/120‑1 (MQ=255)
                                                    aTGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                                                       <  2:868811/139‑1 (MQ=255)
                                                     tGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTAt                                                                                      <  1:806458/139‑1 (MQ=255)
                                                      gTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTAtt                                                                                     <  1:459164/139‑1 (MQ=255)
                                                        aCTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACGAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          >  2:982359/1‑116 (MQ=255)
                                                        aCTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtggc                                                                                   <  1:651386/139‑4 (MQ=255)
                                                         cTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          >  1:686473/1‑115 (MQ=255)
                                                          ttttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                 <  1:1283560/139‑1 (MQ=255)
                                                            ttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             <  2:909344/109‑1 (MQ=255)
                                                            ttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             >  1:909344/1‑109 (MQ=255)
                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACTAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGaaa                                                                            >  2:986722/4‑139 (MQ=255)
                                                               atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGc                                                                                                        >  2:369927/4‑111 (MQ=255)
                                                                 tCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAacc                                                                          >  1:1472712/1‑139 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCcac                                                                                                                            <  1:269052/88‑1 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCCCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAacca                                                                         >  2:1200981/2‑139 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   >  2:1375234/2‑97 (MQ=255)
                                                                 cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  1:1375234/96‑1 (MQ=255)
                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTAt                                                                                      >  2:1419316/1‑124 (MQ=255)
                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTAt                                                                                      <  1:1419316/124‑1 (MQ=255)
                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCAt                                                                                                     <  2:24916/109‑1 (MQ=255)
                                                                    gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                       >  1:39039/1‑139 (MQ=255)
                                                                     cTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGACAAAGTACCAGGCATCTTTTCGCGt                                                                                                                                   >  2:1445937/1‑78 (MQ=255)
                                                                     cTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGACAAAGTACCAGGCATCTTTTCGCGt                                                                                                                                   <  1:1445937/78‑1 (MQ=255)
                                                                        aaTCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGc                                                                   <  1:960279/139‑1 (MQ=255)
                                                                                    aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                 >  1:313884/1‑113 (MQ=255)
                                                                                    aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCt                                                                                 <  2:313884/113‑1 (MQ=255)
                                                                                          gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATcc                                                 >  2:1496169/1‑139 (MQ=255)
                                                                                            cTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCg                                               <  1:234180/139‑1 (MQ=255)
                                                                                            cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCg                                               >  2:530175/1‑139 (MQ=255)
                                                                                             ttCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCgg                                              <  2:1472712/139‑1 (MQ=255)
                                                                                               cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa                                            >  2:275293/1‑139 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                                                                   <  2:209831/52‑1 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   <  1:1319274/68‑1 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                                                                   >  2:1319274/1‑68 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCAcc                                                                                                             >  1:317245/1‑74 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaact                                                            <  2:856156/123‑3 (MQ=255)
                                                                                               cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaact                                                            >  1:856156/1‑121 (MQ=255)
                                                                                                  ccTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGCGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGt                                         <  2:995108/138‑1 (MQ=255)
                                                                                                      catacgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                            >  1:683536/5‑132 (MQ=255)
                                                                                                      cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCttt                                                                                                                                         <  2:81916/39‑1 (MQ=255)
                                                                                                      cATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                                                                        >  1:1328573/1‑104 (MQ=255)
                                                                                                      cATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                            <  2:683536/132‑1 (MQ=255)
                                                                                                             gggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGc                              <  1:1496169/139‑1 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              >  2:699885/1‑42 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                                                              <  1:699885/42‑1 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGct                                                                                              <  1:706364/74‑2 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaac                                                             <  2:81650/107‑2 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGaac                                                             >  1:81650/1‑106 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                <  1:1204500/120‑1 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                >  2:1204500/1‑120 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTggg                                                                                                         >  2:829651/1‑63 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTggg                                                                                                         <  1:829651/63‑1 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                             >  1:339718/1‑139 (MQ=255)
                                                                                                              ggTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGcc                             >  1:256104/1‑139 (MQ=255)
                                                                                                                  aTAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAg                                                                                                >  2:204770/1‑68 (MQ=255)
                                                                                                                  aTAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAg                                                                                                <  1:204770/68‑1 (MQ=255)
                                                                                                                  aTAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                >  1:722163/1‑116 (MQ=255)
                                                                                                                  aTAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATccc                                                <  2:722163/116‑1 (MQ=255)
                                                                                                                     gTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCg                      <  2:339718/139‑1 (MQ=255)
                                                                                                                      ttGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGt                     <  2:256104/139‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACTCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                <  1:716349/93‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGt                                                                                                   >  1:158/1‑58 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAg                                                                                                >  2:1344085/1‑61 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGa                                                                              <  2:302871/79‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCgg                                                                >  2:716349/1‑93 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCATATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                       <  1:829255/102‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACt                                                                                                                  <  2:67066/43‑1 (MQ=255)
                                                                                                                         gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGcc                                                                  >  2:1422650/1‑91 (MQ=255)
                                                                                                                           aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                       <  2:664633/68‑1 (MQ=255)
                                                                                                                           aaaGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTa                                                                                       >  1:664633/1‑68 (MQ=255)
                                                                                                                             aGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                                                       >  2:132722/1‑66 (MQ=255)
                                                                                                                             aGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTa                                                                                       <  1:132722/66‑1 (MQ=255)
                                                                                                                              gTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                                                                   <  2:55968/69‑1 (MQ=255)
                                                                                                                               tACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          >  2:1040029/1‑45 (MQ=255)
                                                                                                                               tACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTgg                                                                                                          <  1:1040029/45‑1 (MQ=255)
                                                                                                                                aCCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACTGAGCTGGAACTTCATGAAATCCCGGGAGGTaa                                       <  2:1006739/111‑2 (MQ=255)
                                                                                                                                aCCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTGATGAAATCCCGGGAGGTaa                                       >  1:1006739/1‑110 (MQ=255)
                                                                                                                                 ccAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                  <  2:703929/115‑7 (MQ=255)
                                                                                                                                         cTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCag  >  2:1403120/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAG  >  NC_000913/566140‑566415

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