Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 604,824 C→T 23.8% intergenic (+191/+91) insL1 → / ← entD IS186 transposase/phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273604,8240CT23.8% 508.6 / 166.1 244intergenic (+191/+91)insL1/entDIS186 transposase/phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base C (93/93);  new base T (22/36);  total (115/129)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 249 total aligned reads displayed.

AGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGCACAGC  >  CP009273/604675‑604971
                                                                                                                                                     |                                                                                                                                                   
aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATcc                                                                                                                                                     >  2:1988606/1‑150 (MQ=255)
aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATcc                                                                                                                                                     <  2:3220937/150‑1 (MQ=255)
aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCCTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATcc                                                                                                                                                     >  2:4067012/1‑150 (MQ=255)
    ggCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGtt                                                                                                                                                 >  1:672700/1‑150 (MQ=255)
       cATTCGTGGGACAGTATCACTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGAATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCTCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                                              <  2:2036001/150‑1 (MQ=255)
        aTTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtca                                                                                                                                             >  1:2377751/1‑150 (MQ=255)
             tGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtca                                                                                                                                        >  2:4783072/1‑150 (MQ=255)
             tGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtca                                                                                                                                        >  2:4782654/1‑150 (MQ=255)
             tGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtca                                                                                                                                        >  1:3653849/1‑150 (MQ=255)
              gggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTgtcaggtcag                                                                                                                                       >  1:1903233/1‑150 (MQ=255)
              gggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTgtcaggtcag                                                                                                                                       <  2:3398539/150‑1 (MQ=255)
              gggACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggtcag                                                                                                                                       >  2:1713628/1‑150 (MQ=255)
                gACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGa                                                                                                                                     <  1:10801/150‑1 (MQ=255)
                gACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGa                                                                                                                                     <  1:2495491/150‑1 (MQ=255)
                gACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGa                                                                                                                                     <  1:2438221/150‑1 (MQ=255)
                      aTCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                                                                                                               <  2:399684/150‑1 (MQ=255)
                       tCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAg                                                                                                                              <  2:3722232/150‑1 (MQ=255)
                        cccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAgg                                                                                                                             <  1:1238930/150‑1 (MQ=255)
                         ccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGc                                                                                                                            >  2:202841/1‑150 (MQ=255)
                            gACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCAcc                                                                                                                         >  2:773841/1‑150 (MQ=255)
                              cAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTg                                                                                                                       <  2:14863/150‑1 (MQ=255)
                                  cTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCttt                                                                                                                   <  2:1385411/150‑1 (MQ=255)
                                    aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTac                                                                                                                 <  2:2812628/150‑1 (MQ=255)
                                    aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTac                                                                                                                 <  2:178656/150‑1 (MQ=255)
                                     cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTaca                                                                                                                >  2:346756/1‑150 (MQ=255)
                                     cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTaca                                                                                                                >  2:1366186/1‑150 (MQ=255)
                                       aaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACt                                                                                                              >  2:2312122/1‑150 (MQ=255)
                                       aaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACt                                                                                                              >  2:1051603/1‑150 (MQ=255)
                                        aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACtt                                                                                                             <  2:3002485/150‑1 (MQ=255)
                                        aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACtt                                                                                                             >  2:1860637/1‑150 (MQ=255)
                                        aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACtt                                                                                                             <  2:3760168/150‑1 (MQ=255)
                                        aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACtt                                                                                                             >  2:823997/1‑150 (MQ=255)
                                         aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACttt                                                                                                            >  2:1285901/1‑150 (MQ=255)
                                          cGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTc                                                                                                           <  2:290199/150‑1 (MQ=255)
                                           gCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCg                                                                                                          >  1:4222714/1‑150 (MQ=255)
                                               ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTCAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCg                                                                                                      >  2:3209182/1‑150 (MQ=255)
                                               ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCg                                                                                                      >  2:4473141/1‑150 (MQ=255)
                                                tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGt                                                                                                     >  2:204335/1‑150 (MQ=255)
                                                tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGt                                                                                                     <  1:2121037/150‑1 (MQ=255)
                                                 gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                                    >  1:3381685/1‑150 (MQ=255)
                                                 gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                                    <  1:196624/150‑1 (MQ=255)
                                                  aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTgg                                                                                                   >  1:1163095/1‑150 (MQ=255)
                                                   agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGt                                                                                                  <  2:2155040/150‑1 (MQ=255)
                                                    gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTc                                                                                                 <  2:2711985/150‑1 (MQ=255)
                                                    gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTc                                                                                                 <  2:979116/150‑1 (MQ=255)
                                                         cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGa                                                                                            <  1:1064674/150‑1 (MQ=255)
                                                                   tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTAAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCggg                                                                                  >  2:3039623/1‑150 (MQ=255)
                                                                    cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGc                                                                                 <  2:907253/150‑1 (MQ=255)
                                                                     ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                                                <  2:1562799/150‑1 (MQ=255)
                                                                     ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                                                >  1:3522327/1‑150 (MQ=255)
                                                                     ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                                                <  1:184858/150‑1 (MQ=255)
                                                                      tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGa                                                                               <  2:3178102/150‑1 (MQ=255)
                                                                         cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATc                                                                            <  1:2095691/150‑1 (MQ=255)
                                                                            ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTa                                                                         >  1:2746446/1‑150 (MQ=255)
                                                                            ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTa                                                                         <  1:1466532/150‑1 (MQ=255)
                                                                                   aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtg                                                                  <  1:890424/150‑1 (MQ=255)
                                                                                   aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtg                                                                  >  1:50909/1‑150 (MQ=255)
                                                                                      cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgttg                                                               >  1:967530/1‑150 (MQ=255)
                                                                                       gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgttgt                                                              >  1:835179/1‑150 (MQ=255)
                                                                                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCaa                                                           >  2:1953814/1‑150 (MQ=255)
                                                                                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCaa                                                           >  1:3726152/1‑150 (MQ=255)
                                                                                          gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCaa                                                           <  2:504124/150‑1 (MQ=255)
                                                                                           cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAt                                                          <  1:4592232/150‑1 (MQ=255)
                                                                                            aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAtt                                                         <  1:1247146/150‑1 (MQ=255)
                                                                                            aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAtt                                                         <  1:3202130/150‑1 (MQ=255)
                                                                                            aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAtt                                                         >  1:1015610/1‑150 (MQ=255)
                                                                                            aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAAtt                                                         <  1:39134/150‑1 (MQ=255)
                                                                                             ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTa                                                        >  1:1216130/1‑150 (MQ=255)
                                                                                                 tcAGATGGTTCGACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCg                                                    >  2:1826388/1‑150 (MQ=255)
                                                                                                 tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCg                                                    <  1:3166480/150‑1 (MQ=255)
                                                                                                 tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCg                                                    <  1:3050128/150‑1 (MQ=255)
                                                                                                   aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCgtg                                                  >  1:4370810/1‑150 (MQ=255)
                                                                                                    gATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCgtgt                                                 <  2:3339998/150‑1 (MQ=255)
                                                                                                      tGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTg                                               >  2:433281/1‑150 (MQ=255)
                                                                                                        gTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGc                                             <  1:4143708/150‑1 (MQ=255)
                                                                                                        gTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGc                                             >  1:2958914/1‑150 (MQ=255)
                                                                                                         ttCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGca                                            >  1:3115274/1‑150 (MQ=255)
                                                                                                              cccTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCg                                       <  1:1745862/150‑1 (MQ=255)
                                                                                                                 ttAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTa                                    <  2:4653580/150‑1 (MQ=255)
                                                                                                                  tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTAt                                   >  1:543554/1‑150 (MQ=255)
                                                                                                                  tAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTAt                                   <  2:4562824/150‑1 (MQ=255)
                                                                                                                   aaGTTAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATg                                  >  2:2428043/1‑150 (MQ=255)
                                                                                                                       tAGCGCTTATGGGAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTa                              >  1:1604087/1‑150 (MQ=255)
                                                                                                                                tGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtcttttctt                     >  1:4182941/1‑150 (MQ=255)
                                                                                                                                 gggAGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtcttttcttt                    <  2:1917407/150‑1 (MQ=255)
                                                                                                                                  ggAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATcttttctttt                   <  1:202900/150‑1 (MQ=255)
                                                                                                                                   gAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGGGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTa                  >  2:745305/1‑150 (MQ=255)
                                                                                                                                    aGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAt                 <  2:84911/150‑1 (MQ=255)
                                                                                                                                    aGTAATCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAt                 <  1:2719533/150‑1 (MQ=255)
                                                                                                                                    aGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTAt                 >  1:1227045/1‑150 (MQ=255)
                                                                                                                                       aaTCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTg              >  2:1165099/1‑150 (MQ=255)
                                                                                                                                         tCCCCGCATATCTGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGcc            >  1:2893926/1‑150 (MQ=255)
                                                                                                                                         tCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGcc            >  2:2736475/1‑150 (MQ=255)
                                                                                                                                             cGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTg        <  2:50925/150‑1 (MQ=255)
                                                                                                                                              gCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGc       <  1:4640351/150‑1 (MQ=255)
                                                                                                                                               cATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGca      >  2:2595873/1‑150 (MQ=255)
                                                                                                                                                atatCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGcac     >  1:4444032/1‑150 (MQ=255)
                                                                                                                                                 tatCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGcaca    <  1:4782654/150‑1 (MQ=255)
                                                                                                                                                 tatCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGcaca    <  1:4783072/150‑1 (MQ=255)
                                                                                                                                                   tCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGCACAGc  >  1:1512996/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCTTTTATCTGCCAGTGCACAGC  >  CP009273/604675‑604971

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: