Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,154,402 A→T 40.9% E279V (GAA→GTA)  btuB → vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,154,4020AT40.9% 162.7 / 309.2 232E279V (GAA→GTA) btuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (49/88);  new base T (38/57);  total (87/145)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 236 total aligned reads displayed.

TACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCGTTGGTCA  >  CP009273/4154256‑4154549
                                                                                                                                                  |                                                                                                                                                   
tACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGaa                                                                                                                                                    <  1:591047/148‑1 (MQ=255)
   tGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTg                                                                                                                                                 <  2:2878841/148‑1 (MQ=255)
    gATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGa                                                                                                                                                <  2:3658326/148‑1 (MQ=255)
      tGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGAtt                                                                                                                                              <  1:2187293/148‑1 (MQ=255)
          tGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAAt                                                                                                                                          <  2:3384357/148‑1 (MQ=255)
           ggAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATc                                                                                                                                         >  2:1570473/1‑148 (MQ=255)
              gCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATcaca                                                                                                                                      >  1:173328/1‑148 (MQ=255)
              gCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATcaca                                                                                                                                      <  1:621855/148‑1 (MQ=255)
               cGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACaa                                                                                                                                     >  2:2631045/1‑148 (MQ=255)
               cGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACaa                                                                                                                                     <  2:591252/148‑1 (MQ=255)
                 gCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACt                                                                                                                                   >  1:1495048/1‑148 (MQ=255)
                 gCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACt                                                                                                                                   >  1:2700689/1‑148 (MQ=255)
                   tttGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCa                                                                                                                                 <  1:2999877/148‑1 (MQ=255)
                        gcgcGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTAcc                                                                                                                            <  2:1212986/148‑1 (MQ=255)
                        gcgcGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTAcc                                                                                                                            <  2:788711/148‑1 (MQ=255)
                        gcgcGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTAcc                                                                                                                            <  1:1892830/148‑1 (MQ=255)
                        gcgcGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTAcc                                                                                                                            <  2:2705344/148‑1 (MQ=255)
                          gcgGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAg                                                                                                                          <  1:3387708/148‑1 (MQ=255)
                           cgGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGc                                                                                                                         >  1:51623/1‑148 (MQ=255)
                              ctatggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCtat                                                                                                                      <  1:2935323/148‑1 (MQ=255)
                              ctatggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCtat                                                                                                                      <  1:3672936/148‑1 (MQ=255)
                               tatggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCtata                                                                                                                     >  2:1358212/1‑148 (MQ=255)
                               tatggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCtata                                                                                                                     >  2:3444069/1‑148 (MQ=255)
                                atggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAg                                                                                                                    >  1:3761931/1‑148 (MQ=255)
                                 tggctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGc                                                                                                                   <  1:247337/148‑1 (MQ=255)
                                   gctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCa                                                                                                                 <  1:864500/148‑1 (MQ=255)
                                    ctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCAt                                                                                                                <  2:3063783/148‑1 (MQ=255)
                                    ctatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCAt                                                                                                                >  2:2077093/1‑148 (MQ=255)
                                     tatgATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATa                                                                                                               <  1:1624904/148‑1 (MQ=255)
                                         aTAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCaa                                                                                                           >  1:602312/1‑148 (MQ=255)
                                             ccGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGAt                                                                                                       >  2:853390/1‑148 (MQ=255)
                                             ccGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGAt                                                                                                       <  1:1378127/148‑1 (MQ=255)
                                              cGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGAtt                                                                                                      >  2:1279475/1‑148 (MQ=255)
                                               gTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTa                                                                                                     <  1:2346620/148‑1 (MQ=255)
                                                  ccAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGATGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACaa                                                                                                  >  2:1018195/1‑148 (MQ=255)
                                                   cAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAAc                                                                                                 <  1:1172788/148‑1 (MQ=255)
                                                    aaTTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACt                                                                                                <  2:1976525/148‑1 (MQ=255)
                                                     aTTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTa                                                                                               <  1:3092828/148‑1 (MQ=255)
                                                     aTTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTa                                                                                               <  2:1774791/148‑1 (MQ=255)
                                                     aTTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTa                                                                                               <  2:3303859/148‑1 (MQ=255)
                                                     aTTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTa                                                                                               <  2:2612452/148‑1 (MQ=255)
                                                        aTGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGa                                                                                            <  2:3812640/148‑1 (MQ=255)
                                                        aTGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGa                                                                                            <  2:3739332/148‑1 (MQ=255)
                                                             gcgTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATcccc                                                                                       <  1:1713698/148‑1 (MQ=255)
                                                             gcgTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATcccc                                                                                       <  1:1511455/148‑1 (MQ=255)
                                                             gcgTATTAGTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATcccc                                                                                       <  1:2889072/148‑1 (MQ=255)
                                                                tattatTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCAtt                                                                                    <  1:3683959/148‑1 (MQ=255)
                                                                tattatTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAg                                                                                                         >  2:324890/1‑127 (MQ=255)
                                                                tattatTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCAtt                                                                                    <  2:1251695/148‑1 (MQ=255)
                                                                 attattCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTa                                                                                   >  1:401737/1‑148 (MQ=255)
                                                                 attattCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTa                                                                                   >  1:2304084/1‑148 (MQ=255)
                                                                   tattCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATg                                                                                 <  1:3125239/148‑1 (MQ=255)
                                                                   tattCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCAGTATg                                                                                 >  1:3180187/1‑148 (MQ=255)
                                                                   tattCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATg                                                                                 <  1:3316485/148‑1 (MQ=255)
                                                                     ttCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGt                                                                               <  1:707482/148‑1 (MQ=255)
                                                                       ctcCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCg                                                                             <  2:513382/148‑1 (MQ=255)
                                                                       ctcCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCg                                                                             <  2:2849881/148‑1 (MQ=255)
                                                                            tgTTCACCGTTGCTCGATACCAGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATg                                                                        <  2:2756760/147‑1 (MQ=255)
                                                                            ggTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATg                                                                        <  2:1560589/148‑1 (MQ=255)
                                                                            ggTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATg                                                                        <  1:3290220/148‑1 (MQ=255)
                                                                              ttCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGAt                                                                      <  1:1863968/148‑1 (MQ=255)
                                                                              ttCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGAt                                                                      <  2:2776566/148‑1 (MQ=255)
                                                                               tCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGAtt                                                                     <  1:1661262/148‑1 (MQ=255)
                                                                               tCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGAtt                                                                     >  2:661367/1‑148 (MQ=255)
                                                                                 aCCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATtcg                                                                   >  1:2251266/1‑148 (MQ=255)
                                                                                 aCCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATtcg                                                                   >  1:2055935/1‑148 (MQ=255)
                                                                                 aCCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATtcg                                                                   >  1:2022131/1‑148 (MQ=255)
                                                                                  ccGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATtcgt                                                                  <  1:7148/148‑1 (MQ=255)
                                                                                         tCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGAc                                                           >  1:1405034/1‑148 (MQ=255)
                                                                                            aTACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCt                                                        <  2:2455130/148‑1 (MQ=255)
                                                                                              aCCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCg                                                      >  2:1748586/1‑148 (MQ=255)
                                                                                                 cGTAATCTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCgatg                                                   <  2:1844452/148‑1 (MQ=255)
                                                                                                   tAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCgatgag                                                 <  2:1460026/148‑1 (MQ=255)
                                                                                                     aaCTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCgatgagat                                               >  2:3457966/1‑148 (MQ=255)
                                                                                                       ctctATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCgatgagatga                                             <  1:1071368/148‑1 (MQ=255)
                                                                                                         ctATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAg                                           <  1:3483503/148‑1 (MQ=255)
                                                                                                           ataGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCa                                         >  2:2234922/1‑148 (MQ=255)
                                                                                                           ataGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCa                                         <  1:1254387/148‑1 (MQ=255)
                                                                                                               ccAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATac                                     >  1:1952799/1‑148 (MQ=255)
                                                                                                                cAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATaca                                    <  2:2423287/148‑1 (MQ=255)
                                                                                                                cAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATaca                                    >  1:624451/1‑148 (MQ=255)
                                                                                                                cAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATaca                                    <  1:1570473/148‑1 (MQ=255)
                                                                                                                 aaaGTTGGGACGCCGGGATGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATacac                                   <  1:1018195/148‑1 (MQ=255)
                                                                                                                      tGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTcc                              >  2:942107/1‑148 (MQ=255)
                                                                                                                       gggACGCCGGGCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCa                             <  2:1024444/148‑1 (MQ=255)
                                                                                                                        ggACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAg                            <  1:2352380/148‑1 (MQ=255)
                                                                                                                               gggCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCaa                     <  1:2820367/148‑1 (MQ=255)
                                                                                                                               gggCTGCGCTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCaa                     <  1:3802477/148‑1 (MQ=255)
                                                                                                                                    gcgcTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAAt                >  1:1217979/1‑148 (MQ=255)
                                                                                                                                    gcgcTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAAt                <  2:2300776/148‑1 (MQ=255)
                                                                                                                                     cgcTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATg               <  2:3611446/148‑1 (MQ=255)
                                                                                                                                       cTATAACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTc             <  2:615176/148‑1 (MQ=255)
                                                                                                                                        tataACGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCa            >  2:3015506/1‑148 (MQ=255)
                                                                                                                                        tataACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCa            <  2:602166/148‑1 (MQ=255)
                                                                                                                                         ataACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCAt           >  2:3763811/1‑148 (MQ=255)
                                                                                                                                         ataACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCAt           <  1:1814986/148‑1 (MQ=255)
                                                                                                                                           aaCGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCg         >  2:766987/1‑148 (MQ=255)
                                                                                                                                            aCGGCGTACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCGt        <  1:3106341/148‑1 (MQ=255)
                                                                                                                                                cGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCGTTGGt    >  1:3611190/1‑148 (MQ=255)
                                                                                                                                                  aaCTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCGTTGGTCa  <  1:992909/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATGATAACCGTACCAATTATGACGCGTATTATTCTCCCGGTTCACCGTTGCTCGATACCCGTAAACTCTATAGCCAAAGTTGGGACGCCGGGCTGCGCTATAACGGCGAACTGATTAAATCACAACTCATTACCAGCTATAGCCATAGCAAAGATTACAACTACGATCCCCATTATGGTCGTTATGATTCGTCGGCGACGCTCGATGAGATGAAGCAATACACCGTCCAGTGGGCAAACAATGTCATCGTTGGTCA  >  CP009273/4154256‑4154549

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: