Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,7830CA48.2% 405.0 / inf 1116I69I (ATC→ATAybcXhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (324/253);  new base A (263/274);  total (588/528)
Fisher's exact test for biased strand distribution p-value = 1.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1173 total aligned reads displayed.

TACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAC  >  NC_000913/580688‑580878
                                                                                               |                                                                                               
  cTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCg                                                                                               >  2:1520015/1‑96 (MQ=255)
  cTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACg                                                                                               >  1:8215698/1‑96 (MQ=255)
   tATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGt                                                                                              >  1:8960014/1‑96 (MQ=255)
      ctcCCCAAAGCGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAg                                                                                           <  1:2300645/96‑1 (MQ=255)
        ccccAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGt                                                                                         <  1:3480076/96‑1 (MQ=255)
        ccccAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGt                                                                                         <  2:2145494/96‑1 (MQ=255)
        cccGAAAAAATCCGCATACCAGGAAGGGCGCGGGGAACCCCTGCCCTTCCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGt                                                                                         <  1:5385461/96‑1 (MQ=255)
          ccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                       <  2:5338145/94‑3 (MQ=255)
          ccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                       >  1:3872567/3‑94 (MQ=255)
           caaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaa                                                                                      >  2:8536943/2‑96 (MQ=255)
            aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                                     <  1:10001121/96‑1 (MQ=255)
            aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                                     <  2:9606221/96‑1 (MQ=255)
             aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaacg                                                                                    >  2:11327340/1‑94 (MQ=255)
             aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaacg                                                                                    >  1:6617076/1‑94 (MQ=255)
              agagTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAAcgt                                                                                   <  1:9503158/96‑1 (MQ=255)
              aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtgt                                                                                   <  2:11950217/96‑1 (MQ=255)
                  tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCTTTTCAGCGGGCCATCAGGAATGCGATGGGCTGCGACGACATCCGTGCGGGGAAtgtgagg                                                                               >  1:2993419/1‑93 (MQ=255)
                     gcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                            <  2:3393359/93‑1 (MQ=255)
                     gCATACCAGGGAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGGGAGGTGAATGGGGTGCAg                                                                            >  1:5987946/1‑96 (MQ=255)
                     gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                            >  2:9026390/1‑96 (MQ=255)
                     gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                            >  1:3200851/1‑96 (MQ=255)
                      cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                                           >  2:6547044/1‑96 (MQ=255)
                      cATACCAGGAAGGCGGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTTCATCCGTGAGGTGAATGTGGTGAAGt                                                                           <  1:6946119/96‑1 (MQ=255)
                        tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTGGATAAGTCt                                                                         >  1:1647335/1‑96 (MQ=255)
                         aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                        <  2:13083703/96‑1 (MQ=255)
                         aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                        <  2:3983769/96‑1 (MQ=255)
                          ccAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                                                       >  2:6961916/1‑96 (MQ=255)
                            aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGccc                                                                     <  2:976547/96‑1 (MQ=255)
                            aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcac                                                                     <  1:11517192/96‑3 (MQ=255)
                             ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                    >  2:8048284/1‑96 (MQ=255)
                              gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                                   >  1:12116568/1‑96 (MQ=255)
                              gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                                   <  1:12395971/96‑1 (MQ=255)
                              gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                                   <  2:12266103/96‑1 (MQ=255)
                               aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtg                                                                  <  1:11223244/96‑1 (MQ=255)
                                 gggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgtt                                                                <  1:5649533/96‑2 (MQ=255)
                                  ggCGCTGGAAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                               <  2:10750777/96‑1 (MQ=255)
                                   gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATCCGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttgn                                                              <  1:5157700/96‑4 (MQ=255)
                                   gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                                                              <  1:7611826/96‑1 (MQ=255)
                                     gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                            >  2:6164220/1‑96 (MQ=255)
                                     gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                            >  1:2512224/1‑96 (MQ=255)
                                      cagggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                           >  1:10834577/3‑96 (MQ=255)
                                      cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTa                                                           >  2:13527785/1‑96 (MQ=255)
                                       tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAt                                                          >  1:573345/1‑96 (MQ=255)
                                        gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAtt                                                         <  1:7877597/96‑1 (MQ=255)
                                        gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTAtt                                                         >  2:2001965/1‑96 (MQ=255)
                                         ggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTc                                                        >  2:9379685/1‑96 (MQ=255)
                                          gAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGGAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTca                                                       >  2:7745046/1‑95 (MQ=255)
                                             acgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCaaaa                                                    <  2:2561039/92‑1 (MQ=255)
                                             acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                                    >  2:12666211/1‑96 (MQ=255)
                                              cgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                   >  2:7290943/4‑96 (MQ=255)
                                               acTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                                                  >  1:8601813/1‑96 (MQ=255)
                                                  gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt                                               <  2:829520/96‑3 (MQ=255)
                                                  gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgtt                                               <  1:1359279/96‑3 (MQ=255)
                                                  gCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                                               >  2:133227/1‑96 (MQ=255)
                                                    cgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTACTCAAAAAtgttgc                                             <  1:13894209/94‑1 (MQ=255)
                                                    ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGATATTCCAAAAtgctgc                                             >  1:12800362/1‑96 (MQ=255)
                                                     cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                            <  1:2736507/96‑1 (MQ=255)
                                                     cTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                            >  1:13350285/1‑96 (MQ=255)
                                                       ttCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTgg                                          >  2:14012127/1‑96 (MQ=255)
                                                       ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTgg                                          >  1:711327/1‑96 (MQ=255)
                                                        tCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg                                         <  2:14161879/96‑1 (MQ=255)
                                                          aGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGCGGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtg                                       >  1:13018370/1‑96 (MQ=255)
                                                            cgagCCATCATGAATGCGATGGGCAGCTACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGATGCTGGGTGtt                                     <  2:4881663/93‑1 (MQ=255)
                                                            cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                                     <  1:4530978/96‑1 (MQ=255)
                                                            cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGCGAGGTGGATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                                     >  2:2368764/1‑96 (MQ=255)
                                                              ggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                   <  2:8857618/96‑1 (MQ=255)
                                                               gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                                  <  2:13731132/93‑1 (MQ=255)
                                                                cgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                                 <  1:1903695/94‑1 (MQ=255)
                                                                cgatcatGAATGCGATGGGTAACGATTACAGACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                                 <  1:10293032/94‑1 (MQ=255)
                                                                    catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTa                             <  2:10461273/96‑1 (MQ=255)
                                                                     atGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTat                            >  1:3694594/1‑95 (MQ=255)
                                                                     atGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTat                            >  1:9140795/1‑95 (MQ=255)
                                                                       gAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTAttt                          <  1:4750991/96‑1 (MQ=255)
                                                                        aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACttt                         <  2:5838167/96‑1 (MQ=255)
                                                                         aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                        <  1:2114074/96‑1 (MQ=255)
                                                                           gCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtatt                      >  2:4470787/1‑95 (MQ=255)
                                                                           gCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                      <  2:2839462/96‑1 (MQ=255)
                                                                             gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga                    >  1:14157406/1‑96 (MQ=255)
                                                                              aTGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattcag                   <  2:1709041/96‑5 (MQ=255)
                                                                                 ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                >  2:9179787/1‑96 (MQ=255)
                                                                                 ggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGGCGGTTATTCCAAAATGCTGCTGGGTGTTTATGGCTACTTTATAGAGCAt                >  1:7332448/1‑96 (MQ=255)
                                                                                  gCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa               >  2:11566962/1‑96 (MQ=255)
                                                                                  gCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa               <  2:6318338/96‑1 (MQ=255)
                                                                                  gCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa               <  2:3768090/96‑1 (MQ=255)
                                                                                   taacGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACaa              <  2:4253265/93‑1 (MQ=255)
                                                                                   taacGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACaa              >  2:1108399/4‑96 (MQ=255)
                                                                                    aGCGACTACATCCGTGAGGCGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAag             >  1:8423376/1‑96 (MQ=255)
                                                                                     acGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc            <  2:10250141/95‑1 (MQ=255)
                                                                                      cGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca           >  1:6242636/1‑96 (MQ=255)
                                                                                      cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca           >  1:4225619/1‑96 (MQ=255)
                                                                                      cGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAacca           >  1:11871784/1‑96 (MQ=255)
                                                                                       gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcag          >  1:14412865/4‑96 (MQ=255)
                                                                                       gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcag          >  2:3168319/4‑96 (MQ=255)
                                                                                        aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc         <  1:9766239/96‑1 (MQ=255)
                                                                                        aCTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcaac         <  1:10577674/96‑3 (MQ=255)
                                                                                         tgaCATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAAGCAgcg        >  1:9367563/3‑96 (MQ=255)
                                                                                         cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg        >  1:4971093/1‑96 (MQ=255)
                                                                                         cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg        >  2:3035088/1‑96 (MQ=255)
                                                                                         cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcg        >  2:5552018/1‑96 (MQ=255)
                                                                                          tACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcgc       <  2:7507918/96‑1 (MQ=255)
                                                                                               |                                                                                               
TACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAC  >  NC_000913/580688‑580878

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