Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,136,368 G→T 20.9% 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,3680GT20.9% 724.0 / 134.4 262R386S (CGT→AGT) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (88/119);  new base T (31/24);  total (119/143)
Fisher's exact test for biased strand distribution p-value = 6.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 264 total aligned reads displayed.

CTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATG  >  CP009273/3136222‑3136517
                                                                                                                                                  |                                                                                                                                                     
 ttCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTa                                                                                                                                                   >  2:409688/1‑150 (MQ=255)
 ttCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTa                                                                                                                                                   >  2:2610831/1‑150 (MQ=255)
   ccATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATg                                                                                                                                                 >  2:3574005/1‑150 (MQ=255)
      tGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATc                                                                                                                                              <  1:3897639/150‑1 (MQ=255)
      tGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATc                                                                                                                                              >  1:1928694/1‑150 (MQ=255)
         cgcgCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgtt                                                                                                                                           >  1:1645241/1‑150 (MQ=255)
         cgcgCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtt                                                                                                                                           <  2:774463/150‑1 (MQ=255)
           cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                         <  2:1952446/150‑1 (MQ=255)
             ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCgttttgt                                                                                                                                       >  1:851830/1‑150 (MQ=255)
             ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgt                                                                                                                                       <  2:698778/150‑1 (MQ=255)
              cacaTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgtt                                                                                                                                      >  2:919852/1‑150 (MQ=255)
              cacaTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgtt                                                                                                                                      <  1:4351080/150‑1 (MQ=255)
                 aTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAATTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTa                                                                                                                                   <  1:3040893/150‑1 (MQ=255)
                 aTGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTa                                                                                                                                   <  1:2855753/150‑1 (MQ=255)
                  tGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaa                                                                                                                                  <  2:2714662/150‑1 (MQ=255)
                  tGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaa                                                                                                                                  >  2:4249376/1‑150 (MQ=255)
                   ggTACCGGACGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                                 <  1:3564064/150‑1 (MQ=255)
                   ggTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                                 >  1:1756285/1‑150 (MQ=255)
                        cGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAaggagg                                                                                                                            <  2:1078637/150‑1 (MQ=255)
                          gAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGaa                                                                                                                          <  1:4428999/150‑1 (MQ=255)
                             gCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTa                                                                                                                       >  1:2554995/1‑150 (MQ=255)
                              cccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                                      <  1:2308696/150‑1 (MQ=255)
                                 aTGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCgg                                                                                                                   >  1:698098/1‑150 (MQ=255)
                                  tGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGc                                                                                                                  >  2:497419/1‑150 (MQ=255)
                                   gCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCa                                                                                                                 >  2:1537404/1‑150 (MQ=255)
                                   gCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATAAGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCa                                                                                                                 <  2:2841480/150‑1 (MQ=255)
                                    cTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAg                                                                                                                >  1:4614222/1‑150 (MQ=255)
                                    cTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAg                                                                                                                >  1:2248719/1‑150 (MQ=255)
                                     tttCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGt                                                                                                               >  1:625397/1‑150 (MQ=255)
                                     tttCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGt                                                                                                               >  1:1452952/1‑150 (MQ=255)
                                       tCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAg                                                                                                             <  2:1714086/150‑1 (MQ=255)
                                       tCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAg                                                                                                             >  2:3660076/1‑150 (MQ=255)
                                       tCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAg                                                                                                             <  2:2405117/150‑1 (MQ=255)
                                           aaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGa                                                                                                         >  1:1064050/1‑150 (MQ=255)
                                           aaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGa                                                                                                         >  1:3381424/1‑150 (MQ=255)
                                             cgacgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAAt                                                                                                       >  1:1342147/1‑150 (MQ=255)
                                              gacgCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATg                                                                                                      <  2:1795201/150‑1 (MQ=255)
                                                   cATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTAc                                                                                                 <  2:204714/150‑1 (MQ=255)
                                                   cATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTAc                                                                                                 >  2:2323505/1‑150 (MQ=255)
                                                      tttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGAt                                                                                              >  1:3929034/1‑150 (MQ=255)
                                                       ttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATg                                                                                             >  1:873770/1‑150 (MQ=255)
                                                         cGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTa                                                                                           <  2:3767425/150‑1 (MQ=255)
                                                          gATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTaa                                                                                          <  2:2979735/150‑1 (MQ=255)
                                                             ttCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcg                                                                                       <  2:4349339/150‑1 (MQ=255)
                                                             ttCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcg                                                                                       <  2:1893698/150‑1 (MQ=255)
                                                              tCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcgc                                                                                      >  1:512269/1‑150 (MQ=255)
                                                                 cAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                                   >  2:4553841/1‑150 (MQ=255)
                                                                  aaTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                                                                  <  1:1615477/150‑1 (MQ=255)
                                                                        tGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAAt                                                                            <  2:1646926/150‑1 (MQ=255)
                                                                         gCGTAACTGCCCCTCAATTCGGGTTACCAGTTCATTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATa                                                                           <  1:3153045/150‑1 (MQ=255)
                                                                          cGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAg                                                                          <  1:821488/150‑1 (MQ=255)
                                                                              aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcac                                                                      <  2:1611753/150‑1 (MQ=255)
                                                                                   cccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAc                                                                 <  2:1311635/150‑1 (MQ=255)
                                                                                   cccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAc                                                                 <  2:664985/150‑1 (MQ=255)
                                                                                    ccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAcc                                                                <  1:1631250/150‑1 (MQ=255)
                                                                                    ccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAcc                                                                >  1:2767416/1‑150 (MQ=255)
                                                                                     cTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCa                                                               >  1:3249565/1‑150 (MQ=255)
                                                                                      tCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGACGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCaa                                                              >  2:4534345/1‑150 (MQ=255)
                                                                                         aTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAg                                                           <  1:1847845/150‑1 (MQ=255)
                                                                                          tACGTGTTACCGGATCAATAGTAATTCTCTGGCTCTTGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGa                                                          <  2:1100259/150‑1 (MQ=255)
                                                                                           aCGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGaa                                                         >  2:370945/1‑150 (MQ=255)
                                                                                           aCGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGaa                                                         <  1:2916310/150‑1 (MQ=255)
                                                                                              ggTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATc                                                      >  1:1821360/1‑150 (MQ=255)
                                                                                                 tACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAAc                                                   <  2:473001/150‑1 (MQ=255)
                                                                                                   ccGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAg                                                 >  1:2182544/1‑150 (MQ=255)
                                                                                                    cGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGt                                                <  1:459639/150‑1 (MQ=255)
                                                                                                     ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                               <  2:2438286/150‑1 (MQ=255)
                                                                                                     ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTACGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                               >  2:527263/1‑150 (MQ=255)
                                                                                                        tCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttctt                                            <  1:3857869/150‑1 (MQ=255)
                                                                                                         cAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttcttc                                           <  1:3380022/150‑1 (MQ=255)
                                                                                                         cAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttcttc                                           >  2:1164142/1‑150 (MQ=255)
                                                                                                          aaTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcc                                          >  1:3706852/1‑150 (MQ=255)
                                                                                                          aaTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTACGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcc                                          >  2:1258999/1‑150 (MQ=255)
                                                                                                            tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                        >  2:3109138/1‑150 (MQ=255)
                                                                                                            tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAATGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcaca                                        >  2:3171666/1‑150 (MQ=255)
                                                                                                             aGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAg                                       <  2:3825404/150‑1 (MQ=255)
                                                                                                                aagtctcTGGCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATTGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTc                                    <  2:1930784/147‑1 (MQ=255)
                                                                                                                  ttCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCgg                                  <  1:448154/150‑1 (MQ=255)
                                                                                                                   tctcTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                                 >  2:1898602/1‑150 (MQ=255)
                                                                                                                     tctGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                               <  1:3995849/150‑1 (MQ=255)
                                                                                                                        ggCTCATGCTTTGCTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                            >  2:4587374/1‑150 (MQ=255)
                                                                                                                         gCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                           <  2:1756285/150‑1 (MQ=255)
                                                                                                                          cTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAgg                          <  1:4084288/150‑1 (MQ=255)
                                                                                                                                cTTTGCTCGCCTCATGACGATTATGATCGATTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACACTTCTTCCCAGGTCGGGAAGTAGGCGTAAc                    <  1:2470343/150‑1 (MQ=255)
                                                                                                                                 tttGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAcc                   >  1:1977824/1‑150 (MQ=255)
                                                                                                                                 tttGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAcc                   >  1:2554831/1‑150 (MQ=255)
                                                                                                                                 tttGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAcc                   <  1:4404399/150‑1 (MQ=255)
                                                                                                                                     cTCGCCTCATGACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAccgccg               <  1:1023037/150‑1 (MQ=255)
                                                                                                                                     cTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAccgccg               >  2:3484685/1‑150 (MQ=255)
                                                                                                                                       cGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGACGTAGGCGTAACCgccgcc             >  1:3672624/1‑150 (MQ=255)
                                                                                                                                        gCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGccc            <  1:2633882/150‑1 (MQ=255)
                                                                                                                                           tCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTTTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCggg         >  2:4190716/1‑150 (MQ=255)
                                                                                                                                           tCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTTTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCggg         >  2:87606/1‑150 (MQ=255)
                                                                                                                                            cATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt        >  1:4460188/1‑150 (MQ=255)
                                                                                                                                             aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt       <  1:4220984/150‑1 (MQ=255)
                                                                                                                                               gACTATTATGATCGTTTTGTTTTAACGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCACTAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGa     <  2:589750/150‑1 (MQ=255)
                                                                                                                                               gACTATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGa     >  2:680308/1‑150 (MQ=255)
                                                                                                                                               gACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGa     <  2:4443369/150‑1 (MQ=255)
                                                                                                                                                 cGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt   >  2:3552576/1‑150 (MQ=255)
                                                                                                                                                 cGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt   <  2:3337829/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATG  >  CP009273/3136222‑3136517

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: