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

Only 100 of 109 total aligned reads displayed.

TCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGA  >  NC_000913/566111‑566373
                                                                                                                                      |                                                                                                                                
tCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGcgcg                                                                                                                              >  1:1339898/1‑139 (MQ=255)
   ccAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGt                                                                                                                           <  1:622955/139‑1 (MQ=255)
         aCAAAACCAATTTCCGCCAGTGCCTCGTCCATTGTTTCGATGAACGCCGGCACCATCTCGTCAAAACGCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGt                                                                                                                     <  1:1184492/139‑1 (MQ=255)
            aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgg                                                                                                                  >  2:502675/1‑139 (MQ=255)
            aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgg                                                                                                                  >  2:1018636/1‑139 (MQ=255)
            aaaCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTgg                                                                                                                  <  1:693423/139‑1 (MQ=255)
              aCCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCa                                                                                                                <  2:1339898/139‑1 (MQ=255)
                  aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCACTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCa                                                                                                                <  1:1359128/135‑1 (MQ=255)
                  aTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCACTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCa                                                                                                                >  2:1359128/1‑135 (MQ=255)
                        gCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAgg                                                                                                      >  2:618441/1‑139 (MQ=255)
                         ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCATGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                >  1:1132200/1‑112 (MQ=255)
                         ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCggg                                                                                                                            <  2:1292128/116‑1 (MQ=255)
                         ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCggg                                                                                                                            >  1:1292128/1‑116 (MQ=255)
                         ccAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGc                                                                                                     >  1:493666/1‑139 (MQ=255)
                         ccAGTGCCTCGTCAATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACACGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                <  2:1132200/112‑1 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacgc                                                                                                                               >  1:659049/1‑109 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacgc                                                                                                                               <  2:659049/109‑1 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCAttt                                                                                                 >  1:334362/1‑137 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCa                                                                                                    <  1:1207140/136‑1 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCa                                                                                                    <  1:1329157/136‑1 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCa                                                                                                    >  2:1207140/1‑136 (MQ=255)
                             tGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCa                                                                                                    >  2:1329157/1‑136 (MQ=255)
                                  cGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcg                                                                                            <  1:794543/139‑1 (MQ=255)
                                   gTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                >  1:835423/1‑102 (MQ=255)
                                   gTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                <  2:835423/102‑1 (MQ=255)
                                     ccATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                         <  1:1018636/139‑1 (MQ=255)
                                       aTTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCaaa                                                                                                              >  1:1502050/1‑116 (MQ=255)
                                       aTTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCaaa                                                                                                              <  2:1502050/116‑1 (MQ=255)
                                               gATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               <  2:836702/139‑1 (MQ=255)
                                               gATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               <  2:707349/139‑1 (MQ=255)
                                               gATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCAGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               <  1:884310/139‑1 (MQ=255)
                                                 tGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCt                                                                             >  2:885764/1‑139 (MQ=255)
                                                     cTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                         <  2:334362/139‑1 (MQ=255)
                                                     cTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTAc                                                                         <  2:493666/139‑1 (MQ=255)
                                                               aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                                                                                                >  2:591910/1‑74 (MQ=255)
                                                               aTCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                                                                                                <  1:591910/74‑1 (MQ=255)
                                                                     tCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTa                                                                                                           <  1:1178027/89‑1 (MQ=255)
                                                                     tCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTa                                                                                                           >  2:1178027/1‑89 (MQ=255)
                                                                         aaCTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                     <  2:660129/139‑1 (MQ=255)
                                                                          aCTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGct                                                    <  1:502675/139‑2 (MQ=255)
                                                                                        tttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATgcg                                                                                                                                <  2:417750/49‑1 (MQ=255)
                                                                                         ttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                    <  1:1201290/92‑1 (MQ=255)
                                                                                         ttCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCAccc                                                                                    >  2:1201290/1‑92 (MQ=255)
                                                                                         ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGtt                                                                                                                    <  2:856141/60‑1 (MQ=255)
                                                                                         ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGtt                                                                                                                    >  1:856141/1‑60 (MQ=255)
                                                                                            atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                >  1:316572/4‑45 (MQ=38)
                                                                                            atcccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                >  1:526573/4‑45 (MQ=38)
                                                                                             ttCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               >  2:1485704/1‑93 (MQ=255)
                                                                                              cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTc                                                                                                                          <  1:1136251/48‑1 (MQ=39)
                                                                                              cccGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTc                                                                                                                          >  2:1136251/2‑49 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGTCCTTCACGCTTCATCCGCGGGTCATAGTTGGCAAAGTACCCGGCATTTTTTCGCGTCACCCACACGCTCTACTGCACCTGGGGCATGTAAGCTGACTTTATGTCCTCGCAaccacc                             >  1:35274/1‑139 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATacg                                                                                                                                >  1:431284/1‑40 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                     <  1:1410629/115‑1 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGc                                                     >  2:1410629/1‑115 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                              <  1:1152243/138‑1 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccac                              >  2:1152243/1‑138 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                    >  1:258399/1‑84 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCAccc                                                                                    <  2:258399/84‑1 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                             >  1:1344016/1‑139 (MQ=255)
                                                                                                 gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                             <  1:1043731/139‑1 (MQ=255)
                                                                                                      aTCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGcc                        >  1:1125613/1‑139 (MQ=255)
                                                                                                           gACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                         <  2:723700/69‑1 (MQ=255)
                                                                                                           gACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGt                                                                                         >  1:723700/1‑69 (MQ=255)
                                                                                                            aCATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                         <  2:1048297/116‑1 (MQ=255)
                                                                                                            aCATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGc                                         >  1:1048297/1‑116 (MQ=255)
                                                                                                                 aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGt                                                                                         <  1:332389/63‑1 (MQ=255)
                                                                                                                 aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGt                                                                                         >  2:332389/1‑63 (MQ=255)
                                                                                                                 aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtg                                           >  2:1062602/1‑108 (MQ=255)
                                                                                                                 aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg             <  1:618441/139‑1 (MQ=255)
                                                                                                                 aTGCAGGCCGTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtg                                           <  1:1062602/109‑2 (MQ=255)
                                                                                                                    cAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               <  1:1533711/70‑1 (MQ=255)
                                                                                                                     aGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGc                                                                                                     >  2:797051/1‑47 (MQ=255)
                                                                                                                     aGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGc                                                                                                     <  1:797051/47‑1 (MQ=255)
                                                                                                                     aGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAt         >  2:1576315/1‑139 (MQ=255)
                                                                                                                     aGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAAt         >  2:20156/1‑139 (MQ=255)
                                                                                                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               >  2:1266940/1‑67 (MQ=255)
                                                                                                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATg                                                                               <  1:1266940/67‑1 (MQ=255)
                                                                                                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATg                                                          >  2:812533/1‑88 (MQ=255)
                                                                                                                       gCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATg                                                          <  1:812533/88‑1 (MQ=255)
                                                                                                                       gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATcc       >  1:423749/1‑139 (MQ=255)
                                                                                                                        ccTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATccc      >  2:1486924/1‑139 (MQ=255)
                                                                                                                         cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCt                                                                  <  1:1471619/78‑1 (MQ=255)
                                                                                                                         cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCt                                                                  <  1:404941/78‑1 (MQ=255)
                                                                                                                         cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCt                                                                  >  2:404941/1‑78 (MQ=255)
                                                                                                                         cTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCt                                                                  >  2:1471619/1‑78 (MQ=255)
                                                                                                                           tCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTACGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCggg   >  1:591495/1‑139 (MQ=255)
                                                                                                                            cACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGa  >  2:1130000/1‑139 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                         <  1:1003968/68‑1 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTAc                                                                         >  2:1003968/1‑68 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa  >  2:1521835/1‑139 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa  >  1:274815/1‑139 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa  >  1:1036020/1‑139 (MQ=255)
                                                                                                                            cACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa  >  1:402318/1‑139 (MQ=255)
                                                                                                                                 ttCATGCGCGGGTCATAGTTGGCAAAGTTCCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGa                     >  1:937651/1‑115 (MQ=255)
                                                                                                                                 ttCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGa                     <  2:1391468/115‑1 (MQ=255)
                                                                                                                                 ttCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGa                     >  1:1391468/1‑115 (MQ=255)
                                                                                                                                     tGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGt                                                                                         >  1:435337/1‑43 (MQ=255)
                                                                                                                                     tGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGt                                                                                         <  2:435337/43‑1 (MQ=255)
                                                                                                                                     tGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCt                                                                             >  2:1486595/1‑55 (MQ=255)
                                                                                                                                     tGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCt                                                                             <  1:1486595/55‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
TCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGA  >  NC_000913/566111‑566373

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