Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 93,391 G→T 100% T301T (ACG→ACT mraY → phospho‑N‑acetylmuramoyl‑pentapeptide transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927393,3910GT97.1% 540.1 / 6.5 174T301T (ACG→ACTmraYphospho‑N‑acetylmuramoyl‑pentapeptide transferase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (76/93);  total (79/95)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 174 total aligned reads displayed.

CCTGTGGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACGCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTG  >  CP009273/93247‑93535
                                                                                                                                                |                                                                                                                                                
ccTGTGGTTTAACACCTATCCGGCGCATGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGAGTGTTCGTGGTAGAAACTCtt                                                                                                                                               <  2:248747/148‑1 (MQ=255)
 cTGTGGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCttt                                                                                                                                              >  2:58327/1‑148 (MQ=255)
  tgtgGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTc                                                                                                                                             <  2:684117/148‑1 (MQ=255)
  tgtgGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTc                                                                                                                                             <  2:645137/148‑1 (MQ=255)
   gtgGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCt                                                                                                                                            >  2:1110268/1‑148 (MQ=255)
   gtgGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCt                                                                                                                                            >  2:2100817/1‑148 (MQ=255)
    tgGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTg                                                                                                                                           >  2:2781156/1‑148 (MQ=255)
      gTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTc                                                                                                                                         >  1:1971826/1‑148 (MQ=255)
      gTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTc                                                                                                                                         >  1:199095/1‑148 (MQ=255)
        ttAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCAt                                                                                                                                       <  1:2595550/148‑1 (MQ=255)
        ttAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCAt                                                                                                                                       <  1:80945/148‑1 (MQ=255)
         tAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATc                                                                                                                                      <  2:2852775/148‑1 (MQ=255)
           acacCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCt                                                                                                                                    >  1:3493842/1‑148 (MQ=255)
             acCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGc                                                                                                                                  <  1:2366879/148‑1 (MQ=255)
             acCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGc                                                                                                                                  >  1:434219/1‑148 (MQ=255)
              ccTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCa                                                                                                                                 <  2:643609/148‑1 (MQ=255)
                tATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTTCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAgg                                                                                                                               <  1:1093657/148‑1 (MQ=255)
                 aTCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGt                                                                                                                              >  2:746824/1‑148 (MQ=255)
                 aTCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGt                                                                                                                              <  2:2029290/148‑1 (MQ=255)
                   ccGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCg                                                                                                                            >  2:2322951/1‑148 (MQ=255)
                    cGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCgg                                                                                                                           <  1:1111785/148‑1 (MQ=255)
                        gcAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTcc                                                                                                                       <  1:1707697/148‑1 (MQ=255)
                         cAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCt                                                                                                                      <  2:391159/148‑1 (MQ=255)
                          aGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCtt                                                                                                                     <  1:1788645/148‑1 (MQ=255)
                          aGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCtt                                                                                                                     <  1:1787065/148‑1 (MQ=255)
                            gTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTa                                                                                                                   >  2:1624192/1‑148 (MQ=255)
                              cTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTaaa                                                                                                                 <  1:2725590/148‑1 (MQ=255)
                               tttATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAAc                                                                                                                <  1:349855/148‑1 (MQ=255)
                                ttATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACt                                                                                                               >  2:1692876/1‑148 (MQ=255)
                                ttATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACt                                                                                                               >  2:1315414/1‑148 (MQ=255)
                                    gggCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTgcgc                                                                                                           >  2:3226235/1‑148 (MQ=255)
                                        gATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGAc                                                                                                       <  1:3481550/148‑1 (MQ=255)
                                        gATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGAc                                                                                                       <  1:3481452/148‑1 (MQ=255)
                                         aTGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACa                                                                                                      <  1:664047/148‑1 (MQ=255)
                                          tGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACaa                                                                                                     <  1:2715832/148‑1 (MQ=255)
                                          tGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACaa                                                                                                     <  1:2207491/148‑1 (MQ=255)
                                          tGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACaa                                                                                                     <  1:3498015/148‑1 (MQ=255)
                                              ggTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACGCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTa                                                                                                 >  1:299734/1‑148 (MQ=255)
                                              ggTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTAGGGGGGGCGTGTTCGTGGTAGAAAATCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTa                                                                                                 >  1:2304005/1‑148 (MQ=255)
                                               gTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTAt                                                                                                >  2:1982947/1‑148 (MQ=255)
                                               gTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTAt                                                                                                >  1:509499/1‑148 (MQ=255)
                                                      ggCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGc                                                                                         <  1:1260877/148‑1 (MQ=255)
                                                        cGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCAt                                                                                       >  1:3410236/1‑148 (MQ=255)
                                                          ttAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATgg                                                                                     <  2:1185746/148‑1 (MQ=255)
                                                           tAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGc                                                                                    >  2:2751530/1‑148 (MQ=255)
                                                              gtggtgCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCAcc                                                                                 >  2:1531471/1‑148 (MQ=255)
                                                                 gtgCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACa                                                                                                      <  2:1950583/124‑1 (MQ=255)
                                                                 gtgCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACa                                                                                                      >  1:1950583/1‑124 (MQ=255)
                                                                 gtgCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCGGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGAt                                                                              <  1:853525/148‑1 (MQ=255)
                                                                  tgCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGAtt                                                                             <  1:1054499/148‑1 (MQ=255)
                                                                   gCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATtc                                                                            <  2:2222089/148‑1 (MQ=255)
                                                                     gTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATtcat                                                                          <  1:1584484/148‑1 (MQ=255)
                                                                           cATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATcaccac                                                                    >  2:215147/1‑148 (MQ=255)
                                                                           cATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCGTTCTGTCATCCTGCATGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGAGTCATcaccac                                                                    >  2:1262639/1‑148 (MQ=255)
                                                                             ttATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTa                                                                  >  1:114177/1‑148 (MQ=255)
                                                                                 cGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGaa                                                              >  2:84651/1‑148 (MQ=255)
                                                                                  gCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAAc                                                             <  2:308155/148‑1 (MQ=255)
                                                                                     gTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCCTCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGa                                                          >  2:1638211/1‑148 (MQ=255)
                                                                                      tACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGaa                                                         <  2:3271351/148‑1 (MQ=255)
                                                                                      tACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGaa                                                         <  1:3441372/148‑1 (MQ=255)
                                                                                      tACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGaa                                                         <  2:291496/148‑1 (MQ=255)
                                                                                      tACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGaa                                                         <  2:2571272/148‑1 (MQ=255)
                                                                                       aCTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGaaa                                                        <  1:2726465/148‑1 (MQ=255)
                                                                                        cTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGTCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAg                                                       <  2:159280/148‑1 (MQ=255)
                                                                                        cTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAg                                                       <  2:2117811/148‑1 (MQ=255)
                                                                                         tGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAgg                                                      <  1:584405/148‑1 (MQ=255)
                                                                                           cTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCt                                                    >  1:3446426/1‑148 (MQ=255)
                                                                                           cTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCt                                                    <  1:242522/148‑1 (MQ=255)
                                                                                             aCGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTgg                                                  >  1:3151133/1‑148 (MQ=255)
                                                                                             aCGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTgg                                                  >  1:667288/1‑148 (MQ=255)
                                                                                                tCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCg                                               >  2:1968652/1‑148 (MQ=255)
                                                                                                  aGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGa                                             >  2:1422898/1‑148 (MQ=255)
                                                                                                        tCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACcgcg                                       <  1:2386657/148‑1 (MQ=255)
                                                                                                           tgctgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGt                                    >  1:484386/1‑148 (MQ=255)
                                                                                                            gctgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTc                                   >  1:1329644/1‑148 (MQ=255)
                                                                                                            gctgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTc                                   >  2:2626934/1‑148 (MQ=255)
                                                                                                             ctgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCa                                  <  1:1220913/148‑1 (MQ=255)
                                                                                                              tgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCAt                                 <  1:990417/148‑1 (MQ=255)
                                                                                                              tgGTGATTATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCAt                                 <  1:990923/148‑1 (MQ=255)
                                                                                                                     tATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCAtt                                >  2:1620180/1‑142 (MQ=255)
                                                                                                                     tATGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCAtt                                <  1:1620180/142‑1 (MQ=255)
                                                                                                                       tGGGGGGCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGttt                        <  1:1673775/148‑1 (MQ=255)
                                                                                                                         gggggCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCt                      >  1:2088417/1‑148 (MQ=255)
                                                                                                                           gggCGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTgg                    <  1:1679579/148‑1 (MQ=255)
                                                                                                                              cGTGTTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGatt                 <  1:2854283/148‑1 (MQ=255)
                                                                                                                                tgtTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGattat               <  2:434219/148‑1 (MQ=255)
                                                                                                                                 gtTCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGattatt              >  2:1783989/1‑148 (MQ=255)
                                                                                                                                  ttCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTAttt             >  2:3333566/1‑148 (MQ=255)
                                                                                                                                   tCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTc            >  1:1391209/1‑148 (MQ=255)
                                                                                                                                   tCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTc            >  2:527019/1‑148 (MQ=255)
                                                                                                                                   tCGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTc            >  1:1528875/1‑148 (MQ=255)
                                                                                                                                    cGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCg           <  2:2465089/148‑1 (MQ=255)
                                                                                                                                    cGTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGACTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCg           <  1:1394524/148‑1 (MQ=255)
                                                                                                                                     gTGGTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGc          <  1:2675416/148‑1 (MQ=255)
                                                                                                                                        gTAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGa       >  1:732211/1‑148 (MQ=255)
                                                                                                                                         tAGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGTCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGAt      >  2:2295666/1‑148 (MQ=255)
                                                                                                                                          aGAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATg     >  2:2283445/1‑148 (MQ=255)
                                                                                                                                           gAAACTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGc    >  2:808193/1‑148 (MQ=255)
                                                                                                                                            aaaCTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGACCCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCt   >  2:376347/1‑148 (MQ=255)
                                                                                                                                            aaaCTCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCt   >  2:375923/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCTGTGGTTTAACACCTATCCGGCGCAGGTCTTTATGGGCGATGTAGGTTCGCTGGCGTTAGGTGGTGCGTTAGGCATTATCGCCGTACTGCTACGTCAGGAATTCCTGCTGGTGATTATGGGGGGCGTGTTCGTGGTAGAAACGCTTTCTGTCATCCTGCAGGTCGGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTG  >  CP009273/93247‑93535

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: