Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,434,260 C→A 92.4% P111H (CCT→CAT)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,434,2600CA92.4% 438.4 / 27.0 158P111H (CCT→CAT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base A (89/57);  total (93/65)
Fisher's exact test for biased strand distribution p-value = 7.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 160 total aligned reads displayed.

TTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCT  >  CP009273/4434116‑4434403
                                                                                                                                                |                                                                                                                                               
ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                              >  2:1562504/1‑148 (MQ=255)
  cgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCtt                                                                                                                                            <  2:3048373/148‑1 (MQ=255)
          gcagcCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAg                                                                                                                                    <  2:1803306/148‑1 (MQ=255)
              ccAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGAAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTg                                                                                                                                <  2:1673553/148‑1 (MQ=255)
              ccAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGATTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTg                                                                                                                                <  2:1673459/148‑1 (MQ=255)
                aGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaa                                                                                                                              >  1:1401055/1‑148 (MQ=255)
                 ggCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaag                                                                                                                             >  2:1012500/1‑148 (MQ=255)
                  gCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaaga                                                                                                                            >  1:1567947/1‑148 (MQ=255)
                  gCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaaga                                                                                                                            >  1:2721808/1‑148 (MQ=255)
                   cGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaagaa                                                                                                                           >  1:118448/1‑148 (MQ=255)
                   cGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaagaa                                                                                                                           >  1:1794549/1‑148 (MQ=255)
                      ggCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGATTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTAAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAagag                                                                                                                        <  1:1914270/148‑1 (MQ=255)
                      ggCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAagag                                                                                                                        <  1:1914304/148‑1 (MQ=255)
                          gTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAg                                                                                                                    >  1:1053355/1‑148 (MQ=255)
                            aaaaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAgaagatgaagat                                                                                                                  >  2:1769695/1‑148 (MQ=255)
                            aaaaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGAt                                                                                                                  >  2:94627/1‑148 (MQ=255)
                              aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAg                                                                                                                >  1:727424/1‑148 (MQ=255)
                              aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAg                                                                                                                <  2:1091608/148‑1 (MQ=255)
                              aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAg                                                                                                                <  1:118486/148‑1 (MQ=255)
                                gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                              >  2:1588385/1‑148 (MQ=255)
                                gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                              >  1:686652/1‑148 (MQ=255)
                                gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                              >  1:686641/1‑148 (MQ=255)
                                gCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCg                                                                                                              <  1:1689349/148‑1 (MQ=255)
                                 cGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCgg                                                                                                             >  2:1954675/1‑148 (MQ=255)
                                 cGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCgg                                                                                                             >  2:3040353/1‑148 (MQ=255)
                                  gggCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGt                                                                                                            <  2:118448/148‑1 (MQ=255)
                                  gggCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGt                                                                                                            <  2:2425029/148‑1 (MQ=255)
                                    gCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTcc                                                                                                          >  2:78741/1‑148 (MQ=255)
                                     cAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCg                                                                                                         >  1:1565149/1‑148 (MQ=255)
                                               ctggctggCGGGCTTTAGTGCCTGTGGAACAGCAGTTTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtctct                                                                                               <  1:2826881/148‑1 (MQ=255)
                                                 ggCTGTCGGGCTGGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTcc                                                                                             <  2:148570/148‑1 (MQ=255)
                                                      tCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGcc                                                                                        >  2:1956795/1‑148 (MQ=255)
                                                          gCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGGCATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCAGCTGGATCTCTCCACGCCGTAt                                                                                    >  1:3098592/1‑148 (MQ=255)
                                                          gCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTAt                                                                                    >  2:360879/1‑148 (MQ=255)
                                                           gttGAGTGCCTGTGTAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGAAATTCAGGTAAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATc                                                                                   <  1:882947/147‑1 (MQ=255)
                                                           cTTGAGTGCCTGTGGAACAGCAGTTTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATc                                                                                   <  1:883184/148‑1 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               >  2:315107/1‑148 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               <  1:2062472/148‑1 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               >  1:3196227/1‑148 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               <  1:2607581/148‑1 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               >  2:1060001/1‑148 (MQ=255)
                                                               aGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCAt                                                                               >  2:361129/1‑148 (MQ=255)
                                                                 tGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCa                                                                             >  1:2126832/1‑148 (MQ=255)
                                                                  gCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAc                                                                            >  1:793939/1‑148 (MQ=255)
                                                                   ccTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAcc                                                                           <  1:1052647/148‑1 (MQ=255)
                                                                   ccTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAcc                                                                           <  1:2469885/148‑1 (MQ=255)
                                                                    cTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAccc                                                                          <  2:2997121/148‑1 (MQ=255)
                                                                     tgtgGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCt                                                                         >  2:1647531/1‑148 (MQ=255)
                                                                      gtgGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTg                                                                        >  1:703359/1‑148 (MQ=255)
                                                                            cagcagTGTTGGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACgg                                                                  <  1:1720604/148‑1 (MQ=255)
                                                                              gcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGt                                                                <  2:86024/148‑1 (MQ=255)
                                                                              gcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGt                                                                >  2:1226231/1‑148 (MQ=255)
                                                                                  tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCa                                                            <  1:80539/148‑1 (MQ=255)
                                                                                  tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCa                                                            >  1:56166/1‑148 (MQ=255)
                                                                                  tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTGGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCa                                                            >  2:2485579/1‑148 (MQ=255)
                                                                                    tttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACt                                                          >  2:574643/1‑148 (MQ=255)
                                                                                    tttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACt                                                          >  1:330455/1‑148 (MQ=255)
                                                                                    tttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACt                                                          <  2:2721808/148‑1 (MQ=255)
                                                                                     ttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTg                                                         >  1:2167637/1‑148 (MQ=255)
                                                                                     ttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTg                                                         >  1:3003828/1‑148 (MQ=255)
                                                                                     ttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTg                                                         >  1:2599783/1‑148 (MQ=255)
                                                                                           tAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAAc                                                   >  1:2940398/1‑148 (MQ=255)
                                                                                               gATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCa                                               <  2:540570/148‑1 (MQ=255)
                                                                                                      cgcTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTaa                                        <  1:2703379/148‑1 (MQ=255)
                                                                                                       gcTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAg                                       <  2:1565585/148‑1 (MQ=255)
                                                                                                       gcTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAg                                       <  2:2386561/148‑1 (MQ=255)
                                                                                                        cTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGa                                      <  1:1632325/148‑1 (MQ=255)
                                                                                                        cTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGa                                      <  1:1632344/148‑1 (MQ=255)
                                                                                                          gAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGAtt                                    >  1:865380/1‑148 (MQ=255)
                                                                                                            aaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtga                                  >  2:3096242/1‑148 (MQ=255)
                                                                                                            aaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtga                                  >  2:3096266/1‑148 (MQ=255)
                                                                                                              gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatg                                <  2:865380/148‑1 (MQ=255)
                                                                                                              gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatg                                <  2:925090/148‑1 (MQ=255)
                                                                                                               aCATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatga                               <  1:1381963/148‑1 (MQ=255)
                                                                                                               aCATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatga                               <  2:2342111/148‑1 (MQ=255)
                                                                                                                  ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGacac                            <  1:2912562/148‑1 (MQ=255)
                                                                                                                    cAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacca                          <  1:1541502/148‑1 (MQ=255)
                                                                                                                    cAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacca                          >  2:197884/1‑148 (MQ=255)
                                                                                                                       gTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                       >  2:2151032/1‑148 (MQ=255)
                                                                                                                       gTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                       >  1:1389059/1‑148 (MQ=255)
                                                                                                                       gTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                       >  1:546873/1‑148 (MQ=255)
                                                                                                                       gTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                       >  1:1389232/1‑148 (MQ=255)
                                                                                                                        tCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtg                                   >  2:2353778/1‑135 (MQ=255)
                                                                                                                        tCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGt                      >  1:3298137/1‑148 (MQ=255)
                                                                                                                         cAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACACTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTa                     >  2:920469/1‑148 (MQ=255)
                                                                                                                          aaCATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTAt                    <  1:2403503/148‑1 (MQ=255)
                                                                                                                           aCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATc                   >  1:2410863/1‑148 (MQ=255)
                                                                                                                           aCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATc                   >  1:751049/1‑148 (MQ=255)
                                                                                                                           aCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGAACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATc                   >  1:375996/1‑148 (MQ=255)
                                                                                                                               cGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg               >  2:1479441/1‑148 (MQ=255)
                                                                                                                                  cAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATg            >  2:2705861/1‑148 (MQ=255)
                                                                                                                                    gTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGa          <  1:2642305/148‑1 (MQ=255)
                                                                                                                                     tAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGAc         >  1:1354505/1‑148 (MQ=255)
                                                                                                                                        aaaaaTGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTc      >  2:1658589/1‑148 (MQ=255)
                                                                                                                                         aaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCa     >  1:2586184/1‑148 (MQ=255)
                                                                                                                                         aaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCa     >  2:1273647/1‑148 (MQ=255)
                                                                                                                                         aaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCa     >  1:930886/1‑148 (MQ=255)
                                                                                                                                          aaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCAc    >  1:9343/1‑148 (MQ=255)
                                                                                                                                            aTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCt  >  2:2611426/1‑148 (MQ=255)
                                                                                                                                            aTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCt  >  2:1821277/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCT  >  CP009273/4434116‑4434403

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: