Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 285,869 G→T 17.2% L50I (CTC→ATC)  argF ← ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273285,8690GT17.2% 855.9 / 110.7 280L50I (CTC→ATC) argFornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (113/118);  new base T (13/36);  total (126/154)
Fisher's exact test for biased strand distribution p-value = 4.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

CCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCGGACAT  >  CP009273/285720‑286016
                                                                                                                                                     |                                                                                                                                                   
    ccGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGcgcg                                                                                                                                                 <  2:4387013/150‑1 (MQ=255)
      gagaACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTACGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGaa                                                                                                                                               >  1:1168282/1‑149 (MQ=255)
       agaACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATg                                                                                                                                              <  2:2486604/150‑1 (MQ=255)
                 cgGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCgg                                                                                                                                    <  1:562206/150‑1 (MQ=255)
                   gtgtCCTTAATTGACTCTTTATGCCCAATCGGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTa                                                                                                                                  <  2:1819943/150‑1 (MQ=255)
                    tgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTaa                                                                                                                                 >  2:1692454/1‑150 (MQ=255)
                    tgtCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTaa                                                                                                                                 <  2:1980301/150‑1 (MQ=255)
                      tCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGc                                                                                                                               <  2:1081456/150‑1 (MQ=255)
                      tCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGc                                                                                                                               >  1:2872711/1‑150 (MQ=255)
                      tCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGc                                                                                                                               >  2:2745594/1‑150 (MQ=255)
                       ccTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCt                                                                                                                              >  2:4392268/1‑150 (MQ=255)
                            aTTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGt                                                                                                                         <  1:3637274/150‑1 (MQ=255)
                              tGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTAc                                                                                                                       >  2:1741898/1‑150 (MQ=255)
                               gACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACt                                                                                                                      >  1:224309/1‑150 (MQ=255)
                                 ctctTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTc                                                                                                                    >  2:1553119/1‑150 (MQ=255)
                                       aTGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGcc                                                                                                              >  1:2285445/1‑150 (MQ=255)
                                            cAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCAtttt                                                                                                         <  1:268443/150‑1 (MQ=255)
                                              aTCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCAtttttt                                                                                                       <  2:2055441/150‑1 (MQ=255)
                                                 tGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTa                                                                                                    >  1:1332917/1‑150 (MQ=255)
                                                  ggCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTAt                                                                                                   <  1:1624438/150‑1 (MQ=255)
                                                    cTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCg                                                                                                 <  2:23196/150‑1 (MQ=255)
                                                     tGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAACGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCgg                                                                                                <  1:4072468/150‑1 (MQ=255)
                                                     tGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCgg                                                                                                <  1:1718952/150‑1 (MQ=255)
                                                          cTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTg                                                                                           >  1:251600/1‑150 (MQ=255)
                                                          cTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTg                                                                                           >  1:3179855/1‑150 (MQ=255)
                                                           tCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGa                                                                                          <  1:2413106/150‑1 (MQ=255)
                                                           tCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGa                                                                                          <  1:2509161/150‑1 (MQ=255)
                                                            cGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAg                                                                                         <  2:2446309/150‑1 (MQ=255)
                                                             ggCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGc                                                                                        >  1:3956335/1‑150 (MQ=255)
                                                                 cTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCtgtg                                                                                    >  1:1414858/1‑150 (MQ=255)
                                                                   aaaTAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCg                                                                                  <  1:2152705/150‑1 (MQ=255)
                                                                        ggTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAg                                                                             <  1:2657629/150‑1 (MQ=255)
                                                                           aaCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGt                                                                          >  2:673355/1‑150 (MQ=255)
                                                                           aaCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGt                                                                          >  1:2946541/1‑150 (MQ=255)
                                                                            aCGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTc                                                                         <  1:912060/150‑1 (MQ=255)
                                                                                gcgcgcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTcagca                                                                     <  1:1604312/150‑1 (MQ=255)
                                                                                gcgcgcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTcagca                                                                     <  1:1296902/150‑1 (MQ=255)
                                                                                    gcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAgaga                                                                 <  2:4575190/150‑1 (MQ=255)
                                                                                    gcCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAgaga                                                                 >  1:659317/1‑150 (MQ=255)
                                                                                       cTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGt                                                              <  1:2967426/150‑1 (MQ=255)
                                                                                       cTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGt                                                              <  1:651301/150‑1 (MQ=255)
                                                                                         ggTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGa                                                            >  1:800158/1‑150 (MQ=255)
                                                                                            cAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCGTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACt                                                         <  1:3823269/150‑1 (MQ=255)
                                                                                              aaTGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACtgt                                                       >  1:2819654/1‑150 (MQ=255)
                                                                                              aaTGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACtgt                                                       >  1:3756828/1‑150 (MQ=255)
                                                                                              aaTGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACtgt                                                       >  2:4148529/1‑150 (MQ=255)
                                                                                               aTGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTTAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACtgtg                                                      <  2:3855289/150‑1 (MQ=255)
                                                                                                 gCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCa                                                    <  1:2089670/150‑1 (MQ=255)
                                                                                                   ggCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAgg                                                  >  2:2110842/1‑150 (MQ=255)
                                                                                                     cAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAgggg                                                <  1:1024382/150‑1 (MQ=255)
                                                                                                      aaCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGt                                               <  2:950363/150‑1 (MQ=255)
                                                                                                       aCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCTCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTTCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTa                                              <  2:2484271/150‑1 (MQ=255)
                                                                                                       aCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTa                                              <  2:2582048/150‑1 (MQ=255)
                                                                                                       aCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTa                                              <  2:1383808/150‑1 (MQ=255)
                                                                                                       aCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTa                                              >  1:3286592/1‑150 (MQ=255)
                                                                                                       aCTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTa                                              <  1:3127483/150‑1 (MQ=255)
                                                                                                        cTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACGTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTaa                                             <  1:2168157/150‑1 (MQ=255)
                                                                                                        cTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTaa                                             <  2:3630648/150‑1 (MQ=255)
                                                                                                            gAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTc                                         <  2:1922379/150‑1 (MQ=255)
                                                                                                             aaaGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCg                                        >  2:2869889/1‑150 (MQ=255)
                                                                                                             aaaGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCg                                        <  2:1025291/150‑1 (MQ=255)
                                                                                                             aaaGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCg                                        <  1:3077768/150‑1 (MQ=255)
                                                                                                               agagCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAg                                      <  2:4068154/150‑1 (MQ=255)
                                                                                                                    aaCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGtt                                 >  2:506938/1‑150 (MQ=255)
                                                                                                                     aCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGttt                                <  2:2813372/150‑1 (MQ=255)
                                                                                                                       gggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCa                              <  1:442547/150‑1 (MQ=255)
                                                                                                                       gggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCa                              >  2:4201266/1‑150 (MQ=255)
                                                                                                                       gggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCa                              >  2:739699/1‑150 (MQ=255)
                                                                                                                        ggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAg                             <  1:4514471/150‑1 (MQ=255)
                                                                                                                           aCGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGaaa                          <  1:1811339/150‑1 (MQ=255)
                                                                                                                           aCGAGTCGAGTCTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACGTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGaaa                          <  1:4602060/150‑1 (MQ=255)
                                                                                                                            cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAg                         >  2:1188772/1‑150 (MQ=255)
                                                                                                                             gagtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAgt                        >  2:336223/1‑150 (MQ=255)
                                                                                                                               gtcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAgtgt                      >  1:3512651/1‑150 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCANGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  1:4276754/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  1:926359/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  2:886577/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  2:688927/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  2:3121639/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  2:1270757/150‑1 (MQ=255)
                                                                                                                                tcgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTATTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGtt                     <  2:810094/150‑1 (MQ=255)
                                                                                                                                 cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGttt                    >  2:2554198/1‑150 (MQ=255)
                                                                                                                                   tgtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGttttt                  <  2:3469827/149‑1 (MQ=255)
                                                                                                                                   agtcTTTTTCGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGttttt                  >  2:3787611/1‑150 (MQ=255)
                                                                                                                                   agtcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGttttt                  <  2:1281201/150‑1 (MQ=255)
                                                                                                                                     tcTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTATGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTg                <  1:4184898/150‑1 (MQ=255)
                                                                                                                                        ttttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCTGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGtat             <  1:76723/150‑1 (MQ=255)
                                                                                                                                         tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGtata            >  1:2886278/1‑150 (MQ=255)
                                                                                                                                         tttCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGtata            >  2:4374402/1‑150 (MQ=255)
                                                                                                                                            cGAAGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAAt         <  1:2715275/150‑1 (MQ=255)
                                                                                                                                            cGAAGATGAGCGCGATGTTTTTCCCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAAt         <  2:4011457/150‑1 (MQ=255)
                                                                                                                                            cGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAAt         >  2:3570740/1‑150 (MQ=255)
                                                                                                                                            cGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAAt         <  1:3406145/150‑1 (MQ=255)
                                                                                                                                              aaGATGATCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCg       <  2:3424510/150‑1 (MQ=255)
                                                                                                                                              aaGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCg       >  2:2293955/1‑150 (MQ=255)
                                                                                                                                              aaGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCg       >  1:3344370/1‑150 (MQ=255)
                                                                                                                                               aGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCgg      <  2:2015818/150‑1 (MQ=255)
                                                                                                                                               aGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCgg      >  2:2601156/1‑150 (MQ=255)
                                                                                                                                                gATGAGCGCGTTGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCGGa     <  1:705085/150‑1 (MQ=255)
                                                                                                                                                 aTGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCGGAc    <  1:3336988/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCCCGAGAACCCGCGCGGTGTCCTTAATTGACTCTTTATGCCCAATCTGGCTGCCGCTCGGCCCTAAATAGGTAACGCGCGCGCCCTGGTCAAATGCGGCAACTTCGAAAGAGCAACGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCGGTAAGCTTCTGTACTTCCTTGCCATTTTTTTTATCGGCTTTGAGCTGTGCGGCAAGGGTCAGCAGAGAAGTGAACTGTGCAGGGGTAAAGTCGAGCAGTTTCAGAAAGTGTTTTTTGTATAAATCGGACAT  >  CP009273/285720‑286016

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: