Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,136,368 G→T 14.0% R386S (CGT→AGT)  hybB ← putative hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,136,3680GT14.0% 919.5 / 89.7 292R386S (CGT→AGT) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (113/138);  new base T (23/18);  total (136/156)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 295 total aligned reads displayed.

TCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAT  >  CP009273/3136221‑3136516
                                                                                                                                                   |                                                                                                                                                    
  ttCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATATTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAtt                                                                                                                                                   >  2:693647/1‑149 (MQ=255)
  gtCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAtt                                                                                                                                                   <  2:4177631/148‑1 (MQ=255)
   tCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTa                                                                                                                                                  <  1:2289701/149‑1 (MQ=255)
    ccATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACtattat                                                                                                                                                 <  1:1575436/149‑1 (MQ=255)
             gcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgtttt                                                                                                                                        >  1:3342666/1‑149 (MQ=255)
             gcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                        <  2:1182739/149‑1 (MQ=255)
              ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTATCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttg                                                                                                                                       >  2:2301684/1‑149 (MQ=255)
                acaTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgtt                                                                                                                                     <  1:3972499/149‑1 (MQ=255)
                  aTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgttttgtttt                                                                                                                                   >  2:4231937/1‑149 (MQ=255)
                     gTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                                >  2:2903854/1‑149 (MQ=255)
                      tACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAag                                                                                                                               >  1:4510707/1‑149 (MQ=255)
                       aCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAagg                                                                                                                              <  1:1981863/149‑1 (MQ=255)
                         cggccgCCCACGCTTTCGAAACGACGCCATTGTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAaggag                                                                                                                            <  2:4730111/144‑1 (MQ=255)
                           gAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAggagga                                                                                                                          <  1:1709187/149‑1 (MQ=255)
                            aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGaa                                                                                                                         <  2:4055319/149‑1 (MQ=255)
                            aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGaa                                                                                                                         <  1:103657/149‑1 (MQ=255)
                            aaGCCCATGCTTTCGAAACGACGCCATTTTCGAGTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGaa                                                                                                                         <  2:3067930/149‑1 (MQ=255)
                             aGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAg                                                                                                                        <  2:734696/149‑1 (MQ=255)
                              gCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGt                                                                                                                       <  2:866932/149‑1 (MQ=255)
                              gCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGt                                                                                                                       <  1:3710328/149‑1 (MQ=255)
                                ccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                                     <  1:177459/149‑1 (MQ=255)
                                 cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATc                                                                                                                    >  1:3811148/1‑149 (MQ=255)
                                 cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATc                                                                                                                    >  1:2866058/1‑149 (MQ=255)
                                   tGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCgg                                                                                                                  >  1:3925484/1‑149 (MQ=255)
                                     cTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCa                                                                                                                <  1:550890/149‑1 (MQ=255)
                                      tttCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAg                                                                                                               >  1:3001739/1‑149 (MQ=255)
                                      tttCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAg                                                                                                               <  1:735310/149‑1 (MQ=255)
                                       ttCGACACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGt                                                                                                              <  2:4519221/149‑1 (MQ=255)
                                       ttCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGt                                                                                                              <  2:2084920/149‑1 (MQ=255)
                                        tCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTa                                                                                                             <  1:614907/149‑1 (MQ=255)
                                           aaaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGAc                                                                                                          >  2:2441302/1‑149 (MQ=255)
                                           aaaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGAc                                                                                                          <  1:4619400/149‑1 (MQ=255)
                                            aaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACg                                                                                                         >  2:2204059/1‑149 (MQ=255)
                                             acgacgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGa                                                                                                        >  2:1007479/1‑149 (MQ=255)
                                                acgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATg                                                                                                     <  1:4311305/149‑1 (MQ=255)
                                                acgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATg                                                                                                     >  2:3758003/1‑149 (MQ=255)
                                                acgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATg                                                                                                     <  2:1281681/149‑1 (MQ=255)
                                                   ccATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGt                                                                                                  >  1:2053857/1‑149 (MQ=255)
                                                    cATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTa                                                                                                 <  1:3513381/149‑1 (MQ=255)
                                                    cATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTa                                                                                                 >  1:2532291/1‑149 (MQ=255)
                                                       tttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGa                                                                                              <  2:2551242/149‑1 (MQ=255)
                                                       tttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGa                                                                                              >  1:2657201/1‑149 (MQ=255)
                                                        ttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGAt                                                                                             <  1:1029743/149‑1 (MQ=255)
                                                           gATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTa                                                                                          >  2:2016628/1‑149 (MQ=255)
                                                           gATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTa                                                                                          >  2:1910068/1‑149 (MQ=255)
                                                            aTTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTaa                                                                                         >  1:921102/1‑149 (MQ=255)
                                                             tttCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAg                                                                                        >  1:1247931/1‑149 (MQ=255)
                                                             tttCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAg                                                                                        >  2:148039/1‑149 (MQ=255)
                                                              ttCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgc                                                                                       <  2:315738/149‑1 (MQ=255)
                                                                   aaTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                                  >  2:3448064/1‑149 (MQ=255)
                                                                     tCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAtc                                                                                >  1:1409006/1‑149 (MQ=255)
                                                                       gATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAtctc                                                                              <  1:1119735/149‑1 (MQ=255)
                                                                        aTGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCa                                                                             <  2:27378/149‑1 (MQ=255)
                                                                         tgagTAAATGCCCCACAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                            <  2:4569155/146‑1 (MQ=255)
                                                                         tGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                            >  2:1318283/1‑149 (MQ=255)
                                                                           cGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATa                                                                          <  2:751576/149‑1 (MQ=255)
                                                                             tAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGc                                                                        >  1:2785822/1‑149 (MQ=255)
                                                                              aaaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                                                                       <  1:4467306/149‑1 (MQ=255)
                                                                              aaaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTACGCGCAAATCTCAATAGcc                                                                       >  2:332177/1‑149 (MQ=255)
                                                                               aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCca                                                                      >  2:3695072/1‑149 (MQ=255)
                                                                               aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCca                                                                      <  2:3068599/149‑1 (MQ=255)
                                                                                 tGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcaca                                                                    <  2:3816999/149‑1 (MQ=255)
                                                                                   ccccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACaaa                                                                  >  1:3870612/1‑149 (MQ=255)
                                                                                       tCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCa                                                              <  1:3263457/149‑1 (MQ=255)
                                                                                        cAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa                                                             <  1:2978357/149‑1 (MQ=255)
                                                                                        cAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa                                                             <  2:3919626/149‑1 (MQ=255)
                                                                                        cAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa                                                             <  1:3490901/149‑1 (MQ=255)
                                                                                         aaTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAAt                                                            >  2:2454108/1‑149 (MQ=255)
                                                                                         aaTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAAt                                                            <  2:592582/149‑1 (MQ=255)
                                                                                         aaTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAAt                                                            >  1:3166336/1‑149 (MQ=255)
                                                                                           tACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAg                                                          >  2:3846807/1‑149 (MQ=255)
                                                                                           tACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAg                                                          <  1:233234/149‑1 (MQ=255)
                                                                                            aCGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGa                                                         <  2:4688798/149‑1 (MQ=255)
                                                                                             cGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGaa                                                        <  1:4154372/149‑1 (MQ=255)
                                                                                                gTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATc                                                     >  2:4376469/1‑149 (MQ=255)
                                                                                                  tACCGGATCAATAGTAATTCTCTGGCTCATGATTTGCTCTCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCTGCAGTAGAAGAATGAGTAAGATGTAAGCGCAAATCTCAATAGCCACAAAACAAATAGAAATCaa                                                   >  2:794221/1‑149 (MQ=255)
                                                                                                    ccGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACa                                                 >  2:4781117/1‑149 (MQ=255)
                                                                                                      ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGt                                               <  2:4549542/149‑1 (MQ=255)
                                                                                                       gATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                              <  1:4607122/149‑1 (MQ=255)
                                                                                                       gATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                              >  1:1538417/1‑149 (MQ=255)
                                                                                                        aTCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttc                                             <  1:4576595/149‑1 (MQ=255)
                                                                                                         tCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttct                                            >  2:3643888/1‑149 (MQ=255)
                                                                                                            aTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcc                                         <  1:3604146/149‑1 (MQ=255)
                                                                                                              aGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                       >  2:153033/1‑149 (MQ=255)
                                                                                                               gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAg                                      >  2:1724169/1‑149 (MQ=255)
                                                                                                                  actctcTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTc                                   <  2:4608539/147‑1 (MQ=255)
                                                                                                                    tctcTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCgg                                 <  2:1365686/149‑1 (MQ=255)
                                                                                                                     ctctGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                                <  2:2871240/149‑1 (MQ=255)
                                                                                                                      tctGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGa                               >  1:2385546/1‑149 (MQ=255)
                                                                                                                      tctGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGa                               >  1:3329596/1‑149 (MQ=255)
                                                                                                                          gCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                           <  1:4730935/149‑1 (MQ=255)
                                                                                                                          gCTCAGGCATTGCGCGCCTCCCGACGCTGATGAACGCTTTGTTTTACCGGAGGAAGTATCGGCAGGAGACGCATGAGTACGATGTAAGCGCAAAGCTCAATAGCCACAAAACCAACAGACATCAACAGTTCTTCCCAGGTCGGGAAGTa                           <  2:3298528/149‑1 (MQ=255)
                                                                                                                                       tCGCCTCATGACGATTATGATCGTTTTGTTGTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAACTCAACAGTTCTTCCCAGGTCGGGACGTAGGCGTAAccgccg              >  2:2819940/1‑149 (MQ=255)
                                                                                                                                        cGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACcgccgc             <  2:2941185/149‑1 (MQ=255)
                                                                                                                                          ccTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGccc           <  2:4510623/149‑1 (MQ=255)
                                                                                                                                             cATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCACTAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCggg        <  2:1888206/149‑1 (MQ=255)
                                                                                                                                              aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCgggg       <  1:487494/149‑2 (MQ=255)
                                                                                                                                              aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt       >  1:2421317/1‑149 (MQ=255)
                                                                                                                                              aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt       >  2:499609/1‑149 (MQ=255)
                                                                                                                                                  cGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGaa   <  1:3108386/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAT  >  CP009273/3136221‑3136516

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