Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,138,777 C→A 19.5% G282C (GGT→TGT)  hybO ← hydrogenase 2, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,138,7770CA19.5% 552.1 / 82.0 199G282C (GGT→TGT) hybOhydrogenase 2, small subunit
Reads supporting (aligned to +/- strand):  ref base C (93/67);  new base A (20/19);  total (113/86)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CGCCCTCTTTAGCGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTT  >  CP009273/3138633‑3138923
                                                                                                                                                |                                                                                                                                                  
cgccCTCTTTAGCGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGcc                                                                                                                                                 >  1:1073686/1‑148 (MQ=255)
  cccTCTTTAGCGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCaa                                                                                                                                               >  2:305515/1‑148 (MQ=255)
       tttAGCGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCg                                                                                                                                          <  1:2733687/148‑1 (MQ=255)
            cGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGaa                                                                                                                                     >  1:3437249/1‑148 (MQ=255)
            cGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGaa                                                                                                                                     >  1:3737072/1‑148 (MQ=255)
              ttAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAAtt                                                                                                                                   >  1:1523998/1‑148 (MQ=255)
              ttAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAAtt                                                                                                                                   >  1:1528025/1‑148 (MQ=255)
              ttAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAAtt                                                                                                                                   >  1:1746711/1‑148 (MQ=255)
               tAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTg                                                                                                                                  <  2:479662/148‑1 (MQ=255)
               tAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTg                                                                                                                                  >  2:2735929/1‑148 (MQ=255)
               tAACATCCGGTGTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTg                                                                                                                                  <  2:1991035/148‑1 (MQ=255)
                aaCATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGc                                                                                                                                 <  1:3810855/148‑1 (MQ=255)
                  cATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAg                                                                                                                               >  2:2316657/1‑148 (MQ=255)
                    tCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAATGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCg                                                                                                                             <  1:2983253/148‑1 (MQ=255)
                    tCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCg                                                                                                                             <  2:776756/148‑1 (MQ=255)
                     ccGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGt                                                                                                                            >  1:1051858/1‑148 (MQ=255)
                       ggTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTg                                                                                                                          >  1:3482489/1‑148 (MQ=255)
                       ggTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTg                                                                                                                          >  2:2980366/1‑148 (MQ=255)
                          ttCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGc                                                                                                                       <  1:772616/148‑1 (MQ=255)
                          ttCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGc                                                                                                                       >  1:455413/1‑148 (MQ=255)
                           tCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCa                                                                                                                      <  2:526073/148‑1 (MQ=255)
                             tgtgAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGt                                                                                                                    <  2:1301534/148‑1 (MQ=255)
                              gtgAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGtt                                                                                                                   >  2:2437587/1‑148 (MQ=255)
                              gtgAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGtt                                                                                                                   >  2:2437545/1‑148 (MQ=255)
                               tgAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTg                                                                                                                  >  2:1317310/1‑148 (MQ=255)
                               cgacctcgGAGTTTGAGGCTCCACCCTGGCTAGATGCTCGATGCCTTACTGGCAGGCGATACCCTCCGCGTTTCCGCCATATCTTTCTTGAGCTCACGCCACAGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTg                                                                                                                  <  2:1068687/142‑1 (MQ=255)
                                      ggAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAg                                                                                                           <  1:2108783/148‑1 (MQ=255)
                                       gAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGt                                                                                                          <  1:290694/148‑1 (MQ=255)
                                       gAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTGACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGt                                                                                                          <  2:1444232/148‑1 (MQ=255)
                                         gTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGttt                                                                                                        <  2:659447/148‑1 (MQ=255)
                                         gTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGttt                                                                                                        <  1:167483/148‑1 (MQ=255)
                                           ttGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGACTTGCAGCGTTGAGCAGTTGCCGTAAGTTTCt                                                                                                      >  2:3280495/1‑148 (MQ=255)
                                             gATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTgg                                                                                                    >  2:1229905/1‑148 (MQ=255)
                                              aTTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGc                                                                                                   >  2:2250193/1‑148 (MQ=255)
                                               ttttCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGcc                                                                                                  >  2:2345634/1‑148 (MQ=255)
                                                tttCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGccc                                                                                                 <  2:2860199/148‑1 (MQ=255)
                                                tttCGACGGTGGCAAGCTGATGGATGCCTTTAGGGAAGCCGAGACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGccc                                                                                                 <  2:2114662/148‑1 (MQ=255)
                                                 ttCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCt                                                                                                >  1:2032343/1‑148 (MQ=255)
                                                  tCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCtt                                                                                               <  2:2273692/148‑1 (MQ=255)
                                                     aCGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTAc                                                                                            >  1:991762/1‑148 (MQ=255)
                                                     aCGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTAc                                                                                            <  2:2673648/148‑1 (MQ=255)
                                                      cGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACa                                                                                           <  2:3317583/148‑1 (MQ=255)
                                                                  gatggatgCCTTTATGGAAGCCGTTACCTTCTGCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACCTCGCAGAATTGCAGCGTGGAGCAGTTGCCGTAAGTGTCTGGCCCTGTACAGCCGAGGTGGTa                                                                               >  2:1278857/1‑148 (MQ=255)
                                                                  gatggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTa                                                                               >  2:3567316/1‑148 (MQ=255)
                                                                   atggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCACTCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTAc                                                                              <  2:1340777/148‑1 (MQ=255)
                                                                   atggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTAc                                                                              >  2:3132219/1‑148 (MQ=255)
                                                                    tggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACa                                                                             >  2:1890263/1‑148 (MQ=255)
                                                                    tggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACa                                                                             >  2:355152/1‑148 (MQ=255)
                                                                    tggatgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACa                                                                             >  1:416709/1‑148 (MQ=255)
                                                                        tgCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCa                                                                         >  1:464653/1‑148 (MQ=255)
                                                                         gCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCAc                                                                        <  2:1073686/148‑1 (MQ=255)
                                                                         gCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCAc                                                                        <  1:735808/148‑1 (MQ=255)
                                                                            tttATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAg                                                                     >  1:562988/1‑148 (MQ=255)
                                                                             ttATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGc                                                                    >  1:1118953/1‑148 (MQ=255)
                                                                                    aGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTcgcg                                                             >  1:1883994/1‑148 (MQ=255)
                                                                                      ccGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGt                                                           >  2:2372233/1‑148 (MQ=255)
                                                                                       cGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTg                                                          <  1:2250193/148‑1 (MQ=255)
                                                                                          tACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGcc                                                       <  1:1229905/148‑1 (MQ=255)
                                                                                          tACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGcc                                                       <  2:289262/148‑1 (MQ=255)
                                                                                           aCCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCt                                                      >  2:1980418/1‑148 (MQ=255)
                                                                                           aCCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCt                                                      >  1:142763/1‑148 (MQ=255)
                                                                                           aCCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCt                                                      <  1:642307/148‑1 (MQ=255)
                                                                                              ttcttcGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTtca                                                   >  1:2418303/1‑148 (MQ=255)
                                                                                               tcttcGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTtcat                                                  <  1:3067200/148‑1 (MQ=255)
                                                                                                     ttACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGa                                            <  2:679086/148‑1 (MQ=255)
                                                                                                       aCAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAAc                                          >  2:387228/1‑148 (MQ=255)
                                                                                                        cAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGAAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAAct                                         <  1:102050/148‑1 (MQ=255)
                                                                                                         aGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAActc                                        <  2:1612260/148‑1 (MQ=255)
                                                                                                          gCCATAGCAAGGGTGACCACTCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAActct                                       <  1:2001882/148‑1 (MQ=255)
                                                                                                          gCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAActct                                       >  2:1803462/1‑148 (MQ=255)
                                                                                                          gCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAActct                                       >  2:3627875/1‑148 (MQ=255)
                                                                                                           ccATAGCAAGGGTGACCAATCTCCACCGGCCACACACAGCGATCATCGCAGAATTGAAGCGTTGAGCAGTTGCCGTAATTTTCTTGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCtt                                      <  2:1786606/148‑1 (MQ=255)
                                                                                                            cATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCttt                                     <  1:2361358/148‑1 (MQ=255)
                                                                                                            cATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCttt                                     <  2:2614319/148‑1 (MQ=255)
                                                                                                             aTAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTg                                    >  2:252331/1‑148 (MQ=255)
                                                                                                                 cAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCaa                                >  1:431026/1‑148 (MQ=255)
                                                                                                                  aaGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCaaa                               <  1:2681592/148‑1 (MQ=255)
                                                                                                                  aaGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCaaa                               <  2:2158003/148‑1 (MQ=255)
                                                                                                                  aaGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCaaa                               <  1:2394733/148‑1 (MQ=255)
                                                                                                                  aaGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCaaa                               <  1:2680831/148‑1 (MQ=255)
                                                                                                                    gggTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAAc                             >  2:1030418/1‑148 (MQ=255)
                                                                                                                    gggTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCAGAGGTTGTACAGGCACCAGCCTTCGAGGTGGCCATCATCACCGAACTCTTTGGCAAAAc                             >  2:3385331/1‑148 (MQ=255)
                                                                                                                     ggTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACg                            >  1:1369820/1‑148 (MQ=255)
                                                                                                                     ggTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACg                            >  1:2949736/1‑148 (MQ=255)
                                                                                                                         aCCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCa                        >  2:903487/1‑148 (MQ=255)
                                                                                                                           cAATCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGc                      <  2:2379840/148‑1 (MQ=255)
                                                                                                                            aaTCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCa                     >  1:887328/1‑148 (MQ=255)
                                                                                                                            aaTCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCa                     >  1:423645/1‑148 (MQ=255)
                                                                                                                            aaTCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCa                     >  1:3145273/1‑148 (MQ=255)
                                                                                                                            aaTCGCCACCGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCa                     >  2:2412720/1‑148 (MQ=255)
                                                                                                                                gCCACCGGCCACACACCGCCAACATCGCAGATTTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGTGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGa                 <  1:3295078/148‑1 (MQ=255)
                                                                                                                                gCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGa                 <  2:505254/148‑1 (MQ=255)
                                                                                                                                  cACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAg               <  2:3468314/148‑1 (MQ=255)
                                                                                                                                   aCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGt              <  1:1468028/148‑1 (MQ=255)
                                                                                                                                    ccGGCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTg             <  2:1173021/148‑1 (MQ=255)
                                                                                                                                      ggCCACACACAGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCg           <  1:3164917/148‑1 (MQ=255)
                                                                                                                                       gCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCgg          <  2:3714732/148‑1 (MQ=255)
                                                                                                                                          acacacAGCCAACATCGCCTATTGGCATCGTTGAGCAGTTGCGGTAAGTTTCGGGCCGTTTACAGCCGTGGTGGTAAAGGCACCAGCCTTCGCGGTGGACTTCATAACCGAACACTTTGGCAAAACGACCAGCATCGAAGTGCGGGCg       <  1:1599396/148‑1 (MQ=255)
                                                                                                                                           cacacCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGa      >  2:909237/1‑148 (MQ=255)
                                                                                                                                             cacCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACg    <  1:3627875/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGCCCTCTTTAGCGTTAACATCCGGTTTCTGTGAACGCGGAGTTTGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACAGGCACCAGCCTTCGCGGTGGCCTTCATCACCGAACTCTTTGGCAAAACGACCAGCATCGAAGTGCGGGCGACGTT  >  CP009273/3138633‑3138923

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: