Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,373,880 G→T 9.8% M75I (ATG→ATT)  yjeM → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,373,8800GT 9.8% 1031.3 / 56.9 284M75I (ATG→ATT) yjeMputative transporter
Reads supporting (aligned to +/- strand):  ref base G (139/117);  new base T (22/6);  total (161/123)
Fisher's exact test for biased strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 284 total aligned reads displayed.

ATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGCACTGGCGTA  >  CP009273/4373732‑4374028
                                                                                                                                                    |                                                                                                                                                    
aTAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGa                                                                                                                                                     <  1:2719503/150‑1 (MQ=255)
aTAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGa                                                                                                                                                     <  1:810979/150‑1 (MQ=255)
aTAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCAGTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGa                                                                                                                                                     <  2:2755315/150‑1 (MQ=255)
  aGCCCATCGGCTTATTACTTACTGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGaat                                                                                                                                                   <  1:2308302/150‑1 (MQ=255)
  aGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGaat                                                                                                                                                   >  2:389873/1‑150 (MQ=255)
  aGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGaat                                                                                                                                                   >  2:2629950/1‑150 (MQ=255)
   gCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGaata                                                                                                                                                  <  2:2571521/150‑1 (MQ=255)
      cATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAataata                                                                                                                                               <  1:2729455/150‑1 (MQ=255)
      cATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAataata                                                                                                                                               <  1:2211408/150‑1 (MQ=255)
       aTCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAg                                                                                                                                              >  2:3588450/1‑150 (MQ=255)
        tCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAgt                                                                                                                                             <  1:3603904/150‑1 (MQ=255)
        tCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAgt                                                                                                                                             <  1:96837/150‑1 (MQ=255)
         cGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTGTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAgtg                                                                                                                                            >  2:3857112/1‑150 (MQ=255)
          ggCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAgtgt                                                                                                                                           >  2:705586/1‑150 (MQ=255)
            cTTATTACTTTATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCATTAATGATGGCTGAATTGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCg                                                                                                                                         <  1:610106/150‑1 (MQ=255)
             ttattaCTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCgg                                                                                                                                        >  1:542005/1‑150 (MQ=255)
              tattaCTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGa                                                                                                                                       <  1:4354511/150‑1 (MQ=255)
                 taCTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCa                                                                                                                                    <  1:1979467/150‑1 (MQ=255)
                  aCTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCAc                                                                                                                                   <  2:1074213/150‑1 (MQ=255)
                   cTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACg                                                                                                                                  >  1:4268489/1‑150 (MQ=255)
                      aaTGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGttt                                                                                                                               >  2:4647/1‑150 (MQ=255)
                      aaTGGGTTAGAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGttt                                                                                                                               <  2:4311652/150‑1 (MQ=255)
                        tGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTg                                                                                                                             <  2:677195/150‑1 (MQ=255)
                           gTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCt                                                                                                                          >  2:3784072/1‑150 (MQ=255)
                               taGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCAt                                                                                                                      >  1:4560898/1‑150 (MQ=255)
                               taGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCAt                                                                                                                      >  1:2156153/1‑150 (MQ=255)
                                   gCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGt                                                                                                                  >  1:1016307/1‑150 (MQ=255)
                                    cGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTa                                                                                                                 >  1:1558740/1‑150 (MQ=255)
                                     gATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTAc                                                                                                                <  1:1407152/150‑1 (MQ=255)
                                      aTTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACg                                                                                                               >  2:1988411/1‑150 (MQ=255)
                                        tCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGtt                                                                                                             >  2:77703/1‑150 (MQ=255)
                                         cccTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGttt                                                                                                            >  2:3062909/1‑150 (MQ=255)
                                         cccTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGttt                                                                                                            >  1:363282/1‑150 (MQ=255)
                                         cccTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGttt                                                                                                            <  1:3776764/150‑1 (MQ=255)
                                         cccTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGttt                                                                                                            >  2:1923054/1‑150 (MQ=255)
                                          ccTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTa                                                                                                           <  2:1692627/150‑1 (MQ=255)
                                             tttATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTAtgt                                                                                                        >  1:1878474/1‑150 (MQ=255)
                                                atatatTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGt                                                                                                     >  1:1103967/1‑150 (MQ=255)
                                                  atatTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCATTCATTGGTACGTTTATGTGGttt                                                                                                   >  1:3597427/1‑150 (MQ=255)
                                                    atTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTc                                                                                                 <  1:1295401/150‑1 (MQ=255)
                                                    atTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTc                                                                                                 >  2:2395723/1‑150 (MQ=255)
                                                    atTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCAGTGGTACGTGTATGTGGTTTTc                                                                                                 >  1:3630018/1‑150 (MQ=255)
                                                      tttCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGTTTTTCct                                                                                               >  2:808417/1‑150 (MQ=255)
                                                      tttCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCct                                                                                               >  1:3392326/1‑150 (MQ=255)
                                                      tttCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCct                                                                                               <  1:2506237/150‑1 (MQ=255)
                                                       ttCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCctc                                                                                              >  1:3351324/1‑150 (MQ=255)
                                                        tCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCctct                                                                                             <  2:1057331/150‑1 (MQ=255)
                                                        tCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCAGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCctct                                                                                             <  1:960603/150‑1 (MQ=255)
                                                          tGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTta                                                                                           <  2:235667/150‑1 (MQ=255)
                                                             aTTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTatat                                                                                        <  2:2959992/150‑1 (MQ=255)
                                                               tGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATatca                                                                                      <  1:1553276/150‑1 (MQ=255)
                                                                 ttATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATatcatc                                                                                    >  2:4090259/1‑150 (MQ=255)
                                                                  tATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCt                                                                                   <  2:4546402/150‑1 (MQ=255)
                                                                  tATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCt                                                                                   <  2:1499412/150‑1 (MQ=255)
                                                                    ttCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTgg                                                                                 >  2:4425285/1‑150 (MQ=255)
                                                                      cTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGAt                                                                               >  1:739071/1‑150 (MQ=255)
                                                                       tttATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATg                                                                              <  2:2875280/150‑1 (MQ=255)
                                                                         tATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGt                                                                            >  2:419511/1‑150 (MQ=255)
                                                                           ttCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGGTTTGCCTTCAGTGGTCCGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGa                                                                          >  2:4222862/1‑150 (MQ=255)
                                                                             ccATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTAGTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGt                                                                        <  2:1761143/150‑1 (MQ=255)
                                                                              cATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTa                                                                       <  1:724941/150‑1 (MQ=255)
                                                                               aTTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTAc                                                                      >  2:1488438/1‑150 (MQ=255)
                                                                                  cGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTc                                                                   >  1:1266371/1‑150 (MQ=255)
                                                                                   gCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTcc                                                                  >  2:2019413/1‑150 (MQ=255)
                                                                                   gCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTcc                                                                  >  2:4005778/1‑150 (MQ=255)
                                                                                        aaTGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGaa                                                             >  1:455692/1‑150 (MQ=255)
                                                                                         atgatgGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGaaa                                                            >  2:3630839/1‑150 (MQ=255)
                                                                                           gatgGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGt                                                          >  2:2143160/1‑150 (MQ=255)
                                                                                                 tGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGt                                                    <  2:479396/150‑1 (MQ=255)
                                                                                                   aaaTGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTAc                                                  >  1:1966385/1‑150 (MQ=255)
                                                                                                    aaTGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCGTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTAcc                                                 >  1:2091521/1‑150 (MQ=255)
                                                                                                       gggAGCTGCTTATCGCAAAGAAGAAGGCGGTCTCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGtt                                              >  2:56527/1‑150 (MQ=255)
                                                                                                         gAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTtct                                            <  2:2947854/150‑1 (MQ=255)
                                                                                                         gAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTtct                                            >  1:3954368/1‑150 (MQ=255)
                                                                                                          aGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAACGTTTGGGTACCGTtctc                                           >  1:1330185/1‑150 (MQ=255)
                                                                                                          aGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATAACTTGATGAATAATAGTGTAGTACCACTTTTTTCCTTCATTGGTACGTTTATGTTGTTTTACTCTTATATCATCAGGATGGTGAGTACCTACGCGAAAGTTTGGGTACCTAtctc                                           >  2:2649784/1‑150 (MQ=255)
                                                                                                           gctgctTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCa                                          >  2:3717035/1‑150 (MQ=255)
                                                                                                           gctgctTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCa                                          >  2:3176309/1‑150 (MQ=255)
                                                                                                              gctTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACa                                       >  2:3172007/1‑150 (MQ=255)
                                                                                                                 tATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTc                                    <  2:1433419/150‑1 (MQ=255)
                                                                                                                 tATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTc                                    >  1:3875194/1‑150 (MQ=255)
                                                                                                                    cGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCctc                                 >  1:3711789/1‑150 (MQ=255)
                                                                                                                     gCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCctct                                <  1:3929081/150‑1 (MQ=255)
                                                                                                                        aagaagaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATg                             >  2:3319344/1‑150 (MQ=255)
                                                                                                                         agaagaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATgg                            >  2:836073/1‑150 (MQ=255)
                                                                                                                         agaagaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATgg                            >  2:1136635/1‑150 (MQ=255)
                                                                                                                          gaagaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGt                           >  1:2547630/1‑150 (MQ=255)
                                                                                                                           aagaagGCGGTATCTATTCCTGGATTAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTa                          >  1:356360/1‑150 (MQ=255)
                                                                                                                           aagaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTa                          <  2:1611162/150‑1 (MQ=255)
                                                                                                                           aagaagGCGGTATCTATTCCTGGATGAATAATAGTGTAGGACCACGTTATGCCTTCATTGGTACGATTATGTGGTATTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACAGGTCTCAACATTCCTAAATGGTa                          >  2:1981189/1‑150 (MQ=255)
                                                                                                                             gaagGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGc                        >  1:4094885/1‑150 (MQ=255)
                                                                                                                               agGCGGTATCTATTCCTGGATGACTAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCGTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGa                      >  2:2263403/1‑150 (MQ=255)
                                                                                                                                ggCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGAc                     <  2:2398185/150‑1 (MQ=255)
                                                                                                                                      aTCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGAcc               <  2:3011880/150‑1 (MQ=255)
                                                                                                                                      aTCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGAcc               <  2:1681226/150‑1 (MQ=255)
                                                                                                                                          aTTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGc           >  1:3002284/1‑150 (MQ=255)
                                                                                                                                          aTTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGc           <  2:2089345/150‑1 (MQ=255)
                                                                                                                                              cTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGCACTg       >  2:3402360/1‑150 (MQ=255)
                                                                                                                                               tGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGCACTgg      >  1:234301/1‑150 (MQ=255)
                                                                                                                                                 gATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGCACTGGCg    <  2:1283290/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTTATTCTTTATTCCATTCGCCTTAATGATGGCTGAAATGGGAGCTGCTTATCGCAAAGAAGAAGGCGGTATCTATTCCTGGATGAATAATAGTGTCGGACCACGTTTTGCCTTCATTGGTACGTTTATGTGGTTTTCCTCTTATATCATCTGGATGGTGAGTACCTCCGCGAAAGTTTGGGTACCGTTCTCAACATTCCTCTATGGTAGCGACATGACCCAGCACTGGCGTA  >  CP009273/4373732‑4374028

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: ‑