Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000961,708,9300CA23.2% 175.9 / 20.8 82P671Q (CCG→CAG) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base C (43/20);  new base A (11/8);  total (54/28)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTG  >  U00096/1708792‑1709076
                                                                                                                                          |                                                                                                                                                  
ccGCTCTAGCCCGTGCCACAGCCCGCAAGCTGGAACAGCAACAGTATAATGCAGAGCCAGAACATCTGCTCGATCGGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCg                                                                                                                                          <  1:1268898/149‑1 (MQ=14)
  gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCAAAAGCCCGCAAACGCGAACAGCAACAGGCTAATGCGGa                                                                                                                                        >  2:1269030/1‑149 (MQ=14)
    tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCAAAAGCCCGCAAACGGGAACAGCACCAGGCTAATGCGGAGc                                                                                                                                      >  2:1269022/1‑149 (MQ=14)
        gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACACCAGGTCGATCCGCGCATAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                                                                                                                                  <  2:1996974/149‑1 (MQ=255)
          ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCACAAGCGCGCAAACGGGAACAGCAGCCGGCTAATGCGGAGCCagaag                                                                                                                                >  2:1269020/1‑149 (MQ=255)
                      aCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCGAAGGCCCGAAAACCGGAACAGCAGCAGGCTATAGCGGAGCCAGAAGACGAGGTTGCTc                                                                                                                    >  2:1269028/1‑149 (MQ=14)
                                     gCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGCCCAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGcag                                                                                                     >  2:1269048/1‑147 (MQ=255)
                                           ggCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGCTATAGCACGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGTGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGtggccg                                                                                               >  2:1268972/1‑144 (MQ=17)
                                            gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGc                                                                                                                                           <  1:1268916/104‑1 (MQ=255)
                                            gCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGc                                                                                                                                           >  2:1268916/1‑104 (MQ=255)
                                             cTAATGCCGAGCCAGAACAACAGGTTGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGCGCACAAGCCCGCAAACGCGACCAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGAAAAGCTGCAGtggcaccg                                                                                             >  2:1269051/1‑142 (MQ=14)
                                                           gAACAACAGGTCGAACCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCAAGAGCCAAAGCGCGAAATCTGGAGCAGCAAAATGCTAATGCGGAGCCAGAAGAACATGTAGATAAGCGCAAAGCGGCCGTCGAAGCGGCTATTGCaggagc                                                                               >  2:1268981/1‑144 (MQ=11)
                                                             agaacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCACGCAAACGGGAACAGCAACAGGCTACTGCAGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCGGCCGTGGCCGCGGCTATTGCACgagcga                                                                             >  2:1269059/3‑147 (MQ=11)
                                                               aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAACGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGACCAGGTTGATCCGCGCAAAGCTGCAGTAGACGCGGCTATTgcgcgcgcacag                                                                           >  2:1269018/1‑145 (MQ=14)
                                                                 caGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTGGCCGCGGCTATTGCCCgcgcgcaggg                                                                         >  2:1269058/1‑143 (MQ=14)
                                                                  aGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCAAAAGCCAGCAAGCAGGAACAGCAACAGGCTAATGAGGATCCAGAAGAACAGGTGGATCCGCGCAAAGCGGCAGTGGAAGAGGATATTGAACGCGCGAAGGcc                                                                        >  2:1268995/1‑148 (MQ=11)
                                                                       gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCgg                                                                                                                                         <  1:1268920/79‑1 (MQ=255)
                                                                       gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCgg                                                                                                                                         >  2:1268920/1‑79 (MQ=255)
                                                                         tCCGCGCAACGCCGCCGTCGAAGCTTCTATTTCTGTTGCCTCTCCCCGCTATCTGTAACAGCTACCGGATAATGACGAGCTAGATCTATAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGc                                                                 <  2:1268985/149‑1 (MQ=255)
                                                                           cgcgCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGTCCCAGAAGAACAGGTTGATCCGCGCAAAg                                                                                                           >  1:1268990/1‑105 (MQ=255)
                                                                              gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCACCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGc                                                                                                          >  1:1269004/1‑103 (MQ=255)
                                                                                 aaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAACGCACAAAAGCAGGAACa                                                          >  1:1269011/1‑148 (MQ=17)
                                                                                 aaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGACGCGGCTATTGCACGCGCCAACGCACTCAAGCTGGACCAg                                                         >  2:1269086/1‑149 (MQ=17)
                                                                                 aaGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCAGGAACAg                                                         >  2:1269015/1‑149 (MQ=255)
                                                                                      gccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACa                                                    >  2:1269067/1‑149 (MQ=255)
                                                                                      gccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGAAGCGGCTATTGCACGCGCCAACGCCCGCAACCAGGACCAGCATCa                                                    >  2:1269083/1‑149 (MQ=17)
                                                                                       ccgTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCACGCAAGCAGGAAcagcagcag                                                   >  2:1268999/1‑149 (MQ=255)
                                                                                          tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAAt                                                                                                                                             <  1:1268931/56‑1 (MQ=255)
                                                                                          tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                         >  1:1269003/1‑76 (MQ=255)
                                                                                          tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGACCAGCAGCAGGCt                                                >  1:1269029/1‑149 (MQ=255)
                                                                                          tCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCACGCGCAAAAGCCCGCAAGCTGGAACAGCAACAGGCt                                                >  1:1269089/1‑149 (MQ=255)
                                                                                           cGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCa                                                       >  1:1269063/1‑141 (MQ=255)
                                                                                            gAAGCAGCTATCGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGc                                                                                                          <  1:1268932/89‑1 (MQ=18)
                                                                                                ccGCTATTGCCCGTGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc                                                                                                                             <  2:1268937/66‑1 (MQ=255)
                                                                                                 cGCTATTGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                                                                                                                                  <  1:1268936/60‑1 (MQ=255)
                                                                                                  gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCgg                                                                                                                                         <  1:1268940/52‑1 (MQ=255)
                                                                                                  gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaaa                                                                                                            >  1:1269057/1‑81 (MQ=255)
                                                                                                  gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGAAGCGGCTATTGCACGCGCAAAAGCACTCAAGCTGGAACAGCAGCAGGCTAATGCgga                                        >  2:1269049/1‑148 (MQ=255)
                                                                                                   cTATTGCCCGCGCCAAAGCTCGCAATCTTTAACAGCAACCGGCTAATGCCGAGCCAGAAGATCAGGTTTATCCTCGCTAAGCCGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTa                                       <  2:1268945/149‑1 (MQ=14)
                                                                                                   cTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTaa                                                                                                                                              <  2:1268897/46‑1 (MQ=255)
                                                                                                    tATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCAAAAGCACGAAAGCAGGAACAGCAACAGGATAAtcaggaag                                      >  2:1269037/1‑142 (MQ=25)
                                                                                                    tATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGCAGCGGCTATTGCACGCGCCAAGGCACAAAAGCTGGCACAGCAACCGGATATTGAGGTAc                                      >  2:1269060/1‑149 (MQ=21)
                                                                                                     aTTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTAcc                                     >  2:1430156/1‑149 (MQ=255)
                                                                                                     aTTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGCCGCGGCTATTGCACgcgc                                                                               >  1:1268987/1‑107 (MQ=255)
                                                                                                       tGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAc                                                                                   <  2:1268947/101‑1 (MQ=11)
                                                                                                       tGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGc                                                                                                                                           <  1:1268948/45‑1 (MQ=255)
                                                                                                       tGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACgcgc                                                                               >  1:1269019/1‑105 (MQ=255)
                                                                                                        gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                                                                                                                 <  1:1268949/38‑1 (MQ=255)
                                                                                                         cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTAGCCGCGGCTATTGCCCGCGCGAAAGCACAAAAGCACGACCAGCACCAGGCTTAATCGGTACCagaa                                 >  2:1269017/1‑149 (MQ=17)
                                                                                                          ccGCGCCAAAGCGCGCAATCTGGAACAGCAACCGGCTAATGCGGa                                                                                                                                        <  2:1268950/45‑1 (MQ=255)
                                                                                                           cGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCAGCAGGCTAATGCGGTACCaga                                  >  2:1269088/1‑146 (MQ=255)
                                                                                                            gAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAgaagaa                              >  2:1269045/1‑149 (MQ=255)
                                                                                                             aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGc                                                                                     >  1:1269035/1‑93 (MQ=255)
                                                                                                               ccAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGc                                                                                     >  1:1269000/1‑91 (MQ=255)
                                                                                                                 aaaGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGt                         >  2:1269079/1‑149 (MQ=255)
                                                                                                                   aGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCgcgc                                                                               >  1:1269087/1‑93 (MQ=255)
                                                                                                                      gcgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATc                    >  2:1269005/1‑149 (MQ=255)
                                                                                                                         cAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACgcgc                                                                               >  1:1269006/1‑87 (MQ=255)
                                                                                                                         cAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCagc                 >  2:1269097/1‑149 (MQ=255)
                                                                                                                          aaaCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGATAATGAGGTACCAGAAGTACAGGTAGATCcgcg                >  2:1269075/1‑149 (MQ=18)
                                                                                                                           agctggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAc                                                                                   <  2:1268958/77‑1 (MQ=11)
                                                                                                                            tctggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAg                                                                                               <  1:1268965/65‑1 (MQ=255)
                                                                                                                              gggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGt                                                                                                    >  1:1269109/1‑61 (MQ=255)
                                                                                                                               ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAg                                                                                                                           <  1:1268957/37‑1 (MQ=255)
                                                                                                                               ggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGa                                                                                                                      <  1:1268962/42‑1 (MQ=255)
                                                                                                                               ggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGcc                                                                                                      <  1:1268956/58‑1 (MQ=255)
                                                                                                                               ggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGc                                                                                              <  1:1268963/66‑1 (MQ=255)
                                                                                                                               ggAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTg                                                               <  1:1268955/97‑1 (MQ=255)
                                                                                                                                gAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTAGAAAACGCTATAGCGACATCGAAAGAAAGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGc          >  2:1269038/1‑149 (MQ=11)
                                                                                                                                 aaCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACgcgc                                                                               >  1:1269100/1‑79 (MQ=255)
                                                                                                                                  accgCAACCGGCCAATGCCGAGCCAGAACAACATGTTGATCCTCTCAAAGCGGCCGTCGAATCCGCTAATCCACGCGCCAAATCTCCCAACCTGGTACAGCAACAGGCTAATGCGGTTCCAGAAGAACAGGTTGATCCGCGCAAAGCgg        <  2:1268966/146‑1 (MQ=14)
                                                                                                                                  aCAGCAACCGGCTAATGCGGATCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACgcgc                                                                               >  1:1269062/1‑78 (MQ=255)
                                                                                                                                  aCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATc                                                                                                                    <  2:1268907/41‑1 (MQ=255)
                                                                                                                                  aCAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGACGCGGCTATTGCACgcgc                                                                               >  1:1269071/1‑78 (MQ=255)
                                                                                                                                   cAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGAAGCGGCTATTGCACGCGCCAAAGCCCGAAAGCAGGAACAGCAACAGGCTAATGAGGTACCAGAAGAACAGGTTGATCCTCTCAAAGAGGc       >  2:1269046/1‑149 (MQ=17)
                                                                                                                                   cAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                                                        <  1:1268967/36‑1 (MQ=255)
                                                                                                                                   cAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                                                        >  2:1268967/1‑36 (MQ=255)
                                                                                                                                    aGCAACCGGCTAATGCGGAGCCAGAAGCACAGGTTGATGCGCGCAAAGCTGCCGTCGATGCGGCTATTGCACGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCa      <  2:1879705/149‑1 (MQ=17)
                                                                                                                                    aGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                                                                             <  1:2763510/46‑1 (MQ=255)
                                                                                                                                       aaCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAACGCTGCCGTCGAAGCGGCTATTGCACGCGCCAACGCACGCAAGCGGGACCAGCAACAGGCTAATTCGGTACCAGAAGAACAGGTGGCTCCGCGCAGAGCGGCAGtt   >  2:1269073/1‑149 (MQ=38)
                                                                                                                                       aaCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGtt   >  1:1269002/1‑149 (MQ=255)
                                                                                                                                        aCCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCGGGAACAGCAACGGGCTACTGCGGTACCAGAAGAACAGGTTTATCCGCGCAAAGCGGCAGGTg  >  1:2363932/1‑149 (MQ=38)
                                                                                                                                          |                                                                                                                                                  
CCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTG  >  U00096/1708792‑1709076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑