Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7980GT68.4% 903.7 / inf 812V74V (GTG→GTTybcXhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (132/124);  new base T (256/299);  total (388/424)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 856 total aligned reads displayed.

TCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGC  >  NC_000913/580706‑580890
                                                                                            |                                                                                            
  ggcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgt                                                                                            >  1:858724/5‑93 (MQ=255)
  cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggt                                                                                            <  1:4004202/93‑1 (MQ=255)
  cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggt                                                                                            <  1:4122857/93‑1 (MQ=255)
   gcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtt                                                                                           >  2:1484801/4‑92 (MQ=255)
    cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                                          >  1:1389489/1‑93 (MQ=255)
      tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                                        >  2:5216231/1‑93 (MQ=255)
       aCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                                                                       >  2:262400/1‑93 (MQ=255)
        ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                                                                      <  1:2870047/93‑1 (MQ=255)
        ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                                                                      <  2:1356246/93‑1 (MQ=255)
         cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATTGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                                     <  2:2276532/93‑1 (MQ=255)
           ggAAGGGCGCTGGGAAACGTTGCGGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                                                                   <  2:4195661/93‑1 (MQ=255)
            gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                                  >  2:5670098/1‑93 (MQ=255)
            gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                                  <  1:929888/93‑1 (MQ=255)
             aaGGGCGCTGGGTAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                                                                 <  1:5479056/93‑1 (MQ=255)
             aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                                                                 >  1:49985/1‑93 (MQ=255)
             aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                                                                 >  1:4440769/1‑93 (MQ=255)
              aGGGCGCTGGTAAACGTTGCCGTTTCAACGCGCAATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                                <  2:527531/93‑1 (MQ=255)
               gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                               >  2:4891833/1‑93 (MQ=255)
                ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtg                                                                              >  2:3608721/1‑93 (MQ=255)
                 gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                                             <  2:1620366/93‑1 (MQ=255)
                 gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                             <  2:915693/93‑1 (MQ=255)
                  cgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                                                                            <  2:393758/93‑2 (MQ=255)
                  cgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCATGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                                                            <  1:2357797/93‑1 (MQ=255)
                   gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                                           <  2:5031466/93‑1 (MQ=255)
                    cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                                          <  1:247608/93‑1 (MQ=255)
                    cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                                          <  2:1273341/93‑1 (MQ=255)
                     tGGGAACCGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                         <  2:4323062/93‑2 (MQ=255)
                     tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                         <  1:1081075/93‑2 (MQ=255)
                     tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                         >  2:4748449/1‑92 (MQ=255)
                     tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                         >  1:3758554/1‑92 (MQ=255)
                      gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                                        >  1:264827/1‑91 (MQ=255)
                      gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                                        >  1:2793673/1‑91 (MQ=255)
                      gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggct                                                                        >  2:295913/1‑91 (MQ=255)
                           acgcTGCCCTTTCAGCGGACCATCATGAATGAGATGGGCAGCGACTGCATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                                                                   <  2:3068764/90‑1 (MQ=255)
                             acTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTTTTccca                                                                 >  2:4414186/1‑91 (MQ=255)
                               tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                               >  2:3778797/1‑93 (MQ=255)
                               tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                               <  2:1844500/93‑1 (MQ=255)
                               tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                               >  1:1994607/1‑93 (MQ=255)
                                gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                                              <  2:6087248/93‑1 (MQ=255)
                                gCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                                                              <  1:704211/93‑1 (MQ=255)
                                 ccgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                                                             <  1:5382829/90‑1 (MQ=255)
                                  cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgt                                                            >  2:3190180/3‑92 (MQ=255)
                                  ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgc                                                            >  2:2430491/1‑93 (MQ=255)
                                    tttCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                                          >  1:3110614/1‑93 (MQ=255)
                                    tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                          <  2:1910281/93‑1 (MQ=255)
                                    tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                          <  1:1518604/93‑1 (MQ=255)
                                    tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                          <  1:4343522/93‑1 (MQ=255)
                                    tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTGAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                          >  2:5935055/1‑93 (MQ=255)
                                     ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                         <  1:3855858/93‑1 (MQ=255)
                                     ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                         <  2:5858789/93‑1 (MQ=255)
                                     ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                         <  1:2290077/93‑1 (MQ=255)
                                      tCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                                        <  2:1170766/93‑1 (MQ=255)
                                      tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                                        <  1:450599/93‑1 (MQ=255)
                                          cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgt                                                    <  1:783919/93‑1 (MQ=255)
                                            cgcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGgtgt                                                  <  1:4858610/89‑1 (MQ=255)
                                             gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                                 >  1:3948452/4‑93 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               <  1:1275846/92‑1 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               >  1:3232376/2‑93 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               <  2:3831498/92‑1 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               <  1:1764934/92‑1 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               <  2:5792067/92‑1 (MQ=255)
                                               gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                               <  1:5751229/92‑1 (MQ=255)
                                                atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                              >  1:272634/1‑93 (MQ=255)
                                                 tcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                                             <  2:2071994/93‑1 (MQ=255)
                                                 tcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTTTTTATg                                             >  2:5998980/1‑93 (MQ=255)
                                                  catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGATGCTGGGTGTTTATGc                                            >  2:5493171/1‑93 (MQ=255)
                                                   atGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcg                                           <  2:4075620/93‑2 (MQ=255)
                                                   atGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc                                           <  1:4594017/93‑1 (MQ=255)
                                                    tGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcgt                                          >  1:1634679/1‑91 (MQ=255)
                                                    tGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcgt                                          <  1:5841695/93‑3 (MQ=255)
                                                      aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                        <  2:535964/93‑1 (MQ=255)
                                                           gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtat                                   <  1:3420956/93‑1 (MQ=255)
                                                           gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                                   <  2:3539704/93‑1 (MQ=255)
                                                            aTGGGTAACGATTACATACGTGAGCTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATCCAAAAATGTTGCTGGGTGTTTATGCGTATTTtatt                                  <  1:6010989/93‑2 (MQ=255)
                                                             tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                                 >  2:3420056/1‑91 (MQ=255)
                                                              gggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga                                >  1:523771/1‑93 (MQ=255)
                                                                 taacGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCa                             <  2:5132580/90‑1 (MQ=255)
                                                                 cAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCa                             <  2:5602137/93‑1 (MQ=255)
                                                                    cGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACaa                          <  2:1690561/93‑1 (MQ=255)
                                                                    cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATaa                          <  1:3585972/93‑1 (MQ=255)
                                                                     gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAag                         >  1:5595329/4‑93 (MQ=255)
                                                                     gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAag                         >  1:5095103/4‑93 (MQ=255)
                                                                      aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                        >  2:779360/1‑93 (MQ=255)
                                                                      aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                        <  2:2800710/93‑1 (MQ=255)
                                                                       ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                       <  1:6139567/92‑1 (MQ=255)
                                                                       cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca                       <  1:254673/93‑1 (MQ=255)
                                                                        tACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcag                      >  1:1880126/1‑93 (MQ=255)
                                                                         aCATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcagc                     <  2:2722107/93‑1 (MQ=255)
                                                                          cATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcg                    >  2:1801179/1‑93 (MQ=255)
                                                                            tacGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcgaa                  >  2:4601256/3‑91 (MQ=255)
                                                                            tCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTGCTTTATAGCGCATAAGCAGCCCa                  >  2:1022172/1‑93 (MQ=255)
                                                                             acGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaa                 <  1:216584/92‑1 (MQ=255)
                                                                             acGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaa                 <  2:330024/92‑1 (MQ=255)
                                                                             acGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaa                 <  1:4715231/92‑1 (MQ=255)
                                                                             acGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGAATTTTATTCAGCACAAGCAGCGaaa                 <  1:4229147/92‑1 (MQ=255)
                                                                               gTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaatt               >  2:743906/1‑91 (MQ=255)
                                                                                   ggTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt           >  1:303802/1‑93 (MQ=255)
                                                                                    gTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCtt          <  1:520536/93‑1 (MQ=255)
                                                                                        aTGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTc      <  1:3338197/93‑1 (MQ=255)
                                                                                         tgtgGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTcc     >  1:5918772/1‑93 (MQ=255)
                                                                                            |                                                                                            
TCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGC  >  NC_000913/580706‑580890

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