Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,4250TA47.7% 169.2 / 255.4 155G123G (GGT→GGA) nohDDLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base A (35/37);  total (77/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 183 total aligned reads displayed.

CGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAA  >  NC_000913/580325‑580511
                                                                                                    |                                                                                      
nnngTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                                        <  2:10028977/98‑1 (MQ=17)
cgcgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                                        >  2:5294960/1‑101 (MQ=32)
cgcgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                                        <  1:9099997/101‑1 (MQ=32)
 gcgTGCGCAGTCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg                                                                                       <  1:10565678/101‑1 (MQ=17)
 gcgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg                                                                                       <  1:10141900/101‑1 (MQ=32)
 gcgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg                                                                                       <  2:1833470/101‑1 (MQ=32)
  ggTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa                                                                                      <  1:6705094/100‑1 (MQ=17)
  cgagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga                                                                                      <  1:14849563/98‑1 (MQ=2)
  cgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa                                                                                      >  1:5364050/1‑101 (MQ=33)
  cgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGGGa                                                                                      >  1:11589570/1‑101 (MQ=18)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGGGaa                                                                                     >  2:5119321/3‑101 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     <  2:6885149/99‑1 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     <  2:2472814/99‑1 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     <  1:4651240/99‑1 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     >  2:1789826/3‑101 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     >  1:6239744/3‑101 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     >  2:876441/3‑101 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     >  2:362630/3‑101 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     <  1:11191606/99‑1 (MQ=2)
   gagcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTAGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                                     <  1:12205823/99‑1 (MQ=2)
   gTGCGCAGGCCGATGCACAGGAGCTTAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCTGATCGCAGGTGaa                                                                                     >  1:6089014/1‑101 (MQ=12)
   gTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGtnnn                                                                                     >  2:9222576/1‑98 (MQ=18)
   gTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGaa                                                                                     <  1:4321418/101‑1 (MQ=32)
    agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaag                                                                                    <  2:6709630/100‑2 (MQ=2)
     nnncAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt                                                                                   <  2:9436404/98‑3 (MQ=2)
     gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                                                   >  1:11884436/1‑101 (MQ=32)
     gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGGGAAAt                                                                                   >  1:522527/1‑101 (MQ=17)
     gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGAGCTGTCGCAGATGGCCGgagannn                                                                                   >  2:10499295/1‑98 (MQ=2)
      cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGaannn                                                                                  >  2:14297658/1‑98 (MQ=17)
      cgcAGGCAGATTCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                                                  <  1:12146386/101‑4 (MQ=2)
      cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                                                  <  1:2976614/101‑4 (MQ=2)
      cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                                                  <  1:11450929/101‑4 (MQ=2)
       nnnggCTGATGCACAGGAGCTGAATAATGCCATAGACTCCGCTGAAGTGGTTGAATCCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg                                                                                 <  2:9900535/98‑1 (MQ=11)
       gtaGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCCCGGATCGCCGGTGaaannn                                                                                 >  2:12008314/3‑98 (MQ=11)
       gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag                                                                                 <  1:6745130/101‑5 (MQ=2)
       gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag                                                                                 >  1:13216960/1‑97 (MQ=2)
       gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag                                                                                 >  1:14433622/1‑97 (MQ=2)
       gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCTCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg                                                                                 <  2:441586/101‑1 (MQ=17)
        nnngCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTgt                                                                                <  2:14146999/98‑2 (MQ=17)
        nnngCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc                                                                                <  2:8192321/98‑1 (MQ=17)
        cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc                                                                                >  1:6267506/1‑101 (MQ=32)
        cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc                                                                                >  1:3057199/1‑101 (MQ=32)
        cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGGGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGc                                                                                >  1:3868840/1‑101 (MQ=17)
        cAGGCAGATGCCCCGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtnnn                                                                                >  2:13835375/1‑96 (MQ=2)
        cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtagc                                                                                >  1:311630/1‑96 (MQ=2)
         aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAttnnn                                                                               >  2:14431444/1‑98 (MQ=17)
         aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGGGAAATTGcc                                                                               >  2:5919710/1‑101 (MQ=17)
         aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGcc                                                                               >  1:13031362/1‑101 (MQ=2)
          nnncGATGCACAGGAGCTGAAAAATGCGAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTCAGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                              <  2:13368971/98‑1 (MQ=11)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCTCAGGGGAAATTGccg                                                                              >  2:1436330/1‑100 (MQ=12)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTgnnn                                                                              >  2:12305750/1‑98 (MQ=17)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                              >  2:2340343/1‑101 (MQ=32)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                              >  2:84092/1‑101 (MQ=32)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                              >  1:94793/1‑101 (MQ=32)
          ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCa                                                                              <  1:14707542/101‑1 (MQ=32)
          ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCa                                                                              <  1:4266643/101‑1 (MQ=2)
          ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCa                                                                              <  1:14060479/101‑1 (MQ=2)
          ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCa                                                                              >  1:12764906/1‑101 (MQ=2)
          ggCAGATGCACAGTAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCa                                                                              <  2:4218336/101‑1 (MQ=2)
           gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAg                                                                             <  1:7958113/98‑1 (MQ=2)
           gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGcnnn                                                                             >  2:7326467/1‑98 (MQ=17)
           gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAg                                                                             <  1:8061578/101‑1 (MQ=32)
           gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAg                                                                             >  1:2479230/1‑101 (MQ=32)
           gCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAg                                                                             >  2:3015329/1‑101 (MQ=32)
            nnnaTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAgc                                                                            <  2:9543134/98‑2 (MQ=2)
            nnnaTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGCAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGt                                                                            <  2:8559692/98‑1 (MQ=9)
            ctgATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGccnnn                                                                            >  2:14690368/3‑98 (MQ=2)
            ccGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGccnnn                                                                            >  2:7926912/1‑98 (MQ=17)
            ccGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGt                                                                            >  1:7706752/1‑101 (MQ=32)
            cagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAgc                                                                            >  2:541596/3‑100 (MQ=2)
            cagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAgc                                                                            >  1:5547348/3‑100 (MQ=2)
             cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTa                                                                           <  1:10638377/101‑1 (MQ=32)
             cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTa                                                                           >  2:5821795/1‑101 (MQ=32)
             cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGGGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTa                                                                           >  1:9335693/1‑101 (MQ=17)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAgnnn                                                                          >  2:12068084/1‑98 (MQ=2)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAt                                                                          >  1:10203466/1‑101 (MQ=2)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAt                                                                          >  1:10258909/1‑101 (MQ=2)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAt                                                                          <  1:12902840/101‑1 (MQ=2)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAt                                                                          <  1:3179911/101‑1 (MQ=2)
              gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTATGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAt                                                                          >  1:8532097/1‑101 (MQ=2)
              gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                          >  1:6129444/1‑101 (MQ=32)
              gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                          <  1:2052300/101‑1 (MQ=32)
              gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                          <  2:4191965/101‑1 (MQ=32)
              gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                          >  2:4649297/1‑101 (MQ=32)
              gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAt                                                                          <  2:2557158/101‑1 (MQ=32)
               aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAtt                                                                         >  1:11732739/1‑101 (MQ=2)
               aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCAtt                                                                         <  1:1996131/101‑1 (MQ=2)
               aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAtt                                                                         <  1:7326467/101‑1 (MQ=31)
               aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTAtt                                                                         >  1:7610062/1‑101 (MQ=31)
               aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGAGAAATTGCCAGTAtt                                                                         >  2:601216/1‑101 (MQ=14)
                tGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtc                                                                        >  2:6063666/1‑101 (MQ=2)
                tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGCAATTGCCAGTATtc                                                                        >  1:13999911/1‑101 (MQ=11)
                tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCATTATtc                                                                        >  2:2195465/1‑101 (MQ=11)
                tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATtc                                                                        >  1:11140228/1‑101 (MQ=11)
                tGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATtc                                                                        <  1:9641718/101‑1 (MQ=11)
                 gcccaGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct                                                                       <  1:2234593/98‑1 (MQ=0)
                 gcccaGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGAGGATGGCCGGAGAAGTAGCCAGCATtct                                                                       <  2:6176097/98‑1 (MQ=0)
                      ggtgCTGAAAAATGCCAGAGCCTCCGCTGAAGTGGTGTAAACCGCATTCTGTTCTTTCGTGCGGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACg                                                                  <  1:1859171/98‑1 (MQ=1)
                                                                        gcacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACAt                <  1:3538181/96‑1 (MQ=0)
                                                                          acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACAtat              >  2:38569/4‑99 (MQ=0)
                                                                                                    |                                                                                      
CGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAA  >  NC_000913/580325‑580511

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