Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7950TC68.2% 1050.9 / inf 837N73N (AAT→AACybcXhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (135/131);  new base C (268/303);  total (403/434)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 864 total aligned reads displayed.

GAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCC  >  NC_000913/580703‑580887
                                                                                            |                                                                                            
 aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                             >  1:1470383/1‑93 (MQ=255)
    ccGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                                                          <  1:30611/93‑1 (MQ=255)
      gcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtt                                                                                        >  2:1484801/4‑92 (MQ=255)
       cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                                       >  1:1389489/1‑93 (MQ=255)
         tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                                     >  1:330287/1‑93 (MQ=255)
         tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                                     <  1:2418539/93‑1 (MQ=255)
         tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                                     >  2:5216231/1‑93 (MQ=255)
          aCCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGt                                                                                    <  1:927241/93‑1 (MQ=255)
           ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                                                                   <  1:690084/93‑1 (MQ=255)
           ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                                                                   <  1:4433234/93‑1 (MQ=255)
             aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTg                                                                                 <  1:3225282/93‑1 (MQ=255)
             aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTg                                                                                 >  1:5683718/1‑93 (MQ=255)
              ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                                <  1:2761775/93‑1 (MQ=255)
               gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                               <  2:1266038/93‑1 (MQ=255)
                 aGGGCGCTGGTAAACGTTGCCGTTTCAACGCGCAATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                             <  2:527531/93‑1 (MQ=255)
                 aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                             >  1:2516674/1‑93 (MQ=255)
                  gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                            >  2:4891833/1‑93 (MQ=255)
                   ggCGCTGGGAAGCGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCtgtg                                                                           <  1:5576743/93‑1 (MQ=255)
                   ggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                                                           <  1:5366537/93‑1 (MQ=255)
                   ggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                                                           <  1:5796065/93‑1 (MQ=255)
                    gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                          <  2:915693/93‑1 (MQ=255)
                     cgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGCGAATGTGGTGAAGTCTGCCCGGGTc                                                                         >  1:1634186/1‑93 (MQ=255)
                      gcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttg                                                                        >  1:3825998/1‑91 (MQ=255)
                      gcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttg                                                                        <  1:749350/93‑3 (MQ=255)
                      gcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttg                                                                        >  1:4833824/1‑91 (MQ=255)
                       cTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                                                                       <  1:4705273/93‑1 (MQ=255)
                        tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                      <  1:1081075/93‑2 (MQ=255)
                        tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                                      >  2:4748449/1‑92 (MQ=255)
                          ggNAACACTGCCCTTTCAGCGNNCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                                    <  2:1153337/93‑1 (MQ=255)
                              acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTcc                                                                >  2:4523633/1‑93 (MQ=255)
                                acTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTTTTccca                                                              >  2:4414186/1‑91 (MQ=255)
                                 ttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCaaaa                                                             >  1:3571209/2‑93 (MQ=255)
                                  tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                            <  1:270906/93‑1 (MQ=255)
                                  tGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                            >  2:5935558/1‑93 (MQ=255)
                                  ctcccGTTCAGCGGGCCGTCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                            <  1:4119087/91‑1 (MQ=255)
                                   gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                                           >  2:3929716/1‑93 (MQ=255)
                                     cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgt                                                         <  1:5437726/91‑2 (MQ=255)
                                     ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGGCCGTGTCGGTTATTCCAAAAtgc                                                         >  2:2430490/1‑93 (MQ=255)
                                      gtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt                                                        >  2:3828195/2‑91 (MQ=255)
                                       tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAatgatg                                                       >  1:399076/1‑93 (MQ=255)
                                       tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                       <  2:1910281/93‑1 (MQ=255)
                                       tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGCCAACGTTGTTAAGCCTGCCCGTGTTGGCTATTCAAAAAtgttg                                                       >  1:5225574/1‑93 (MQ=255)
                                        ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                      <  2:1619890/93‑1 (MQ=255)
                                        ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                      <  2:744819/93‑1 (MQ=255)
                                        ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                                      <  2:5858789/93‑1 (MQ=255)
                                         tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                                     >  1:4266941/1‑93 (MQ=255)
                                         tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTATGCCCGTGTTGGCTATTCAAAAATgttgct                                                     >  1:5893127/1‑93 (MQ=255)
                                                 cgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGttt                                             <  2:419430/91‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  2:5792067/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  2:2770265/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  2:1003599/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  1:2320615/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  2:1339602/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  2:4071916/92‑1 (MQ=255)
                                                  gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                            <  1:1275846/92‑1 (MQ=255)
                                                   atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                           >  2:2898803/1‑93 (MQ=255)
                                                    tcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                                          <  1:4941477/93‑1 (MQ=255)
                                                    tcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTTTTTATg                                          >  2:5998980/1‑93 (MQ=255)
                                                     catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGc                                         <  2:3401540/93‑1 (MQ=255)
                                                     catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGATGCTGGGTGTTTATGc                                         >  2:5493171/1‑93 (MQ=255)
                                                      atGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcg                                        >  1:3935126/1‑92 (MQ=255)
                                                        gAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTa                                      <  1:1845718/93‑1 (MQ=255)
                                                         aaTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTat                                     <  2:2703552/93‑2 (MQ=255)
                                                         aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                     <  2:535964/93‑1 (MQ=255)
                                                          aTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTatt                                    <  1:5060412/93‑3 (MQ=255)
                                                          aTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTatt                                    <  1:5165832/93‑3 (MQ=255)
                                                          aTGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTatt                                    <  2:3048347/93‑3 (MQ=255)
                                                          aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                    >  1:6068609/1‑93 (MQ=255)
                                                          aTGCGATGGGCAGCGACTACATCCGTGAGCTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                    <  2:3131321/93‑1 (MQ=255)
                                                             cGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTta                                 >  1:2636891/1‑93 (MQ=255)
                                                              gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                                >  2:1839642/1‑93 (MQ=255)
                                                                tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                              <  1:3737322/93‑3 (MQ=255)
                                                                  ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATagag                            <  1:2951307/93‑1 (MQ=255)
                                                                    taacGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCa                          <  2:5132580/90‑1 (MQ=255)
                                                                    taacGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCa                          >  1:4210676/4‑93 (MQ=255)
                                                                    cAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCa                          <  2:5602137/93‑1 (MQ=255)
                                                                      gCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa                        <  2:5014904/93‑1 (MQ=255)
                                                                       cGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACaa                       <  2:1690561/93‑1 (MQ=255)
                                                                        gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAag                      >  2:4293142/4‑93 (MQ=255)
                                                                         attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc                     <  2:2584158/91‑1 (MQ=255)
                                                                         attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc                     <  1:839402/91‑1 (MQ=255)
                                                                         aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                     <  1:1968964/93‑1 (MQ=255)
                                                                         aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                     >  2:779360/1‑93 (MQ=255)
                                                                          ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca                    >  2:6089575/2‑93 (MQ=255)
                                                                          cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca                    >  2:3586016/1‑93 (MQ=255)
                                                                           tACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcag                   >  2:5670425/1‑93 (MQ=255)
                                                                            aCATCCGTGAGGTGAATGTGGTGCAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc                  >  1:1717871/1‑93 (MQ=255)
                                                                            aCATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc                  <  2:3530372/93‑1 (MQ=255)
                                                                            aCATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc                  >  1:2106271/1‑93 (MQ=255)
                                                                             cATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcg                 >  2:1801179/1‑93 (MQ=255)
                                                                              aTCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGATGTTTATGCCTACTTTATAGAGCATAAGCAgcgc                <  1:457389/93‑1 (MQ=255)
                                                                               tacGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcgaa               >  1:2858451/3‑91 (MQ=255)
                                                                               tCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCa               <  2:2067145/93‑1 (MQ=255)
                                                                                ccGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa              <  2:373426/93‑1 (MQ=255)
                                                                                ccGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa              <  1:5968523/93‑1 (MQ=255)
                                                                                acGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaa              <  1:3401499/92‑1 (MQ=255)
                                                                                  gTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaatt            <  2:2272218/93‑3 (MQ=255)
                                                                                  gTGAGGCCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAGCGaaatt            <  1:2230616/93‑3 (MQ=255)
                                                                                      ggggAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt        >  1:4498280/4‑93 (MQ=255)
                                                                                         gAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAt     <  2:324512/93‑1 (MQ=255)
                                                                                            |                                                                                            
GAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCC  >  NC_000913/580703‑580887

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