Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 61,967 A→C 8.6% L433R (CTG→CGG)  rapA ← RNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927361,9670AC 8.6% 432.8 / 33.2 163L433R (CTG→CGG) rapARNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base A (73/76);  new base C (3/11);  total (76/87)
Fisher's exact test for biased strand distribution p-value = 5.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 164 total aligned reads displayed.

TATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACAGCG  >  CP009273/61824‑62114
                                                                                                                                               |                                                                                                                                                   
 aTCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTcg                                                                                                                                                <  2:2784976/148‑1 (MQ=255)
 aTCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTcg                                                                                                                                                >  1:1458872/1‑148 (MQ=255)
   cACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTcgcg                                                                                                                                              >  1:3216612/1‑148 (MQ=255)
    aCCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGt                                                                                                                                             >  1:473927/1‑148 (MQ=255)
     ccTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGtt                                                                                                                                            >  1:1003404/1‑148 (MQ=255)
     ccTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGtt                                                                                                                                            >  2:442447/1‑148 (MQ=255)
     ccTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGtt                                                                                                                                            >  2:442374/1‑148 (MQ=255)
        tCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCg                                                                                                                                         >  2:2099520/1‑148 (MQ=255)
        tCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCg                                                                                                                                         >  2:3191485/1‑148 (MQ=255)
         cAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCgg                                                                                                                                        <  1:1432264/148‑1 (MQ=255)
           aaTTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGa                                                                                                                                      <  1:1108863/148‑1 (MQ=255)
           aaTTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGa                                                                                                                                      >  1:2361092/1‑148 (MQ=255)
            cttCCTGATACATCCGCTCCGGGGAGAGCATCTCGCGAGCACGATCTTCCGCACGTTGACGTGCGCCCATAATGCCGGAGACTTGAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGaa                                                                                                                                     <  2:1901845/147‑1 (MQ=255)
             ttCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAAt                                                                                                                                    <  1:2898209/148‑1 (MQ=255)
             ttCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAAt                                                                                                                                    >  1:1352363/1‑148 (MQ=255)
              tCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATc                                                                                                                                   <  2:3465920/148‑1 (MQ=255)
                  gATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCttt                                                                                                                               <  1:1297185/148‑1 (MQ=255)
                  gATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCttt                                                                                                                               <  1:3455918/148‑1 (MQ=255)
                  gATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCttt                                                                                                                               <  1:3455932/148‑1 (MQ=255)
                    tAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTca                                                                                                                             <  2:1842504/148‑1 (MQ=255)
                        tACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACAcc                                                                                                                         >  2:14257/1‑148 (MQ=255)
                                gggTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTAcgcg                                                                                                                 <  2:2353521/148‑1 (MQ=255)
                                  gTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCATTTCACACCGTTACGCgtg                                                                                                               <  1:956346/148‑1 (MQ=255)
                                    agagCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGtt                                                                                                             <  2:622906/148‑1 (MQ=255)
                                     gagCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTa                                                                                                            <  1:2546249/148‑1 (MQ=255)
                                      agCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTAc                                                                                                           >  2:2103114/1‑148 (MQ=255)
                                         atatCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGa                                                                                                        <  2:846116/148‑1 (MQ=255)
                                          tatCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGaa                                                                                                       >  2:772283/1‑148 (MQ=255)
                                             cgcgAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAg                                                                                                    >  2:1138050/1‑148 (MQ=255)
                                             cgcgAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAg                                                                                                    <  2:78442/148‑1 (MQ=255)
                                               cgAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCa                                                                                                  <  1:2545878/148‑1 (MQ=255)
                                               cgAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGNAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCa                                                                                                  <  1:1801739/148‑1 (MQ=255)
                                                 aGCACGATCTTCCGCACTTTTACGTGCGCCCAGAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCAcg                                                                                                <  2:2323356/148‑1 (MQ=255)
                                                   cACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCAcgcg                                                                                              >  2:480114/1‑148 (MQ=255)
                                                    aCGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCgg                                                                                             <  2:494545/148‑1 (MQ=255)
                                                    aCGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCgg                                                                                             <  1:2830167/148‑1 (MQ=255)
                                                    aCGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCgg                                                                                             <  1:3115149/148‑1 (MQ=255)
                                                         cTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGt                                                                                        >  2:724866/1‑148 (MQ=255)
                                                         cTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGt                                                                                        >  2:2506926/1‑148 (MQ=255)
                                                          ttCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTg                                                                                       <  1:2582106/148‑1 (MQ=255)
                                                          ttCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTg                                                                                       <  1:2768441/148‑1 (MQ=255)
                                                           tCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGc                                                                                      >  2:3225451/1‑148 (MQ=255)
                                                             cGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCg                                                                                    <  2:2192671/148‑1 (MQ=255)
                                                              gCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGt                                                                                   >  1:2031369/1‑148 (MQ=255)
                                                               cACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGACCAGCACGCGGCTGGTGCCGTg                                                                                  >  2:2613287/1‑148 (MQ=255)
                                                               cACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTg                                                                                  <  2:3494163/148‑1 (MQ=255)
                                                               cACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTg                                                                                  >  1:1939168/1‑148 (MQ=255)
                                                               cACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTg                                                                                  <  1:9134/148‑1 (MQ=255)
                                                                  ttttACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCg                                                                               >  2:1324338/1‑148 (MQ=255)
                                                                   tttACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGa                                                                              <  2:2255191/148‑1 (MQ=255)
                                                                     tACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATc                                                                            <  2:2055873/148‑1 (MQ=255)
                                                                     tACGGGCGCCCCTAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATc                                                                            <  1:1984673/148‑1 (MQ=255)
                                                                             cccATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCa                                                                    >  1:1028970/1‑148 (MQ=255)
                                                                              ccATAATTCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTCGCGGCATCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCAt                                                                   <  2:1200682/148‑1 (MQ=255)
                                                                                 tAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGa                                                                <  2:1245679/148‑1 (MQ=255)
                                                                                 tAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGa                                                                <  2:2079712/148‑1 (MQ=255)
                                                                                   aTGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGaaa                                                              >  2:3027448/1‑148 (MQ=255)
                                                                                    tGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAAc                                                             >  1:2983592/1‑148 (MQ=255)
                                                                                    tGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAAc                                                             <  1:1273570/148‑1 (MQ=255)
                                                                                    tGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAAc                                                             <  1:1274156/148‑1 (MQ=255)
                                                                                    tGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAAc                                                             >  1:2145273/1‑148 (MQ=255)
                                                                                         gagaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCt                                                        >  2:659879/1‑148 (MQ=255)
                                                                                         gagaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCt                                                        >  2:659812/1‑148 (MQ=255)
                                                                                          agaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTc                                                       <  2:270427/148‑1 (MQ=255)
                                                                                           gaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTcc                                                      <  2:1003404/148‑1 (MQ=255)
                                                                                             cTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTg                                                    <  1:2593422/148‑1 (MQ=255)
                                                                                                    aGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGgcg                                             <  2:1777243/148‑1 (MQ=255)
                                                                                                      ccGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCTCACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGct                                           >  1:1175528/1‑148 (MQ=255)
                                                                                                      ccGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGct                                           <  2:1028970/148‑1 (MQ=255)
                                                                                                      ccGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGct                                           >  2:3495544/1‑148 (MQ=255)
                                                                                                      ccGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGAGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACTCGTCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGct                                           <  2:3124305/148‑1 (MQ=255)
                                                                                                       cGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGctc                                          >  2:3532066/1‑148 (MQ=255)
                                                                                                         tCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTg                                        >  1:1429697/1‑148 (MQ=255)
                                                                                                         tCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTg                                        >  1:1428318/1‑148 (MQ=255)
                                                                                                         tCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTg                                        >  1:1429042/1‑148 (MQ=255)
                                                                                                          cTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTgg                                       <  1:3100221/148‑1 (MQ=255)
                                                                                                             aTACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCa                                    >  2:486591/1‑148 (MQ=255)
                                                                                                               aCTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATc                                  >  2:1274734/1‑148 (MQ=255)
                                                                                                               aCTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATc                                  >  2:3046203/1‑148 (MQ=255)
                                                                                                                cTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCt                                 <  1:2887127/148‑1 (MQ=255)
                                                                                                                  gCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtc                               >  1:878521/1‑148 (MQ=255)
                                                                                                                  gCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtc                               >  1:878614/1‑148 (MQ=255)
                                                                                                                         tAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctgtcg                        <  1:1315428/148‑1 (MQ=255)
                                                                                                                          aGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctgtcgc                       >  1:1815599/1‑148 (MQ=255)
                                                                                                                          aGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctgtcgc                       >  1:976631/1‑148 (MQ=255)
                                                                                                                            cGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctgtcgctg                     >  2:1391947/1‑148 (MQ=255)
                                                                                                                             ggCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTcgctgtcgctgt                    <  1:659812/148‑1 (MQ=255)
                                                                                                                             ggCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTcgctgtcgctgt                    <  1:659869/148‑1 (MQ=255)
                                                                                                                              gCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGtt                   <  2:66329/148‑1 (MQ=255)
                                                                                                                                aGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtg                 <  1:2726474/148‑1 (MQ=255)
                                                                                                                                 gCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAATCCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgc                >  1:239501/1‑148 (MQ=255)
                                                                                                                                 gCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgc                >  1:3029385/1‑148 (MQ=255)
                                                                                                                                    tAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgctgc             <  1:2622838/148‑1 (MQ=255)
                                                                                                                                      aTGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcct           >  1:1553954/1‑148 (MQ=255)
                                                                                                                                       tGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctg          >  1:2054248/1‑148 (MQ=255)
                                                                                                                                       tGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctg          >  1:3052337/1‑148 (MQ=255)
                                                                                                                                        ggTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctgc         <  2:2983592/148‑1 (MQ=255)
                                                                                                                                           gtgCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAAc      <  2:650586/148‑1 (MQ=255)
                                                                                                                                           gtgCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAAc      >  1:660144/1‑148 (MQ=255)
                                                                                                                                               aGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACAGCg  <  1:1138050/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACAGCG  >  CP009273/61824‑62114

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: