Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,488,443 T→A 11.3% P89P (CCT→CCA ypdI → putative lipoprotein involved in colanic acid biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,488,4430TA11.3% 758.8 / 62.0 263P89P (CCT→CCAypdIputative lipoprotein involved in colanic acid biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (104/127);  new base A (16/15);  total (121/142)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

AAAACATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACAT  >  CP009273/2488294‑2488591
                                                                                                                                                     |                                                                                                                                                    
aaaaCATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAAACTCACAGAAGAAATGACGCCCCCTTTTCCTACAAATGCGAATTATCTTATTCATACt                                                                                                                                                      >  1:3750106/1‑150 (MQ=255)
 aaaCATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCAGTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCtg                                                                                                                                                     >  2:2713263/1‑150 (MQ=255)
  aaCATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCtgt                                                                                                                                                    >  1:179284/1‑150 (MQ=255)
    cATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGc                                                                                                                                                  <  1:3325623/150‑1 (MQ=255)
    cATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGc                                                                                                                                                  <  2:1680976/150‑1 (MQ=255)
     aTCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCg                                                                                                                                                 <  1:4715017/150‑1 (MQ=255)
     aTCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCg                                                                                                                                                 >  1:7118/1‑150 (MQ=255)
         gCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAg                                                                                                                                             <  2:5262563/150‑1 (MQ=255)
         gCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAg                                                                                                                                             <  2:4723820/150‑1 (MQ=255)
          cGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGc                                                                                                                                            <  1:4382766/150‑1 (MQ=255)
              gacgacCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCt                                                                                                                                        <  1:90363/150‑1 (MQ=255)
                 gacCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTc                                                                                                                                     <  1:661490/150‑1 (MQ=255)
                   ccACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTcac                                                                                                                                   <  1:370008/150‑1 (MQ=255)
                     acacacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACAt                                                                                                                                 >  2:4619208/1‑150 (MQ=255)
                     acacacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACAt                                                                                                                                 <  1:4148369/150‑1 (MQ=255)
                      cacacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATg                                                                                                                                <  1:154611/150‑1 (MQ=255)
                      cacacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATg                                                                                                                                >  1:1648999/1‑150 (MQ=255)
                          ctctACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGCTAGCGTCTTTCACATGGCAg                                                                                                                            >  2:1659698/1‑150 (MQ=255)
                           tctACCAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAgg                                                                                                                           <  2:1245383/150‑1 (MQ=255)
                            ctACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGt                                                                                                                          >  2:4872321/1‑150 (MQ=255)
                               cAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGaa                                                                                                                       >  2:383533/1‑150 (MQ=255)
                                aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAAc                                                                                                                      >  2:2204084/1‑150 (MQ=255)
                                 aaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACa                                                                                                                     <  2:110555/150‑1 (MQ=255)
                                 aaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACa                                                                                                                     <  1:1145120/150‑1 (MQ=255)
                                    cTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATc                                                                                                                  <  1:3415926/150‑1 (MQ=255)
                                     ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCa                                                                                                                 <  2:2011696/150‑1 (MQ=255)
                                        aGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGt                                                                                                              >  2:977286/1‑150 (MQ=255)
                                         gCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                             <  1:1353256/150‑1 (MQ=255)
                                          ccATCACAGGAGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTct                                                                                                            <  2:577700/150‑1 (MQ=255)
                                           cATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctg                                                                                                           <  2:2257680/150‑1 (MQ=255)
                                            aTCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctga                                                                                                          <  1:86649/150‑1 (MQ=255)
                                             tCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgat                                                                                                         <  1:2057376/150‑1 (MQ=255)
                                              cacaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatg                                                                                                        <  2:4647560/150‑1 (MQ=255)
                                              cacaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatg                                                                                                        <  1:5043895/150‑1 (MQ=255)
                                               acaGGCGATTTTGTTTAAGATGAGATAAGATGCGCCACCATTAAAACTCACAGAAGAAATGCCGCCCCCATTTACAACAAATGCGAATTATATTATTCATCCATTGCGATAGCGTATTTAAAAAGGCAGGTGAACAATAAATTctgatgc                                                                                                       >  2:4668032/1‑150 (MQ=255)
                                                 aGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                                     >  1:1265639/1‑150 (MQ=255)
                                                 aGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                                     >  2:1458212/1‑150 (MQ=255)
                                                 aGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                                     >  1:240993/1‑150 (MQ=255)
                                                   gCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctgat                                                                                                   <  2:1909148/150‑1 (MQ=255)
                                                      aTTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGAtgtg                                                                                                <  2:1392921/150‑1 (MQ=255)
                                                       ttttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGa                                                                                               >  2:3853817/1‑150 (MQ=255)
                                                                aGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGcc                                                                                      >  2:153559/1‑150 (MQ=255)
                                                                 gATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGccc                                                                                     <  1:3869666/150‑1 (MQ=255)
                                                                  aTGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCa                                                                                    >  1:4839260/1‑150 (MQ=255)
                                                                       agaagaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGa                                                                               <  1:256938/150‑1 (MQ=255)
                                                                       agaagaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGa                                                                               >  2:4536362/1‑150 (MQ=255)
                                                                         aagaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGc                                                                             >  1:3462652/1‑150 (MQ=255)
                                                                         aagaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGc                                                                             >  1:4249715/1‑150 (MQ=255)
                                                                          agaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGcc                                                                            <  1:2995310/150‑1 (MQ=255)
                                                                           gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCa                                                                           <  2:3802518/150‑1 (MQ=255)
                                                                            aTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCAt                                                                          <  2:4202907/150‑1 (MQ=255)
                                                                             tGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATc                                                                         <  2:3693330/150‑1 (MQ=255)
                                                                               cgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATTGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCt                                                                       <  1:2902997/150‑1 (MQ=255)
                                                                                 ccaccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGa                                                                     >  1:3602394/1‑150 (MQ=255)
                                                                                  caccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAc                                                                    >  2:4464533/1‑150 (MQ=255)
                                                                                  caccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAc                                                                    <  2:1577291/150‑1 (MQ=255)
                                                                                  caccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAc                                                                    <  1:775512/150‑1 (MQ=255)
                                                                                    ccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAccc                                                                  >  2:708353/1‑150 (MQ=255)
                                                                                     caTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCa                                                                 >  1:4179896/1‑150 (MQ=255)
                                                                                       ttAAACCTCACAGAAGAAATGCCGCCCCCATTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACt                                                               <  1:1391987/150‑1 (MQ=255)
                                                                                         aaaCCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTAc                                                             >  1:614904/1‑150 (MQ=255)
                                                                                         aaaCCTCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTAc                                                             >  2:2856882/1‑150 (MQ=255)
                                                                                              tCACAGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCt                                                        <  1:2808606/150‑1 (MQ=255)
                                                                                                 caGAAGAAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGcc                                                     <  1:5103086/150‑1 (MQ=255)
                                                                                                  agaagaAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGTTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCa                                                    >  1:5003552/1‑150 (MQ=255)
                                                                                                  agaagaAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCa                                                    >  2:3032720/1‑150 (MQ=255)
                                                                                                  agaagaAATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCa                                                    <  2:240993/150‑1 (MQ=255)
                                                                                                     agaaATGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATAc                                                 <  2:5120909/150‑1 (MQ=255)
                                                                                                       aaaTGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATAccc                                               <  2:179284/150‑1 (MQ=255)
                                                                                                        aaTGCCGCCCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATAcccc                                              >  2:4942275/1‑150 (MQ=255)
                                                                                                           gccgccCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCccacc                                           >  2:1193870/1‑150 (MQ=255)
                                                                                                              gccCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATa                                        <  1:2391294/150‑1 (MQ=255)
                                                                                                              gccCCCTTTTCCTACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATa                                        >  1:896576/1‑150 (MQ=255)
                                                                                                               cccccTTTTCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATaa                                       >  2:3462664/1‑150 (MQ=255)
                                                                                                                      ttCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAat                                <  1:1014728/150‑1 (MQ=255)
                                                                                                                       tCCTACAAAGGCGAATTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAatc                               >  1:527563/1‑150 (MQ=255)
                                                                                                                        cctaCAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAatca                              >  1:4370055/4‑150 (MQ=255)
                                                                                                                        cctaCAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAatca                              >  2:2946676/4‑150 (MQ=255)
                                                                                                                        cctaCAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAatca                              >  2:3610716/4‑150 (MQ=255)
                                                                                                                          taCAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                            <  2:4166562/149‑1 (MQ=255)
                                                                                                                             aaaGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAg                         >  2:90222/1‑150 (MQ=255)
                                                                                                                             aaaGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAg                         <  1:1218534/150‑1 (MQ=255)
                                                                                                                                ggCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                      <  2:740966/150‑1 (MQ=255)
                                                                                                                                  cGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCcaca                    >  2:4337748/1‑150 (MQ=255)
                                                                                                                                   gAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACaa                   <  1:2256103/150‑1 (MQ=255)
                                                                                                                                   gAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACaa                   <  1:1907618/150‑1 (MQ=255)
                                                                                                                                    aaTTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAg                  <  1:4779457/150‑1 (MQ=255)
                                                                                                                                    aaTTATCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAg                  <  1:2674347/150‑1 (MQ=255)
                                                                                                                                      ttATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGcc                <  2:245580/150‑1 (MQ=255)
                                                                                                                                       tATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCg               >  2:5096028/1‑150 (MQ=255)
                                                                                                                                         tCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCggg             >  2:2585325/1‑150 (MQ=255)
                                                                                                                                         tCTTATTCATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCggg             <  2:4101232/150‑1 (MQ=255)
                                                                                                                                          cTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgggg            >  2:3114855/1‑150 (MQ=255)
                                                                                                                                                cATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa      <  1:1458157/150‑1 (MQ=255)
                                                                                                                                                cATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa      <  1:1100405/150‑1 (MQ=255)
                                                                                                                                                cATCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa      >  1:2616898/1‑150 (MQ=255)
                                                                                                                                                 aTCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaaa     >  2:1925484/1‑150 (MQ=255)
                                                                                                                                                 aTCCAGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaaa     <  1:1192600/150‑1 (MQ=255)
                                                                                                                                                  tCCAGTGCGATAGCGTCTTTCACATGGCAGGTCAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCAACCTGACCCACTCCCGCCTGCCATCCCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCCAGAAAc    >  1:4568042/1‑150 (MQ=255)
                                                                                                                                                    cTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACAt  <  2:323535/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACAT  >  CP009273/2488294‑2488591

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: