Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7470GA44.5% 328.0 / inf 512pseudogene (171/309 nt)aaaDDLP12 prophage; predicted tail fiber assembly protein (pseudogene);Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (156/128);  new base A (120/108);  total (276/236)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 538 total aligned reads displayed.

GGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATG  >  NC_000913/580647‑580847
                                                                                                    |                                                                                                    
 gCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGc                                                                                                     <  1:1309877/101‑1 (MQ=255)
  ccGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCg                                                                                                    >  2:5397171/1‑101 (MQ=255)
  caggAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCACCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCg                                                                                                    <  1:5774422/99‑1 (MQ=255)
   ctgAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACgc                                                                                                   >  1:5950674/3‑100 (MQ=255)
            aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATg                                                                                          >  1:4633845/1‑101 (MQ=255)
               cggcggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAAt                                                                                       <  2:3932798/96‑1 (MQ=255)
                ggcggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGCGCGCTGGGAAACGTTGCCGTTTCAACGCGCGCTCATGAATg                                                                                      >  1:780991/5‑101 (MQ=255)
                 gcggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGc                                                                                     <  1:5283116/98‑1 (MQ=255)
                  ccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCg                                                                                    >  1:1513637/1‑101 (MQ=255)
                   ggTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGa                                                                                   >  1:1911673/2‑101 (MQ=255)
                   cgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGa                                                                                   >  2:5920211/1‑101 (MQ=255)
                   cgTTGAATGAGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGa                                                                                   <  1:887219/101‑1 (MQ=255)
                     ttGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCTTTTCAGCGGGCCATCATGAATGCGATg                                                                                 <  1:4101881/101‑1 (MQ=255)
                        taTGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATgggt                                                                              <  2:5756820/100‑2 (MQ=255)
                         aTGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCa                                                                             <  1:2900774/101‑1 (MQ=255)
                          tGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCA‑GAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGc                                                                           <  2:2338168/101‑1 (MQ=255)
                          tGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATgggtaa                                                                            >  1:5486376/1‑98 (MQ=255)
                            ggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCGGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCg                                                                          <  1:7062349/101‑1 (MQ=255)
                               ggATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACt                                                                       >  1:6545942/1‑101 (MQ=255)
                               ggATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACt                                                                       >  1:1044868/1‑101 (MQ=255)
                               ggATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACt                                                                       >  1:3278696/1‑101 (MQ=255)
                                 acgatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTAc                                                                     >  2:2168079/5‑101 (MQ=255)
                                 acgatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTAc                                                                     <  1:276377/97‑1 (MQ=255)
                                 aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTAc                                                                     <  2:602301/101‑1 (MQ=255)
                                   gatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                   >  1:3486646/3‑101 (MQ=255)
                                    cTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATc                                                                  >  2:3799504/1‑101 (MQ=255)
                                    atAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGACGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                  >  1:1643073/2‑100 (MQ=255)
                                     tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                 <  1:2379500/101‑1 (MQ=255)
                                     tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                 >  1:2090375/1‑101 (MQ=255)
                                     tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                 <  1:4220510/101‑1 (MQ=255)
                                     tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                 >  2:7005204/1‑101 (MQ=255)
                                      aaTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCg                                                                >  2:884635/1‑101 (MQ=255)
                                      aaTTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACg                                                                >  1:6656554/1‑101 (MQ=255)
                                       aTTACTATCTCCCGAAAGAATCCGCCTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGt                                                               <  2:1945448/101‑1 (MQ=255)
                                          aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAg                                                            >  1:5602928/1‑101 (MQ=255)
                                           cTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAgg                                                           <  2:2442281/101‑1 (MQ=255)
                                           cTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAgg                                                           <  1:2530141/101‑1 (MQ=255)
                                            tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGt                                                          >  1:5023120/1‑101 (MQ=255)
                                            tATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGGGCTGGGAAACGGTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGt                                                          >  2:6833406/1‑101 (MQ=255)
                                            tATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGt                                                          <  2:4604573/101‑1 (MQ=255)
                                             aTCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTg                                                         >  2:4953979/1‑101 (MQ=255)
                                              tctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTCTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                        <  1:6824644/101‑3 (MQ=255)
                                               ctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGCGGTGaa                                                       >  2:204138/1‑101 (MQ=255)
                                               ctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaa                                                       >  1:1177905/1‑101 (MQ=255)
                                                tcCCGAAATAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                      <  1:6434813/101‑1 (MQ=255)
                                                tcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                      >  2:5602671/1‑101 (MQ=255)
                                                 ccccAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaacg                                                     >  2:772956/1‑99 (MQ=255)
                                                 cccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                     <  1:6279169/101‑1 (MQ=255)
                                                  cccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAAcgt                                                    <  2:4349693/98‑1 (MQ=255)
                                                  cccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAAcgt                                                    >  1:1400496/4‑101 (MQ=255)
                                                   cGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtgtg                                                   <  2:3202138/101‑1 (MQ=255)
                                                      aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtt                                                <  2:2578389/101‑2 (MQ=255)
                                                      aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggtg                                                <  1:4712779/101‑1 (MQ=255)
                                                      aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggtg                                                >  2:146483/1‑101 (MQ=255)
                                                       agagTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtta                                               >  1:4430806/1‑99 (MQ=255)
                                                       aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                               >  1:1984901/1‑101 (MQ=255)
                                                        gAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                              <  1:992455/101‑1 (MQ=255)
                                                        gAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                              >  1:2324030/1‑101 (MQ=255)
                                                         aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                             <  1:564573/101‑1 (MQ=255)
                                                           tCCGCATACCAGGAAGGGCGCTGGGAAACTCTGCCCTTTCAGCGGGCCATCGTGAATGCGATGGGTAGCGACTACATCCGTGAGGTGAATGTGGTTAAGTc                                           <  2:6646921/101‑1 (MQ=255)
                                                             ggcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTg                                         <  2:3107292/97‑1 (MQ=255)
                                                              gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                        >  1:5386665/1‑101 (MQ=255)
                                                               cctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                       >  2:6221849/3‑101 (MQ=255)
                                                               cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                       <  1:6052403/101‑1 (MQ=255)
                                                               cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                       >  2:5539384/1‑101 (MQ=255)
                                                                ctcccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                      <  1:3148167/98‑1 (MQ=255)
                                                                aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTTCATCCGTGAGGTGAATGTGGTGAAGTCTGccc                                      >  1:1954010/1‑101 (MQ=255)
                                                                 tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                     >  2:5211457/1‑101 (MQ=255)
                                                                 tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGTATGTGGTGAAGACTGCCCg                                     >  2:4871184/1‑101 (MQ=255)
                                                                 tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCg                                     <  2:1279075/101‑1 (MQ=255)
                                                                   ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                   <  2:3388504/101‑1 (MQ=255)
                                                                    cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                  >  2:4712346/1‑101 (MQ=255)
                                                                    cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGGCAACGTTGTTAAGTCTGCCCgtgt                                  >  2:17334/1‑101 (MQ=255)
                                                                    cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                  >  2:290259/1‑101 (MQ=255)
                                                                     aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGACTGTGGTGAAGTCTGCCCGTGTc                                 >  1:3693/1‑101 (MQ=255)
                                                                     aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                 >  1:6026586/1‑101 (MQ=255)
                                                                     aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                 <  1:5120532/101‑1 (MQ=255)
                                                                     aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                 <  2:1649876/101‑1 (MQ=255)
                                                                        aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGt                              >  1:3091026/1‑101 (MQ=255)
                                                                        aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGt                              >  2:2804914/1‑101 (MQ=255)
                                                                         aGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                             >  1:6938726/1‑101 (MQ=255)
                                                                              gcTGGGAAACACTGCCCTTCCAGCGGGCAATCATGAACGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAttac                        <  1:3190595/101‑3 (MQ=255)
                                                                               cTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCa                       >  2:3878027/1‑101 (MQ=255)
                                                                                 gggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                     <  1:3417015/101‑1 (MQ=255)
                                                                                   gAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                   >  2:4253837/1‑101 (MQ=255)
                                                                                    aaaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtg                  >  2:3720965/1‑101 (MQ=255)
                                                                                     aaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgc                 <  2:6414509/101‑1 (MQ=255)
                                                                                      acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                <  2:1513637/101‑1 (MQ=255)
                                                                                       cncTNNCCTTTCAGNGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctt               >  2:1136707/1‑100 (MQ=255)
                                                                                       cgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg               <  1:3925958/98‑1 (MQ=255)
                                                                                       cgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg               >  1:3787273/4‑101 (MQ=255)
                                                                                         cTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct             <  1:3366684/101‑1 (MQ=255)
                                                                                          tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg            <  2:2997635/101‑1 (MQ=255)
                                                                                           gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTgg           >  1:4008037/1‑101 (MQ=255)
                                                                                           gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTgg           >  2:3607394/1‑101 (MQ=255)
                                                                                            cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg          >  2:5434260/1‑101 (MQ=255)
                                                                                             cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGgt         <  1:1859829/99‑1 (MQ=255)
                                                                                                ttCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt      <  2:4070283/101‑1 (MQ=255)
                                                                                                ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt      >  1:6179956/1‑101 (MQ=255)
                                                                                                  cAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa    >  2:181813/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATG  >  NC_000913/580647‑580847

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