Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 868,897 G→T 13.4% D155Y (GAC→TAC)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,8970GT13.4% 1011.4 / 87.5 307D155Y (GAC→TAC) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (133/132);  new base T (13/29);  total (146/161)
Fisher's exact test for biased strand distribution p-value = 2.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

ATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGAAC  >  CP009273/868748‑869046
                                                                                                                                                     |                                                                                                                                                     
 taaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGa                                                                                                                                                      <  1:3396640/150‑1 (MQ=255)
  aaataTTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTAc                                                                                                                                                     >  1:4325829/1‑150 (MQ=255)
  aaataTTACCTTTTTTCGCGGAGGGCGAATCATAAAATCGTTGTTTAACATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGAc                                                                                                                                                     <  1:169829/150‑1 (MQ=255)
  aaataTAATCTTTTTTCGCGGCGGGCTAATGATAAAATCGTTGTTTATCATCAGTAACAAGTAAATGGACCGCTTGTCGACGAGCCAGGGCGGGTTATTTTGGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTAc                                                                                                                                                     <  1:186370/150‑1 (MQ=255)
    atatTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGActc                                                                                                                                                   >  2:2989928/1‑150 (MQ=255)
     tattatCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGActct                                                                                                                                                  >  2:936305/1‑150 (MQ=255)
        tatCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATGGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaac                                                                                                                                               <  1:3688989/150‑1 (MQ=255)
         atCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaaca                                                                                                                                              >  1:2195938/1‑150 (MQ=255)
          tCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGTAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaacaa                                                                                                                                             <  1:2953130/150‑1 (MQ=255)
           cTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTaacaac                                                                                                                                            >  1:1046821/1‑150 (MQ=255)
             tttttCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAg                                                                                                                                          >  2:60212/1‑150 (MQ=255)
              ttttCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGc                                                                                                                                         <  2:3394866/150‑1 (MQ=255)
              ttttCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATAGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGc                                                                                                                                         >  2:3994286/1‑150 (MQ=255)
                ttCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTg                                                                                                                                       >  1:3041085/1‑150 (MQ=255)
                  cgcgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGa                                                                                                                                     >  1:696354/1‑150 (MQ=255)
                  cgcgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGa                                                                                                                                     >  1:4448246/1‑150 (MQ=255)
                   gcgGAGGTCTAATGATAAAATCGTTGTTTATCATAAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGaa                                                                                                                                    <  2:3220129/150‑1 (MQ=255)
                    cgGAGGGCTAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAAt                                                                                                                                   <  2:141176/150‑1 (MQ=255)
                    cgGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAAt                                                                                                                                   <  2:141157/150‑1 (MQ=255)
                     ggAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATc                                                                                                                                  >  1:903039/1‑150 (MQ=255)
                     ggAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATc                                                                                                                                  >  1:2979654/1‑150 (MQ=255)
                      gAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCTGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATcc                                                                                                                                 <  2:2135881/150‑1 (MQ=255)
                      gAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACACCAGCTGGAATcc                                                                                                                                 >  2:3550509/1‑150 (MQ=255)
                       aGGGCGAATGATAAAATGGTTGTTTTTAATCAGGAACAATTAAATGGACCGCTTGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCg                                                                                                                                <  1:644431/150‑1 (MQ=255)
                       aGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCg                                                                                                                                <  1:713953/150‑1 (MQ=255)
                         ggCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGc                                                                                                                              <  1:2436192/150‑1 (MQ=255)
                         ggCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGc                                                                                                                              >  2:3487925/1‑150 (MQ=255)
                           cGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCAt                                                                                                                            <  2:23579/150‑1 (MQ=255)
                           cGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCAt                                                                                                                            >  1:1782277/1‑150 (MQ=255)
                            gAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCaga                                                                                                                           >  1:4194992/1‑148 (MQ=255)
                             aaTGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATAc                                                                                                                          <  2:376700/150‑1 (MQ=255)
                             aaTGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATAc                                                                                                                          >  1:3561703/1‑150 (MQ=255)
                                 aTAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAAtt                                                                                                                      <  2:4486865/150‑1 (MQ=255)
                                  tAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTg                                                                                                                     <  2:1272319/150‑1 (MQ=255)
                                   aaaaTCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGc                                                                                                                    >  1:2558592/1‑150 (MQ=255)
                                       tCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTg                                                                                                                <  1:274302/150‑1 (MQ=255)
                                        cGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGa                                                                                                               >  2:1229051/1‑150 (MQ=255)
                                             tttATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTc                                                                                                          <  1:1865898/150‑1 (MQ=255)
                                              ttATCATCAGGAACAATTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTcc                                                                                                         <  1:563613/150‑1 (MQ=255)
                                               tATCATCAGGAACAAGTAAACGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCg                                                                                                        <  2:4880/150‑1 (MQ=255)
                                                  catcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCCTACAATTGCTTTGATCGTCCGGGt                                                                                                     >  1:3116688/1‑150 (MQ=255)
                                                  catcaGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGt                                                                                                     >  1:1206649/1‑150 (MQ=255)
                                                       ggAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTATGCCGATTTTAACCCATCGAAACGACCGTGGTCTGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTaaaaa                                                                                                >  1:2928460/1‑150 (MQ=255)
                                                           cAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGt                                                                                            >  2:2369430/1‑150 (MQ=255)
                                                           cAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGt                                                                                            >  2:66566/1‑150 (MQ=255)
                                                            aaGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTa                                                                                           >  1:3602650/1‑150 (MQ=255)
                                                            aaGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTa                                                                                           >  1:2812185/1‑150 (MQ=255)
                                                            aaGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTa                                                                                           >  1:2915868/1‑150 (MQ=255)
                                                               tAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCTATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTAtct                                                                                        <  1:2344064/150‑1 (MQ=255)
                                                                    ggACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtttt                                                                                   >  2:747597/1‑150 (MQ=255)
                                                                    ggACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCtttt                                                                                   >  2:1059757/1‑150 (MQ=255)
                                                                     gACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTa                                                                                  <  1:3904590/150‑1 (MQ=255)
                                                                       ccGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACg                                                                                <  1:2856027/150‑1 (MQ=255)
                                                                         gCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCt                                                                              >  1:1318134/1‑150 (MQ=255)
                                                                          cTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTa                                                                             <  2:829678/150‑1 (MQ=255)
                                                                          cTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTa                                                                             <  2:315983/150‑1 (MQ=255)
                                                                           ttGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTac                                                                            <  1:3487925/150‑1 (MQ=255)
                                                                            tGTCGACGAGTCAGGGCGGTTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTaca                                                                           <  1:2340101/150‑1 (MQ=255)
                                                                             gTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTacac                                                                          <  2:2668891/150‑1 (MQ=255)
                                                                             gTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTacac                                                                          <  2:372701/150‑1 (MQ=255)
                                                                               cgacgaGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACAt                                                                        >  1:2240555/1‑150 (MQ=255)
                                                                               cgacgaGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACAt                                                                        >  1:2374892/1‑150 (MQ=255)
                                                                               cgacgaGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACAt                                                                        >  1:4224979/1‑150 (MQ=255)
                                                                                    aGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACg                                                                   >  2:203698/1‑150 (MQ=255)
                                                                                     gTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAGCAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACgg                                                                  >  1:3826700/1‑150 (MQ=255)
                                                                                     gTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACgg                                                                  >  2:4324441/1‑150 (MQ=255)
                                                                                         gggCGGTTTATTTTTTCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCaaa                                                              <  1:1213996/150‑1 (MQ=255)
                                                                                           gCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGa                                                            >  2:2362038/1‑150 (MQ=255)
                                                                                            cGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGAt                                                           <  1:3333693/150‑1 (MQ=255)
                                                                                              ggTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCa                                                         >  1:2831919/1‑150 (MQ=255)
                                                                                               gTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCAc                                                        <  2:1105176/150‑1 (MQ=255)
                                                                                                  aTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGAt                                                     <  1:2031701/150‑1 (MQ=255)
                                                                                                  aTTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGAt                                                     <  2:4064919/150‑1 (MQ=255)
                                                                                                  aTTTTTGCCGATTTTAACCAATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGAt                                                     <  1:1590775/150‑1 (MQ=255)
                                                                                                     tttGCCGATTTTAACGCATCGCAACGACCGTTGTCGGTGGCTTCAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTg                                                  >  1:3686615/1‑150 (MQ=255)
                                                                                                     tttGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTg                                                  >  2:1276032/1‑150 (MQ=255)
                                                                                                       tGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGtt                                                >  2:1479791/1‑150 (MQ=255)
                                                                                                         ccGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAg                                              >  1:3274337/1‑150 (MQ=255)
                                                                                                                 aaCCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGAt                                      >  2:4357881/1‑150 (MQ=255)
                                                                                                                     cATCTAAACGACCGTGGTCGGTTGCTTCAGATGACTCTAACAACAGATGGAATCCGGCATACAATTGATTTGATCGTCCGGGTAAAAAATGTATCACTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCCGGGGATaaaa                                  >  2:281278/1‑150 (MQ=255)
                                                                                                                     cATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATaaaa                                  >  2:1598949/1‑150 (MQ=255)
                                                                                                                      aTCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAAt                                 >  2:2580725/1‑150 (MQ=255)
                                                                                                                             cGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTc                          <  2:1750505/150‑1 (MQ=255)
                                                                                                                              gACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCg                         >  2:851296/1‑150 (MQ=255)
                                                                                                                              gACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCg                         >  1:825727/1‑150 (MQ=255)
                                                                                                                              gACCGTAGTCGGTGGCATAAGATGACTCTAACAAAAGCTGGAATACGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATACAATTCATGTCg                         >  2:2900979/1‑150 (MQ=255)
                                                                                                                                 cGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGAtt                      >  1:2469191/1‑150 (MQ=255)
                                                                                                                                  gTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGAttt                     <  2:2209194/150‑1 (MQ=255)
                                                                                                                                       cGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACt                <  1:250927/150‑1 (MQ=255)
                                                                                                                                       cGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACt                <  2:2374938/150‑1 (MQ=255)
                                                                                                                                        ggTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTg               >  2:4190102/1‑150 (MQ=255)
                                                                                                                                         gTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTgg              >  2:3924504/1‑150 (MQ=255)
                                                                                                                                          tGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGc             <  1:502570/150‑1 (MQ=255)
                                                                                                                                          tGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGc             >  1:2925575/1‑150 (MQ=255)
                                                                                                                                            gCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGa           <  1:4019594/150‑1 (MQ=255)
                                                                                                                                             cTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGat          >  1:2705252/1‑150 (MQ=255)
                                                                                                                                                cAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATc       >  1:1492729/1‑150 (MQ=255)
                                                                                                                                                cAGATTACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATc       <  1:4540370/150‑1 (MQ=255)
                                                                                                                                                 aGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCt      <  2:2442401/150‑1 (MQ=255)
                                                                                                                                                     gACTCTAGCAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGAAc  <  2:3826700/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAAATATTATCTTTTTTCGCGGAGGGCGAATGATAAAATCGTTGTTTATCATCAGGAACAAGTAAATGGACCGCTTGTCGACGAGTCAGGGCGGGTTATTTTTGCCGATTTTAACCCATCGAAACGACCGTGGTCGGTGGCTTCAGATGACTCTAACAACAGCTGGAATCCGGCATACAATTGCTTTGATCGTCCGGGTAAAAAATGTATCTCTTTTACGCTACACATCAACGGCAAAGATCACGATTTGTTAGCGGTGGATAAAATTCATGTCGATTTAAACTGGCGATATCTGAAC  >  CP009273/868748‑869046

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