Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7920GC62.4% 449.1 / inf 444V72V (GTG→GTCybcXhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (90/77);  new base C (129/148);  total (219/225)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 458 total aligned reads displayed.

AGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAT  >  NC_000913/580702‑580884
                                                                                          |                                                                                            
aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGaa                                                                                            >  1:4112787/1‑93 (MQ=255)
aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGaa                                                                                            >  2:3874677/1‑93 (MQ=255)
    tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtgtg                                                                                        >  2:3022393/1‑93 (MQ=255)
        cctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtta                                                                                    >  2:562331/3‑91 (MQ=255)
        cctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGATAACGATTACATACGTGAGGTCAACgttgtta                                                                                    <  1:770915/91‑3 (MQ=255)
         ctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTaa                                                                                   >  2:1523234/2‑93 (MQ=255)
         aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                                                                   <  1:4308913/93‑1 (MQ=255)
         aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                                                                   <  1:590949/93‑1 (MQ=255)
         aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGaa                                                                                   >  1:3532399/1‑93 (MQ=255)
          tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                                  >  2:4610706/1‑93 (MQ=255)
          tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                                  >  1:2778570/1‑93 (MQ=255)
          tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGGGAATGTGGCGagg                                                                                  >  2:3751812/1‑91 (MQ=255)
           aCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                                                                 >  2:2214592/1‑93 (MQ=255)
           aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                                                 <  2:295663/93‑1 (MQ=255)
             cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                               <  1:223692/93‑1 (MQ=255)
             cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                               <  2:1236823/93‑1 (MQ=255)
             cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                               >  2:3572030/1‑93 (MQ=255)
             cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                               >  2:663609/1‑93 (MQ=255)
             cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                               <  1:3610743/93‑1 (MQ=255)
               ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                                                             <  1:710341/93‑1 (MQ=255)
               ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                                                             >  2:512432/1‑93 (MQ=255)
               ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGAGAATGTGGTGAAGTCTGc                                                                             >  1:2305847/1‑93 (MQ=255)
                 aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                                                           <  2:1898819/93‑1 (MQ=255)
                  aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                          >  2:3494410/1‑93 (MQ=255)
                  aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                          >  2:3925055/1‑93 (MQ=255)
                  aGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCg                                                                          >  1:4221073/1‑93 (MQ=255)
                   gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                         >  1:1068624/1‑93 (MQ=255)
                   gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                         <  1:2430950/93‑1 (MQ=255)
                     gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                                       <  2:953555/93‑1 (MQ=255)
                     gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                                       <  1:3826544/93‑1 (MQ=255)
                      cgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                                                                      <  1:3760906/93‑2 (MQ=255)
                      cgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                                                                      <  1:149983/93‑2 (MQ=255)
                      cgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                                                                      <  1:4357130/93‑2 (MQ=255)
                       gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                                     >  2:3198044/1‑93 (MQ=255)
                        cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                                    <  1:1688935/93‑1 (MQ=255)
                        cTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                                                                    >  1:141191/1‑93 (MQ=255)
                          gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                                  >  2:2316051/1‑91 (MQ=255)
                          gggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                                  <  1:2937298/93‑1 (MQ=255)
                           ggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTa                                                                 <  1:321317/93‑1 (MQ=255)
                            gAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGATTACATACGTGAGGTCAACGTTGGTAAGTCTGCCCGTGTTGGCTAt                                                                >  1:4709572/1‑93 (MQ=255)
                              aaCGCTGCCCTTTAAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTc                                                              <  1:625631/93‑1 (MQ=255)
                               acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                                                             <  2:2704109/93‑1 (MQ=255)
                                  cTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                                          >  2:602390/1‑93 (MQ=255)
                                   tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                         <  2:1005702/93‑1 (MQ=255)
                                   tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                         <  2:2855755/93‑1 (MQ=255)
                                    gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                                        <  2:1953511/93‑1 (MQ=255)
                                     cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtg                                                       >  1:4175131/1‑93 (MQ=255)
                                      cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgt                                                      <  1:1276951/91‑2 (MQ=255)
                                       cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                                                     <  1:3952690/93‑1 (MQ=255)
                                       cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                                                     <  1:2490313/93‑1 (MQ=255)
                                        tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                    <  2:4674187/93‑1 (MQ=255)
                                          tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                                  <  1:2901969/93‑1 (MQ=255)
                                          tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                                  <  1:4286051/93‑1 (MQ=255)
                                           cAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg                                                 <  1:852035/93‑1 (MQ=255)
                                               gggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtg                                             <  2:3337231/93‑1 (MQ=255)
                                                ggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                            <  1:2805891/93‑1 (MQ=255)
                                                 gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                           >  1:2044297/4‑93 (MQ=255)
                                                 gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                           <  1:3489920/90‑1 (MQ=255)
                                                 gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                           >  2:1721224/4‑93 (MQ=255)
                                                 gCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                                           >  2:3106202/1‑93 (MQ=255)
                                                   gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                         <  1:4424550/92‑1 (MQ=255)
                                                   gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                         <  2:256417/92‑1 (MQ=255)
                                                   gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                         >  2:2770461/2‑93 (MQ=255)
                                                    gtcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                        <  1:4056165/92‑1 (MQ=255)
                                                    atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                        >  2:213562/1‑93 (MQ=255)
                                                    atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                        >  2:2480088/1‑93 (MQ=255)
                                                    atcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTAt                                        >  2:4401400/1‑93 (MQ=255)
                                                     tcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                                       <  2:344253/93‑1 (MQ=255)
                                                     tcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATg                                       >  1:1762410/1‑93 (MQ=255)
                                                        tGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcgt                                    >  2:4044173/1‑91 (MQ=255)
                                                         gAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTa                                   <  1:1216459/93‑1 (MQ=255)
                                                            tGCGATGGGGCACGATTACATACGTGAGGTAAACGTTGTTAAGTCTGCCCGTGTTGGCTACTGCAAAATGTTGCTGGGTGTTTATGCGTAttt                                >  2:1175388/1‑93 (MQ=255)
                                                            tGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt                                <  2:3627498/93‑1 (MQ=255)
                                                              cGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTta                              >  2:2100837/1‑93 (MQ=255)
                                                               gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtat                             >  1:1856414/1‑93 (MQ=255)
                                                               gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtat                             >  1:1945598/1‑93 (MQ=255)
                                                               gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                             >  2:819612/1‑93 (MQ=255)
                                                               gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                             <  1:3373373/93‑1 (MQ=255)
                                                                aTGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtatt                            <  1:211259/93‑2 (MQ=255)
                                                                 tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                           >  2:134447/1‑91 (MQ=255)
                                                                 tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                           <  2:4195555/93‑3 (MQ=255)
                                                                  gggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga                          <  2:4248178/93‑1 (MQ=255)
                                                                   ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATagag                         >  1:3254156/1‑93 (MQ=255)
                                                                   ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATagac                         >  2:1494498/1‑92 (MQ=255)
                                                                     cAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCa                       >  1:1328767/1‑93 (MQ=255)
                                                                       acGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCaca                     >  1:2547871/2‑91 (MQ=255)
                                                                        cGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                         <  1:920134/72‑1 (MQ=255)
                                                                        cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATaa                    >  1:1033020/1‑93 (MQ=255)
                                                                          attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc                  <  1:2638339/91‑1 (MQ=255)
                                                                          attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc                  >  2:1067941/3‑93 (MQ=255)
                                                                           ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                 <  2:2191149/92‑1 (MQ=255)
                                                                           ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                 <  1:498797/92‑1 (MQ=255)
                                                                           ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                 <  1:4008549/92‑1 (MQ=255)
                                                                           ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                 >  1:3662498/2‑93 (MQ=255)
                                                                            tACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcag                <  2:1828506/93‑1 (MQ=255)
                                                                               atacGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcga             >  2:2374236/4‑92 (MQ=255)
                                                                                tacGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcgaa            >  2:1599603/3‑91 (MQ=255)
                                                                                   gTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaatt         >  2:1122976/1‑91 (MQ=255)
                                                                                         tGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTa   >  2:2051433/1‑93 (MQ=255)
                                                                                          gAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAt  >  1:4620032/1‑93 (MQ=255)
                                                                                          |                                                                                            
AGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAT  >  NC_000913/580702‑580884

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