Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,4190CG70.7% 108.9 / 128.6 140I121M (ATC→ATGnohBbacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base C (20/21);  new base G (55/44);  total (75/65)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 149 total aligned reads displayed.

CGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC  >  NC_000913/580327‑580491
                                                                                            |                                                                        
cgTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATc                                                                          <  2:5200327/93‑1 (MQ=33)
 gTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCg                                                                         >  2:5567807/1‑93 (MQ=33)
  tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc                                                                        >  1:4116438/1‑93 (MQ=33)
  tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc                                                                        >  1:5173248/1‑93 (MQ=33)
  tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc                                                                        >  1:1369337/1‑93 (MQ=34)
  tGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGc                                                                        >  2:5834260/1‑93 (MQ=34)
  agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGc                                                                        >  1:2481046/2‑93 (MQ=2)
  agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGc                                                                        <  1:2023388/92‑1 (MQ=2)
  agcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGCTCGTGCTGTCGCGGATGGc                                                                        >  1:2722504/2‑93 (MQ=2)
  agcgcAGGCAGATGCCCAGGAACTAAAAAATGAAAAAGATTCCGCTGAGGTGGTGGAAACCGCATACTGCAGGTTCGTCCTGTCGCGGATGGc                                                                        >  1:1174163/2‑93 (MQ=2)
   gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa                                                                       <  2:1694490/93‑1 (MQ=33)
   gcgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCa                                                                       <  2:2796365/93‑1 (MQ=33)
   gcgcAGGCAGATGCCCAGGACCTAAAAAATGCAAAAGATTCCGCTGAGGCGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc                                                                       <  1:1760922/93‑2 (MQ=2)
   gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc                                                                       <  1:2332955/93‑2 (MQ=2)
   gcgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGcc                                                                       >  2:4665792/1‑92 (MQ=2)
    ggcAGGCCCATGCACAGGAGCTGAAAAATGCCAGAGGCTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg                                                                      <  1:4786897/92‑1 (MQ=255)
    cgcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAg                                                                      <  1:6032636/93‑1 (MQ=33)
    cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGccg                                                                      <  1:6127317/93‑3 (MQ=2)
    cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGccg                                                                      <  1:633751/93‑3 (MQ=2)
    cgcAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGccg                                                                      <  2:1008063/93‑3 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  1:4737405/1‑93 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  2:3267616/1‑93 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  2:1946554/1‑93 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  1:3496974/1‑93 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     <  1:5750297/93‑1 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  2:3141497/1‑93 (MQ=2)
     gcaggcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  2:2296616/1‑93 (MQ=2)
     gcaggcagATCCCCAGGAACTAAAAAATGCAAAAGATTCGGCTGAGGGGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgg                                                                     >  1:1926740/1‑93 (MQ=2)
     gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg                                                                     <  1:3072843/93‑1 (MQ=33)
     gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg                                                                     <  1:5229123/93‑1 (MQ=33)
     gcAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgg                                                                     >  1:2778250/1‑93 (MQ=34)
      cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                    <  2:42033/93‑1 (MQ=33)
      cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                    <  1:2746960/93‑1 (MQ=33)
      cAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGt                                                                    <  2:42034/93‑1 (MQ=33)
      cAGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCgga                                                                    <  1:1464021/93‑2 (MQ=2)
       aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAgggg                                                                   >  2:31587/1‑91 (MQ=21)
       aGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTCGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTg                                                                   <  1:5195327/93‑1 (MQ=18)
       aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag                                                                   >  1:1124569/1‑91 (MQ=2)
       aGGCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCggag                                                                   >  2:1999040/1‑91 (MQ=2)
        ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa                                                                  <  2:1074126/93‑1 (MQ=33)
        ggCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGa                                                                  >  2:4466810/1‑93 (MQ=34)
        ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga                                                                  >  2:3672634/1‑93 (MQ=2)
        ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga                                                                  >  1:1696028/1‑93 (MQ=2)
        ggCAGATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGgaga                                                                  <  2:5556010/93‑1 (MQ=2)
         gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                 <  2:5987613/90‑1 (MQ=2)
         gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                 <  2:3277997/90‑1 (MQ=2)
         gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                 >  2:1170900/4‑93 (MQ=2)
         gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                 <  2:3728962/90‑1 (MQ=2)
         gcagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaa                                                                 <  1:1068175/90‑1 (MQ=2)
         gcagATGCCCAGAAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTAGTGCTGTCGCGGATGGCCGGAGaa                                                                 <  2:5594941/90‑1 (MQ=2)
          cagATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaag                                                                <  1:1634474/91‑2 (MQ=2)
           cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                               >  2:898072/1‑93 (MQ=33)
           cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                               <  1:4186793/93‑1 (MQ=33)
           cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                               <  2:4016569/93‑1 (MQ=33)
           cGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAt                                                               <  2:3388832/93‑1 (MQ=33)
           agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt                                                               <  2:916488/92‑3 (MQ=2)
           agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt                                                               >  1:3775378/2‑91 (MQ=2)
           agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagt                                                               >  1:57807/2‑91 (MQ=2)
           agATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaaga                                                               >  1:5002648/2‑91 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTTGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  2:2278137/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAGCCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  1:3830314/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  2:1500911/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  2:5604658/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  1:2687707/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              <  1:3199024/93‑4 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  2:5294228/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  1:3415153/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  1:1963483/1‑90 (MQ=2)
            gATGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagta                                                              >  1:5021116/1‑90 (MQ=2)
            gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt                                                              >  1:925632/1‑93 (MQ=33)
            gATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAAtt                                                              >  2:4511716/1‑93 (MQ=33)
             aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGGGaagcag                                                             >  1:1796496/1‑89 (MQ=2)
             aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag                                                             >  2:2353909/1‑89 (MQ=2)
             aTGCCCAGGAACTAAAAAATGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGaagtag                                                             >  1:1796495/1‑89 (MQ=2)
             aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg                                                             >  1:2018496/1‑93 (MQ=32)
             aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTg                                                             <  1:3412079/93‑1 (MQ=32)
             aTGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGGGAAATTg                                                             >  2:1326971/1‑93 (MQ=17)
                       actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct                                                   >  1:5555608/5‑93 (MQ=0)
                       actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct                                                   >  1:3554127/5‑93 (MQ=0)
                       actaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtct                                                   >  1:4376165/5‑93 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  >  2:4390699/4‑92 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  >  2:473382/4‑92 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  <  1:4547088/90‑2 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  >  2:542388/4‑92 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  <  2:1227051/90‑2 (MQ=0)
                        ctaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtctt                                                  >  2:5960928/4‑92 (MQ=0)
                         taaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATtcttg                                                 <  2:34308/91‑3 (MQ=0)
                                                gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAg                          >  1:6007665/4‑93 (MQ=0)
                                                   gtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                       <  2:4594462/93‑1 (MQ=0)
                                                      gtcgAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcggcgc                    <  2:5427547/90‑2 (MQ=0)
                                                                       cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa   <  2:2621545/89‑1 (MQ=0)
                                                                       cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGaaaa   <  1:4425721/89‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGGCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  <  2:3922417/90‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  <  1:1075092/90‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  >  2:3429332/4‑93 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  <  2:3401072/90‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  <  1:4579125/90‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  <  2:2459802/90‑1 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAAc  >  2:1418671/4‑93 (MQ=0)
                                                                        acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGAAGTTCCTCTGTCGGTTAAGCGGCGCTTCCCGGAGCTGGAAAAc  >  2:5751229/4‑93 (MQ=0)
                                                                                            |                                                                        
CGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAC  >  NC_000913/580327‑580491

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