Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,521,598 C→A 14.5% Q591K (CAG→AAG)  yrfF → inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,521,5980CA14.5% 963.4 / 94.3 290Q591K (CAG→AAG) yrfFinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (114/134);  new base A (16/26);  total (130/160)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 300 total aligned reads displayed.

CGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCAT  >  CP009273/3521449‑3521747
                                                                                                                                                     |                                                                                                                                                     
cGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATc                                                                                                                                                       <  2:731970/150‑1 (MQ=255)
     aGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGccc                                                                                                                                                  <  2:3825823/150‑1 (MQ=255)
     aGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGccc                                                                                                                                                  <  2:42524/150‑1 (MQ=255)
      gTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCt                                                                                                                                                 <  2:1613168/150‑1 (MQ=255)
      gTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCt                                                                                                                                                 <  2:1273590/150‑1 (MQ=255)
       tGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTg                                                                                                                                                >  1:4041495/1‑150 (MQ=255)
       tGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTg                                                                                                                                                >  2:1917558/1‑150 (MQ=255)
       tGCGGAATCGCTGGATAACCTGGTGGAAACCTCAACCGCGCCGTTTATCACGCATGAAAAGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGAGTCCTGAGTGTCAGCGACGAAGGCAGCGATTTTGATGATCAGCCCTg                                                                                                                                                >  2:2059294/1‑150 (MQ=255)
         cGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGc                                                                                                                                              >  1:2347814/1‑150 (MQ=255)
         cGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGc                                                                                                                                              <  1:1926141/150‑1 (MQ=255)
          ggAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGcc                                                                                                                                             <  2:2929228/150‑1 (MQ=255)
          ggAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGcc                                                                                                                                             >  1:449017/1‑150 (MQ=255)
           gAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCt                                                                                                                                            <  2:490987/150‑1 (MQ=255)
                  tGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCAt                                                                                                                                     >  1:1365037/1‑150 (MQ=255)
                   ggATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATc                                                                                                                                    >  1:2760902/1‑150 (MQ=255)
                     aTAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCAc                                                                                                                                  <  2:775088/150‑1 (MQ=255)
                      tAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACt                                                                                                                                 <  2:1457881/150‑1 (MQ=255)
                       aaCCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACtt                                                                                                                                >  2:3610678/1‑150 (MQ=255)
                          cTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTAc                                                                                                                             <  2:1218716/150‑1 (MQ=255)
                           tggtggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACg                                                                                                                            <  1:2441958/150‑1 (MQ=255)
                           tggtggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACg                                                                                                                            >  1:3080213/1‑150 (MQ=255)
                              tggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACt                                                                                                                         <  1:2979389/150‑1 (MQ=255)
                              tggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACt                                                                                                                         <  1:3674932/150‑1 (MQ=255)
                               ggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTa                                                                                                                        <  1:385169/150‑1 (MQ=255)
                               ggCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTa                                                                                                                        >  1:471680/1‑150 (MQ=255)
                                gCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTAc                                                                                                                       >  2:4562837/1‑150 (MQ=255)
                                  aaCCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTAccc                                                                                                                     <  2:1102614/150‑1 (MQ=255)
                                     cTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACccgcc                                                                                                                  <  2:1852675/150‑1 (MQ=255)
                                        cACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCa                                                                                                               <  2:2278833/150‑1 (MQ=255)
                                         aCCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCaa                                                                                                              >  1:3319333/1‑150 (MQ=255)
                                          ccGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAg                                                                                                             <  2:3220026/150‑1 (MQ=255)
                                           cgcgCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGa                                                                                                            >  1:1246456/1‑150 (MQ=255)
                                           cgcgCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGa                                                                                                            >  1:2827177/1‑150 (MQ=255)
                                             cgcCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAAc                                                                                                          >  2:4270450/1‑150 (MQ=255)
                                              gcCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACa                                                                                                         <  2:878214/150‑1 (MQ=255)
                                              gcCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACa                                                                                                         <  1:3607993/150‑1 (MQ=255)
                                               ccGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAg                                                                                                        >  1:1713192/1‑150 (MQ=255)
                                                 gTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTg                                                                                                      <  2:3714889/150‑1 (MQ=255)
                                                     aTCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAAc                                                                                                  <  1:4140771/150‑1 (MQ=255)
                                                     aTCACGCATGAAACGGCGCGAGCGACACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAAc                                                                                                  <  2:104172/150‑1 (MQ=255)
                                                          gCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCttt                                                                                             >  2:241557/1‑150 (MQ=255)
                                                          gCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCttt                                                                                             <  1:4016452/150‑1 (MQ=255)
                                                           cATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTc                                                                                            <  2:3016673/150‑1 (MQ=255)
                                                            aTGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTTATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTGACGACTACCCGCCGCAAGAACAGTGGAACGCTTTcc                                                                                           <  1:918535/150‑1 (MQ=255)
                                                              gAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAg                                                                                         <  2:4389695/150‑1 (MQ=255)
                                                               aaaCGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGa                                                                                        >  1:1457642/1‑150 (MQ=255)
                                                                  cGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAAc                                                                                     <  2:1340238/150‑1 (MQ=255)
                                                                   ggCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACt                                                                                    <  1:1315716/150‑1 (MQ=255)
                                                                   ggCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACt                                                                                    <  2:2366769/150‑1 (MQ=255)
                                                                        gAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGcaa                                                                               >  2:2749241/1‑149 (MQ=255)
                                                                             gCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATg                                                                          <  2:4042341/150‑1 (MQ=255)
                                                                               acaATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCt                                                                        <  1:2462390/150‑1 (MQ=255)
                                                                                caATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTg                                                                       >  2:1438623/1‑150 (MQ=255)
                                                                                caATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTg                                                                       <  1:3903442/150‑1 (MQ=255)
                                                                                       tAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATa                                                                >  1:3314242/1‑150 (MQ=255)
                                                                                         aaCAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACg                                                              <  2:3610144/150‑1 (MQ=255)
                                                                                           cAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGcc                                                            <  2:1652971/150‑1 (MQ=255)
                                                                                             gTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGt                                                          >  2:1164315/1‑150 (MQ=255)
                                                                                                cGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                       >  1:3505446/1‑150 (MQ=255)
                                                                                                cGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                       <  2:1075332/150‑1 (MQ=255)
                                                                                                cGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTa                                                       >  1:2183561/1‑150 (MQ=255)
                                                                                                 ggCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTaa                                                      <  1:2574981/150‑1 (MQ=255)
                                                                                                  gCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAc                                                     <  1:2736057/150‑1 (MQ=255)
                                                                                                   cccccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACg                                                    >  1:328080/1‑150 (MQ=255)
                                                                                                   cccccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACg                                                    <  1:1046625/150‑1 (MQ=255)
                                                                                                     cccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGcc                                                  <  2:1716857/150‑1 (MQ=255)
                                                                                                     cccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGcc                                                  <  2:1513855/150‑1 (MQ=255)
                                                                                                      ccGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCg                                                 <  2:3287762/150‑1 (MQ=255)
                                                                                                       cggcggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGa                                                <  1:2307557/150‑1 (MQ=255)
                                                                                                           ggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGt                                            >  2:2446132/1‑150 (MQ=255)
                                                                                                           ggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGt                                            <  2:3968339/150‑1 (MQ=255)
                                                                                                           ggATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGt                                            <  2:1764976/150‑1 (MQ=255)
                                                                                                             aTTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCAACAATTTACGACTACCCGCCGCAAGAACAGTGTAACGCTTTCCAGAAACTGGCACAAATGCAGATGCATACGCCGTTAAACCCCGAGGGTAt                                          >  1:1125316/1‑150 (MQ=255)
                                                                                                              ttCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATc                                         <  2:4133233/150‑1 (MQ=255)
                                                                                                              ttCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATc                                         <  2:1121833/150‑1 (MQ=255)
                                                                                                                ccTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGt                                       <  1:2665895/150‑1 (MQ=255)
                                                                                                                  tGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTca                                     >  2:4176737/1‑150 (MQ=255)
                                                                                                                  tGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTca                                     <  1:658476/150‑1 (MQ=255)
                                                                                                                   gATTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTcac                                    >  1:227116/1‑150 (MQ=255)
                                                                                                                    aTTGTCAGCGACGAAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTcaca                                   <  2:178967/150‑1 (MQ=255)
                                                                                                                             gacgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTc                          <  2:3351200/150‑1 (MQ=255)
                                                                                                                             gacgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTc                          >  1:453877/1‑150 (MQ=255)
                                                                                                                              acgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCa                         <  1:466396/150‑1 (MQ=255)
                                                                                                                              acgaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCa                         >  1:2530975/1‑150 (MQ=255)
                                                                                                                                gaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACt                       <  1:3388128/150‑1 (MQ=255)
                                                                                                                                gaAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACt                       >  1:4216979/1‑150 (MQ=255)
                                                                                                                                gaAGGCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACt                       <  2:1373820/150‑1 (MQ=255)
                                                                                                                                    gCAGCGATTTTGTTGATAAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAAAGCCGAGGGTATCGTCACAAAAATCTTCACTGACg                   <  2:399775/150‑1 (MQ=255)
                                                                                                                                     cAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGc                  >  1:540912/1‑150 (MQ=255)
                                                                                                                                      ggcgaggttgttgATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGTACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAATAATCTGCGCTGACGcc                 >  1:1908711/8‑150 (MQ=255)
                                                                                                                                       gCGATTTTGTTGATAAGCCCTGGCCTTCGGCAACACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGATTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCa                <  2:555249/150‑1 (MQ=255)
                                                                                                                                        cGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCaa               >  1:610854/1‑150 (MQ=255)
                                                                                                                                        cGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGAAGCATACGCCGTATAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCaa               >  2:3471934/1‑150 (MQ=255)
                                                                                                                                         gATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAAt              <  2:3186402/150‑1 (MQ=255)
                                                                                                                                          aTTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATg             >  2:2820248/1‑150 (MQ=255)
                                                                                                                                            tttGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGt           <  2:211498/150‑1 (MQ=255)
                                                                                                                                                  gATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACgcag     <  2:2513879/150‑1 (MQ=255)
                                                                                                                                                  gATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACgcag     <  1:2956338/150‑1 (MQ=255)
                                                                                                                                                     cAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCat  >  2:590254/1‑150 (MQ=255)
                                                                                                                                                     cAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCat  <  1:2059273/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTGAGTGCGGAATCGCTGGATAACCTGGTGGCAACCTCCACCGCGCCGTTTATCACGCATGAAACGGCGCGAGCGGCACAATCACTAAACAGTCCGGCCCCCGGCGGATTCCTGATTGTCAGCGACGAAGGCAGCGATTTTGTTGATCAGCCCTGGCCTTCGGCATCACTTTACGACTACCCGCCGCAAGAACAGTGGAACGCTTTCCAGAAACTGGCACAAATGCTGATGCATACGCCGTTTAACGCCGAGGGTATCGTCACAAAAATCTTCACTGACGCCAATGGTACGCAGCAT  >  CP009273/3521449‑3521747

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: