Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,398,168 C→A 9.6% pseudogene (78/300 nt) yrdE → pseudogene, yedZ paralog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,398,1680CA 9.6% 1127.8 / 56.3 302pseudogene (78/300 nt)yrdEpseudogene, yedZ paralog
Reads supporting (aligned to +/- strand):  ref base C (147/126);  new base A (16/13);  total (163/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 304 total aligned reads displayed.

CGCCTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTC  >  CP009273/3398021‑3398317
                                                                                                                                                   |                                                                                                                                                     
cGCCTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCaa                                                                                                                                                     <  1:1593332/150‑1 (MQ=255)
 gCCTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAAt                                                                                                                                                    >  2:3097266/1‑150 (MQ=255)
 gCCTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAAt                                                                                                                                                    >  1:1792298/1‑150 (MQ=255)
  ccTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATc                                                                                                                                                   >  1:4051143/1‑150 (MQ=255)
     ccGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGt                                                                                                                                                <  2:1760781/150‑1 (MQ=255)
     ccGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGt                                                                                                                                                <  1:2555309/150‑1 (MQ=255)
     ccGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGt                                                                                                                                                <  2:2341676/150‑1 (MQ=255)
            cTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAAtttt                                                                                                                                         <  1:367364/150‑1 (MQ=255)
            cTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAAtttt                                                                                                                                         >  2:4233284/1‑150 (MQ=255)
             tGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTg                                                                                                                                        <  1:3819131/150‑1 (MQ=255)
                 aGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTg                                                                                                                                    >  2:194538/1‑150 (MQ=255)
                  ggTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGCGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGc                                                                                                                                   >  1:244542/1‑150 (MQ=255)
                   gTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCt                                                                                                                                  >  1:3925996/1‑150 (MQ=255)
                         gTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCa                                                                                                                            >  1:1853729/1‑150 (MQ=255)
                            gccttGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTc                                                                                                                         <  2:4613103/147‑1 (MQ=255)
                             cTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTca                                                                                                                        >  1:1305090/1‑150 (MQ=255)
                                 gCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGGACCATGGGGTGAATAATCTGATTTTGTTTGACTACACATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCTCATc                                                                                                                    <  2:2087365/150‑1 (MQ=255)
                                  ccTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCt                                                                                                                   <  1:4517239/150‑1 (MQ=255)
                                  ccTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCt                                                                                                                   <  2:340891/150‑1 (MQ=255)
                                  ccTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCt                                                                                                                   >  1:4139890/1‑150 (MQ=255)
                                    tGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTAt                                                                                                                 >  1:2914850/1‑150 (MQ=255)
                                     gggCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATc                                                                                                                >  1:1057002/1‑150 (MQ=255)
                                      ggCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATcc                                                                                                               <  2:1263725/150‑1 (MQ=255)
                                         gACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGa                                                                                                            <  2:853672/150‑1 (MQ=255)
                                         gACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGGCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGa                                                                                                            >  2:2024431/1‑150 (MQ=255)
                                         gACGCTACGCTTAGCCCCTTACATATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCAGTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGa                                                                                                            >  2:2058774/1‑150 (MQ=255)
                                             ctacgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGAt                                                                                                        >  1:1705484/1‑150 (MQ=255)
                                             ctacgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGAt                                                                                                        <  1:1436411/150‑1 (MQ=255)
                                             ctacgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTGTGCGTGCGTTAGCCGTCACAGCGATCCAGACGag                                                                                                        >  1:2861384/1‑149 (MQ=255)
                                              tacgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATg                                                                                                       <  1:79187/150‑1 (MQ=255)
                                                cgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCa                                                                                                     >  2:5676/1‑150 (MQ=255)
                                                cgctTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCa                                                                                                     <  2:728709/150‑1 (MQ=255)
                                                   ttAGCCCCTTACTTATTTCTGGTACCATTGTGTGAATAATATGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTg                                                                                                  <  2:3711014/150‑1 (MQ=255)
                                                   ttAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTg                                                                                                  >  1:2143951/1‑150 (MQ=255)
                                                      gCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGaaa                                                                                               >  1:834585/1‑150 (MQ=255)
                                                        cccTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAAt                                                                                             >  1:285320/1‑150 (MQ=255)
                                                        cccTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAAt                                                                                             <  2:3346983/150‑1 (MQ=255)
                                                         ccTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAAtt                                                                                            <  2:2253473/150‑1 (MQ=255)
                                                         ccTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAAtt                                                                                            >  2:2296600/1‑150 (MQ=255)
                                                          cttacttaTTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTg                                                                                           <  2:3586087/150‑1 (MQ=255)
                                                           ttacttaTTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTgg                                                                                          >  1:1271008/1‑150 (MQ=255)
                                                            tacttaTTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTggg                                                                                         >  2:143318/1‑150 (MQ=255)
                                                                 aTTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATc                                                                                    >  2:4446639/1‑150 (MQ=255)
                                                                  tttCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATcc                                                                                   >  1:1465089/1‑150 (MQ=255)
                                                                       ggTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGc                                                                              >  2:2959719/1‑150 (MQ=255)
                                                                       ggTACCATGGGGTGAATAATATGATATTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGAGTGATTTAGCATTAACATCAATCCAGACGATGCAGTGAAAAATGGGTAAACCCCAGc                                                                              >  2:2735166/1‑150 (MQ=255)
                                                                           ccATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAcc                                                                          <  2:4200926/150‑1 (MQ=255)
                                                                           ccATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAAcc                                                                          <  2:2703538/150‑1 (MQ=255)
                                                                             agggggTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGc                                                                        <  2:187724/148‑1 (MQ=255)
                                                                              tGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCt                                                                       <  2:3274332/150‑1 (MQ=255)
                                                                               ggggTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTg                                                                      >  2:1058219/1‑150 (MQ=255)
                                                                                  gTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                   >  2:3262352/1‑150 (MQ=255)
                                                                                  gTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGt                                                                   >  2:2706564/1‑150 (MQ=255)
                                                                                     aataatCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAAt                                                                <  1:1003749/150‑1 (MQ=255)
                                                                                        aatCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtc                                                             <  2:346470/150‑1 (MQ=255)
                                                                                         atCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtcg                                                            <  1:855246/150‑1 (MQ=255)
                                                                                           cTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATgtcgtc                                                          >  2:514185/1‑150 (MQ=255)
                                                                                             gATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTa                                                        <  1:3705857/150‑1 (MQ=255)
                                                                                             gATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTa                                                        >  1:3309139/1‑150 (MQ=255)
                                                                                                 ttgtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                    <  1:4550805/150‑1 (MQ=255)
                                                                                                   gtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGt                                                  <  1:3608910/150‑1 (MQ=255)
                                                                                                     ttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcg                                                <  2:4047285/150‑1 (MQ=255)
                                                                                                      tgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgc                                               >  2:133760/1‑150 (MQ=255)
                                                                                                         cTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGAt                                            <  2:3034939/150‑1 (MQ=255)
                                                                                                            cAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCt                                         <  2:1698850/150‑1 (MQ=255)
                                                                                                            cAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCt                                         <  2:1397953/150‑1 (MQ=255)
                                                                                                             aaaTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                        <  2:1919638/150‑1 (MQ=255)
                                                                                                             aaaTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                        >  1:612970/1‑150 (MQ=255)
                                                                                                              aaTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTgg                                       >  2:738337/1‑150 (MQ=255)
                                                                                                              aaTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTgg                                       >  2:4122789/1‑150 (MQ=255)
                                                                                                                  aaTCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATc                                   <  1:3022760/150‑1 (MQ=255)
                                                                                                                     cACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCt                                >  1:226567/1‑150 (MQ=255)
                                                                                                                     cACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCt                                <  2:607000/150‑1 (MQ=255)
                                                                                                                           aaCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                          >  2:4121813/1‑150 (MQ=255)
                                                                                                                            aCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                         <  2:1233000/150‑1 (MQ=255)
                                                                                                                            aCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                         <  2:244379/150‑1 (MQ=255)
                                                                                                                            aCCTATTTAATGCTGAGCATTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                         <  2:4585850/150‑1 (MQ=255)
                                                                                                                              cTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                       <  2:632012/150‑1 (MQ=255)
                                                                                                                                aTTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtgt                     >  1:19704/1‑150 (MQ=255)
                                                                                                                                aTTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtgt                     <  2:899382/150‑1 (MQ=255)
                                                                                                                                aTTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtgt                     >  1:2271469/1‑150 (MQ=255)
                                                                                                                                   tAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                  <  1:1223420/150‑1 (MQ=255)
                                                                                                                                   tAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                  >  1:3162514/1‑150 (MQ=255)
                                                                                                                                   tAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGt                  >  1:1665040/1‑150 (MQ=255)
                                                                                                                                    aaTGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATACTGGCATCCCTACATTATTGGTGGTc                 >  2:254314/1‑150 (MQ=255)
                                                                                                                                      tGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtg               >  1:1428014/1‑150 (MQ=255)
                                                                                                                                       gCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgt              <  2:3925996/150‑1 (MQ=255)
                                                                                                                                        cTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg             >  2:3532047/1‑150 (MQ=255)
                                                                                                                                        cTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg             >  1:2392952/1‑150 (MQ=255)
                                                                                                                                        cTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg             <  2:2856854/150‑1 (MQ=255)
                                                                                                                                        cTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg             <  2:2167389/150‑1 (MQ=255)
                                                                                                                                        cTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtgtg             >  2:1525622/1‑150 (MQ=255)
                                                                                                                                         tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGa            >  1:932490/1‑150 (MQ=255)
                                                                                                                                         tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGa            >  2:1085124/1‑150 (MQ=255)
                                                                                                                                          gAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGaa           <  1:41848/150‑1 (MQ=255)
                                                                                                                                             cATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAt        >  1:502299/1‑150 (MQ=255)
                                                                                                                                             cATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAt        <  2:702548/150‑1 (MQ=255)
                                                                                                                                             cATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAt        >  1:2617406/1‑150 (MQ=255)
                                                                                                                                              aTTGTAAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGAtt       >  1:4276275/1‑150 (MQ=255)
                                                                                                                                                tGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTAt     >  2:3782321/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCCTCCGTCGCCTGTTAGGTTTATGTTGCTTTGCCTGGGCGACGCTACGCTTAGCCCCTTACTTATTTCTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTC  >  CP009273/3398021‑3398317

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: ‑