Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000961,708,7090TA21.8% 102.8 / 28.2 89R597R (CGT→CGA) rsxCSoxR [2Fe‑2S] reducing system protein RsxC
Reads supporting (aligned to +/- strand):  ref base T (23/38);  new base A (9/8);  total (32/57)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGA  >  U00096/1708561‑1708852
                                                                                                                                                    |                                                                                                                                               
cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                                                                                                 <  1:1639167/149‑1 (MQ=255)
     aCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                                                                                                            <  1:1639170/149‑1 (MQ=255)
      ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                                                                                                           <  2:1639171/149‑1 (MQ=255)
        gTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgc                                                                                                                                         <  2:3347832/149‑1 (MQ=255)
         ttGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgc                                                                                                                                         <  1:1639172/148‑1 (MQ=255)
            cTGATCCACGTAATACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                                                                     <  2:1639173/149‑1 (MQ=255)
               aTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                                                                                                  <  1:1639174/149‑1 (MQ=255)
                  cACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                                                                               <  2:1639175/149‑1 (MQ=255)
                    cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                                                                             <  2:1639177/149‑1 (MQ=255)
                       aaaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACCACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACa                                                                                                                          <  2:1639178/149‑1 (MQ=255)
                         aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGCACAGCAACAGGCTAATGCGGAACCAGAACACCAGGCCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                        <  1:1639179/149‑1 (MQ=255)
                         aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                                                                        <  1:1639180/149‑1 (MQ=255)
                          aCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                                                                                       <  1:1639181/149‑1 (MQ=255)
                            tGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACACCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                                                                                                     <  1:1639182/149‑1 (MQ=255)
                             gCCGTTGACGCAGCTATCGCCCGCTCCAAAGCGCGCAAGCTTGAACAGCAACAGGCTTATGCGGAACCAGACCAACCGGTCGAGCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                    <  1:1639186/149‑1 (MQ=14)
                             gCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                    <  1:1639187/149‑1 (MQ=255)
                               agTTGCAGCAGCTCGCTCACGCGCCAAATCTCGCAAGCTGGAACAGCAACAGGCTACTCCGGACCCAGAACACCTGGCCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAgg                                                                                                                  <  2:1639197/148‑1 (MQ=14)
                                gTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGc                                                                                                                 <  1:1639192/149‑1 (MQ=255)
                                   gAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                                                              <  2:1639193/149‑1 (MQ=255)
                                               tccgtcgcCAATGCTCGCAAGCGTTAACAGCAACTGTCTAATGCGTATCTCTTACATCAGGGCGATCCGGGAAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga                                                                                                  <  2:4025257/144‑1 (MQ=9)
                                                    cgcCAAAGCGCGCAAGCTGGAACAGCAACAGCCAAATGCGGAACCAGCACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAc                                                                                             <  2:1639199/149‑1 (MQ=14)
                                                     gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                            <  2:1639198/101‑1 (MQ=12)
                                                     gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAGGCAGCCGACGAGCCCGATATTGCCAGTGACTAAGCGCGCAAGCTGGCCCAGCATCAGGCTAAGGCGCAACCAGACGAGCa                                                                                            >  2:4036157/1‑149 (MQ=14)
                                                     gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                            >  1:1639198/1‑101 (MQ=17)
                                                     gcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGAAACCAGAACAACAGCTCGATCCGAAAAAAGCAACCTTCGAAGCCGCTATTGAACGTGAAAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCg                                                                                                          >  2:1639367/1‑135 (MQ=11)
                                                        aaaGCGCTCTAGCTTGATCAGCAACAGGCTAATTCGGAACCAGAACAACAGGCCGATCCTCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                         <  2:1639200/149‑1 (MQ=14)
                                                        aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTAGATCCGCGAAAGGCAGCCGTCGACGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                         >  2:1639295/1‑149 (MQ=11)
                                                         acgctcgcAAGCTGGCACAGCAACAGGCAACTGCGGACCCAGCACCCCAGGCCGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                        <  2:1639203/144‑1 (MQ=14)
                                                         aatctcgcAAGCTGGAACATCATCAGCATAATTCGGAACCAGAACAACAGGCCGATCCTCGCAAATCTGCCGTCGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                        <  2:1639205/144‑1 (MQ=14)
                                                              cggaaGCTGGAAGACCAACAGCAAAATCAGGACGCAGAACCAGAGGACGATCCGCGCAAATCCGCCGTCGAAGCCGCTATGGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATcc                                                                                   <  2:1639215/146‑1 (MQ=14)
                                                                   gCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                            >  1:1639298/1‑87 (MQ=18)
                                                                   ccTGGAACACCTTCAGGAAAATCCGGACGCAGAACACGAGGCCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCa                                                                              <  2:1639218/148‑1 (MQ=14)
                                                                     tGGAACCGCAACATGCTACTGCGCCACCAGACGACCAGGCCGATCCTCGCACAGCCGCCGCCGCCGCCGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa                                                                            <  2:4685744/149‑1 (MQ=11)
                                                                     tGGAACAGCAACAGGCTAATGCGGTACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                                                                                         <  1:1639301/88‑1 (MQ=255)
                                                                     tGGAACAGCAACAGGCTAATGCGGTACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                                                                                         >  2:1639301/1‑88 (MQ=255)
                                                                      ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                                                                                                 >  2:1639251/1‑79 (MQ=34)
                                                                       gAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGGAGAACAGGTCGATCCGCGCAAAgc                                                                          >  1:3745731/1‑149 (MQ=255)
                                                                        aaCATCAACAGGCTAATGCGGAACCAGAACCACAGGCCGATCCCCGCAAACCCGCCGTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                         <  2:1639209/149‑1 (MQ=14)
                                                                        aaCACCAACAGCCTACTCCGGACCCAGAACAACAGGCCGATCCGCGCAAATCAGCCGTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgcc                                                                         <  2:1639208/149‑1 (MQ=14)
                                                                              aaCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                            <  1:1639303/76‑1 (MQ=21)
                                                                              aaCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                                                                            >  2:1639303/1‑76 (MQ=21)
                                                                               aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAg                                                                                                                                          <  1:1639305/77‑1 (MQ=21)
                                                                               aCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAg                                                                                                                                          >  2:1639305/1‑77 (MQ=21)
                                                                                 aGGCTAATGCGCAACCAGAACACCCGGTCGATCCGCGCAAAGCCGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                <  2:2814873/149‑1 (MQ=255)
                                                                                  ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGcc                                                                                                                                              >  2:1639311/1‑70 (MQ=25)
                                                                                  ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCTCGCAAGCTgg                                                                                                                              >  2:1639386/1‑86 (MQ=11)
                                                                                      aaTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                                           >  1:1639313/1‑69 (MQ=25)
                                                                                      aaTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                                           <  2:1639313/69‑1 (MQ=25)
                                                                                       aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                                                                       <  1:1639308/72‑1 (MQ=21)
                                                                                       aTGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                                                                       >  2:1639308/1‑72 (MQ=21)
                                                                                               ccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                                                                                                 >  1:4184744/1‑54 (MQ=31)
                                                                                                  gAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGc                                                                                                           >  2:1639314/1‑89 (MQ=255)
                                                                                                      aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                           <  2:1639212/149‑1 (MQ=255)
                                                                                                         aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                         >  1:1639213/1‑100 (MQ=37)
                                                                                                            ccGATCCCCGCAAACCTGCCGTCGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGc                                     <  2:1639214/148‑1 (MQ=11)
                                                                                                              gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                                                                                              <  1:1639220/74‑1 (MQ=34)
                                                                                                              gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACa                                                                                                                    <  1:1639219/68‑1 (MQ=25)
                                                                                                                tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTaa                                                                                                              >  2:1639319/1‑72 (MQ=21)
                                                                                                                tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGATGGGCAGGTCGATCCGCGCTACGCCGCGGTCGGAGCCGCTATTGCCCGTGCCAACGCGCGCAAGc                                 >  2:1639221/1‑149 (MQ=14)
                                                                                                                       aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATg                                                                                                            >  1:1901684/1‑67 (MQ=35)
                                                                                                                        aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                        >  2:1639332/1‑85 (MQ=255)
                                                                                                                             gccgTTGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                          <  1:1639224/79‑1 (MQ=255)
                                                                                                                             gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATg                                                                                                            <  1:1639228/61‑1 (MQ=35)
                                                                                                                             gccgTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAgg                                                                                          <  1:1639227/79‑1 (MQ=255)
                                                                                                                             gccgTCGAAGCCGCTATCGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                                                                                                        <  1:1639232/65‑1 (MQ=255)
                                                                                                                             gccgTCGAAGCCGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                                                                                                        <  1:1639230/65‑1 (MQ=255)
                                                                                                                             gccgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                                                                        <  1:1639225/81‑1 (MQ=255)
                                                                                                                               cgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                                                                                       <  1:2589891/48‑1 (MQ=38)
                                                                                                                               cgTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                                                                                                               <  1:1639235/56‑1 (MQ=255)
                                                                                                                                gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                               >  2:1639234/1‑119 (MQ=255)
                                                                                                                                gTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                               <  1:1639234/119‑1 (MQ=255)
                                                                                                                                 tggatgCCGCTATTGCCCGCGCCATAGCGCGCAAGCTGGAACATCATCAGGATAATCCCGAACCAGAACAACAGGGCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt                <  1:1639244/144‑1 (MQ=11)
                                                                                                                                 tCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACa                                                                                            >  1:1639341/1‑73 (MQ=255)
                                                                                                                                   gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcgcg                                                                                <  2:1639421/83‑1 (MQ=255)
                                                                                                                                   gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCcgcg                                                                                >  1:1639421/1‑83 (MQ=255)
                                                                                                                                   gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCg                                                                                       <  1:1639240/76‑1 (MQ=255)
                                                                                                                                   gAAGCCGCTATAGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGc                                                                                                                 <  1:1639241/50‑1 (MQ=11)
                                                                                                                                    aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCg                                                                                                          >  1:1639344/1‑56 (MQ=255)
                                                                                                                                    aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCa                                                                              >  2:1639343/1‑84 (MQ=255)
                                                                                                                                    aaGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGt                                                                                         <  1:1639242/73‑1 (MQ=255)
                                                                                                                                      tttgCTATAGCCCGTGCCAAAGCGCGAAAGCTGGAACAGCAACAGGCTAATTCGGAACCAGAAGAACAGTTCTATCCGCTCAAAGCCTCCGCCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGc           <  2:1639246/146‑1 (MQ=11)
                                                                                                                                      gCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAgc                                                                          >  1:1639433/1‑86 (MQ=255)
                                                                                                                                           tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                         <  2:1639354/66‑1 (MQ=255)
                                                                                                                                           tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                         >  1:1639354/1‑66 (MQ=255)
                                                                                                                                            atcgCCTGTGCCAAAGCGCTCATCCAGGACCACCAACACCCCAATGCGGACCCATTCGCACCTCTCGTTCCGCGCCCAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGcc     <  2:4713772/146‑1 (MQ=9)
                                                                                                                                            aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCg                                                                                                          >  1:1639358/1‑48 (MQ=255)
                                                                                                                                            aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                                                                                                        >  1:1639357/1‑50 (MQ=255)
                                                                                                                                            aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCaga                                                                                                  >  1:1639356/1‑56 (MQ=255)
                                                                                                                                            aTTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCAGCCGt                                                                    >  2:1639355/1‑86 (MQ=255)
                                                                                                                                               ggccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa  <  1:4539726/147‑1 (MQ=255)
                                                                                                                                                cccGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGa                                                                                                        >  1:1639359/1‑46 (MQ=38)
                                                                                                                                                  cGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCaaa                                                                            >  1:1639286/1‑72 (MQ=39)
                                                                                                                                                    tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAgg                                                                                          >  1:1639316/1‑56 (MQ=255)
                                                                                                                                                    tGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTAGATCCGCGCAAAGCTGCCGt                                                                    >  2:1639362/1‑78 (MQ=11)
                                                                                                                                                    |                                                                                                                                               
CAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGA  >  U00096/1708561‑1708852

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: ‑