Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,136,365 G→T 5.6% H387N (CAT→AAT)  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,3650GT 5.6% 871.1 / 17.3 231H387N (CAT→AAT) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base G (106/112);  new base T (8/5);  total (114/117)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 232 total aligned reads displayed.

ATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGG  >  CP009273/3136218‑3136510
                                                                                                                                                   |                                                                                                                                                 
aTCTCTTCCATGCTGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATg                                                                                                                                                   <  2:2001354/148‑1 (MQ=255)
  ctctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGAc                                                                                                                                                 >  1:1457105/1‑148 (MQ=255)
   tctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACg                                                                                                                                                <  2:2634661/148‑1 (MQ=255)
    ctTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGa                                                                                                                                               <  1:772586/148‑1 (MQ=255)
      tCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAtt                                                                                                                                             <  1:2456743/148‑1 (MQ=255)
      tCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGAtt                                                                                                                                             <  1:982978/148‑1 (MQ=255)
       ccATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTa                                                                                                                                            >  2:1453600/1‑148 (MQ=255)
             cgcgCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCg                                                                                                                                      <  1:476631/148‑1 (MQ=255)
               cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtt                                                                                                                                    <  1:1187616/148‑1 (MQ=255)
               cgcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtt                                                                                                                                    <  1:3341953/148‑1 (MQ=255)
                gcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttt                                                                                                                                   <  2:1490755/148‑1 (MQ=255)
                 ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                  >  1:125278/1‑148 (MQ=255)
                 ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgtttt                                                                                                                                  >  2:2141018/1‑148 (MQ=255)
                         tACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                          <  1:2798173/148‑1 (MQ=255)
                         tACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTaaa                                                                                                                          <  2:109499/148‑1 (MQ=255)
                          aCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAag                                                                                                                         <  1:109023/148‑1 (MQ=255)
                          aCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAag                                                                                                                         <  2:616247/148‑1 (MQ=255)
                           ccGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAagg                                                                                                                        <  1:2736638/148‑1 (MQ=255)
                           ccGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAagg                                                                                                                        <  1:2707896/148‑1 (MQ=255)
                               aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAggagga                                                                                                                    <  2:2176411/148‑1 (MQ=255)
                               aaGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAggagga                                                                                                                    <  2:2176444/148‑1 (MQ=255)
                                 gCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAg                                                                                                                  <  1:3200944/148‑1 (MQ=255)
                                 gCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTAAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAg                                                                                                                  >  1:965350/1‑148 (MQ=255)
                                  cccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGt                                                                                                                 <  1:1346219/148‑1 (MQ=255)
                                   ccATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATGGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTa                                                                                                                <  1:2596007/148‑1 (MQ=255)
                                    cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                               <  2:1003936/148‑1 (MQ=255)
                                    cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                               >  2:2299460/1‑148 (MQ=255)
                                    cATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTAt                                                                                                               >  2:2663988/1‑148 (MQ=255)
                                     aTGCTTTCGACACGACACCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATc                                                                                                              <  1:629488/148‑1 (MQ=255)
                                     aTGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATc                                                                                                              >  2:2320523/1‑148 (MQ=255)
                                      tGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCg                                                                                                             >  2:2393610/1‑148 (MQ=255)
                                      tGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCg                                                                                                             >  2:1044542/1‑148 (MQ=255)
                                       gCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCgg                                                                                                            >  1:3109435/1‑148 (MQ=255)
                                        cTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGc                                                                                                           >  2:2005405/1‑148 (MQ=255)
                                        cTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGc                                                                                                           >  2:948667/1‑148 (MQ=255)
                                              aaaCGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGa                                                                                                     >  1:1689714/1‑148 (MQ=255)
                                                 ggacgCGATTTTCTTTTTCTCAATCGATGCTTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTTTGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGa                                                                                                  <  1:2752439/147‑1 (MQ=255)
                                                     gCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGa                                                                                              <  1:1919676/148‑1 (MQ=255)
                                                      ccATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAg                                                                                             >  2:2463635/1‑148 (MQ=255)
                                                      ccATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAg                                                                                             <  2:577244/148‑1 (MQ=255)
                                                      ccATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAg                                                                                             >  1:1932958/1‑148 (MQ=255)
                                                       cATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGt                                                                                            >  2:2106703/1‑148 (MQ=255)
                                                        aTTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTa                                                                                           >  1:1519030/1‑148 (MQ=255)
                                                           ttCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGa                                                                                        >  2:669095/1‑148 (MQ=255)
                                                            tCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTAt                                                                                       >  2:1588756/1‑148 (MQ=255)
                                                              gATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGt                                                                                     >  1:2647115/1‑148 (MQ=255)
                                                                   cGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcg                                                                                <  1:3371500/148‑1 (MQ=255)
                                                                   cGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAgcg                                                                                <  1:2921312/148‑1 (MQ=255)
                                                                       aTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                            <  2:2556582/148‑1 (MQ=255)
                                                                       aTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                            <  1:3463088/148‑1 (MQ=255)
                                                                       aTCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCaaa                                                                            >  1:2317958/1‑148 (MQ=255)
                                                                        tCGATGCGTAACTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                                                           <  2:1077759/148‑1 (MQ=255)
                                                                        tCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAAt                                                                           >  1:2225348/1‑148 (MQ=255)
                                                                        tCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAAt                                                                           <  2:3203944/148‑1 (MQ=255)
                                                                             gCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                      <  2:1943579/148‑1 (MQ=255)
                                                                             gCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCaa                                                                      >  2:3174746/1‑148 (MQ=255)
                                                                              cGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAAt                                                                     <  1:374859/148‑1 (MQ=255)
                                                                                 aaaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGc                                                                  >  1:2550670/1‑148 (MQ=255)
                                                                                  aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                                                                 <  2:1113461/148‑1 (MQ=255)
                                                                                  aaTGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGcc                                                                 <  2:2891593/148‑1 (MQ=255)
                                                                                    tGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCcac                                                               >  2:1228590/1‑148 (MQ=255)
                                                                                     gCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCcaca                                                              >  1:1620128/1‑148 (MQ=255)
                                                                                        ccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACaaaa                                                           >  2:378927/1‑148 (MQ=255)
                                                                                        ccTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACaaaa                                                           >  1:3015752/1‑148 (MQ=255)
                                                                                          tCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAcc                                                         >  1:2012853/1‑148 (MQ=255)
                                                                                          tCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAAcc                                                         <  2:2522975/148‑1 (MQ=255)
                                                                                           cAATACGGTTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCCACAAAACCa                                                        <  2:1911778/148‑1 (MQ=255)
                                                                                           cAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCa                                                        >  2:227549/1‑148 (MQ=255)
                                                                                            aaTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCCACAAAACCaa                                                       <  1:2833626/148‑1 (MQ=255)
                                                                                                  ggTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGaaa                                                 >  1:51706/1‑148 (MQ=255)
                                                                                                  ggTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGaaa                                                 >  2:1853344/1‑148 (MQ=255)
                                                                                                       ccGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAAc                                            <  2:1769939/148‑1 (MQ=255)
                                                                                                         ggATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAg                                          >  1:441569/1‑148 (MQ=255)
                                                                                                          gATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGt                                         >  1:312882/1‑148 (MQ=255)
                                                                                                          gATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGt                                         <  2:3557490/148‑1 (MQ=255)
                                                                                                           aTCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGtt                                        <  2:844246/148‑1 (MQ=255)
                                                                                                              aaTAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGttctt                                     <  2:2070035/148‑1 (MQ=255)
                                                                                                                tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcc                                   >  1:3527127/1‑148 (MQ=255)
                                                                                                                tAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTcc                                   <  1:2141018/148‑1 (MQ=255)
                                                                                                                 aGTAATTCTCTGGCTCATGCTGTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGTAGTATCGGCTGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTccc                                  <  2:3391665/148‑1 (MQ=255)
                                                                                                                  gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                 <  1:235927/148‑1 (MQ=255)
                                                                                                                  gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                 <  1:2634040/148‑1 (MQ=255)
                                                                                                                  gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCa                                 <  1:1087409/148‑1 (MQ=255)
                                                                                                                  gTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAAGCAACAGTTCTTCCCa                                 >  1:3319530/1‑148 (MQ=255)
                                                                                                                   tAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAg                                >  2:2482046/1‑148 (MQ=255)
                                                                                                                   tAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAg                                >  1:3414444/1‑148 (MQ=255)
                                                                                                                       tctcTGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCg                            >  1:1631261/1‑148 (MQ=255)
                                                                                                                       tctcTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAAt                                                      <  1:2481289/122‑1 (MQ=255)
                                                                                                                         tctGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCggg                          >  2:3092044/1‑148 (MQ=255)
                                                                                                                          ctGGCTCATGCTTTGCTCGCCTCATTACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACTATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGa                         <  1:1975567/148‑1 (MQ=255)
                                                                                                                            ggCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAg                       >  1:75052/1‑148 (MQ=255)
                                                                                                                             gCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGt                      <  2:1652219/148‑1 (MQ=255)
                                                                                                                              cTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTa                     <  1:1853344/148‑1 (MQ=255)
                                                                                                                               tCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                    <  1:448421/148‑1 (MQ=255)
                                                                                                                                  tGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                 <  2:3145415/148‑1 (MQ=255)
                                                                                                                                  tGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                 <  2:1370753/148‑1 (MQ=255)
                                                                                                                                  tGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                 >  2:708081/1‑148 (MQ=255)
                                                                                                                                           cGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTACCAGGTCGGGAAGTAGGCGTAAccgccg        >  2:3229889/1‑148 (MQ=255)
                                                                                                                                             ccTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGGTCTTCCCAGGTCGGGAAGTAGGCGTAACCgccgcc      <  2:2804181/148‑1 (MQ=255)
                                                                                                                                                 aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCggg  <  1:2044472/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGG  >  CP009273/3136218‑3136510

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: