Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,248 A→T 44.5% intergenic (+53/‑296) metL → / → metF Bifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,122,2480AT44.5% 170.3 / inf 353intergenic (+53/‑296)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (87/109);  new base T (54/103);  total (141/212)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 366 total aligned reads displayed.

CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTT  >  CP009273/4122099‑4122397
                                                                                                                                                     |                                                                                                                                                     
 gATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTaa                                                                                                                                                      >  2:2824306/1‑150 (MQ=255)
   tAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATg                                                                                                                                                    <  2:2398976/150‑1 (MQ=255)
    aaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGa                                                                                                                                                   <  1:1288519/150‑1 (MQ=255)
             gTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTg                                                                                                                                          <  1:3487824/150‑1 (MQ=255)
              tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGa                                                                                                                                         >  2:3116354/1‑150 (MQ=255)
              tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTGTCGGGTTTATGGTTTCTTATGAAATATATTGa                                                                                                                                         >  2:3884107/1‑150 (MQ=255)
                aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAAt                                                                                                                                       <  1:584501/150‑1 (MQ=255)
                aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAAt                                                                                                                                       <  1:1816970/150‑1 (MQ=255)
                 tCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAAtt                                                                                                                                      <  1:475862/150‑1 (MQ=255)
                  ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATACTGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTa                                                                                                                                     <  2:2975763/150‑1 (MQ=255)
                   cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTAt                                                                                                                                    <  1:1742213/150‑1 (MQ=255)
                   cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTAt                                                                                                                                    <  2:37389/150‑1 (MQ=255)
                    gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATc                                                                                                                                   <  2:4995081/150‑1 (MQ=255)
                    gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATc                                                                                                                                   <  2:2654029/150‑1 (MQ=255)
                    gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATc                                                                                                                                   <  2:215303/150‑1 (MQ=255)
                     cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCa                                                                                                                                  <  1:884372/150‑1 (MQ=255)
                      ggACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCAt                                                                                                                                 <  2:3166588/150‑1 (MQ=255)
                        aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAg                                                                                                                               <  2:1980502/150‑1 (MQ=255)
                          cTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGa                                                                                                                             >  1:2478347/1‑150 (MQ=255)
                            ggCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGAtt                                                                                                                           <  1:1435553/150‑1 (MQ=255)
                                       cgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAACTATATTGAATTATCATAGGATTAGGCCGGATTa                                                                                                                >  2:4560124/1‑150 (MQ=255)
                                           gTCACCGCCGGGGCGAGTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCg                                                                                                            <  2:4342871/150‑1 (MQ=255)
                                              aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGttt                                                                                                         <  1:4999227/150‑1 (MQ=255)
                                              aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGttt                                                                                                         <  1:4141948/150‑1 (MQ=255)
                                              aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGttt                                                                                                         <  2:1168741/150‑1 (MQ=255)
                                                cgccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTac                                                                                                       >  1:2832824/1‑150 (MQ=255)
                                                 gccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacg                                                                                                      <  1:3173192/150‑1 (MQ=255)
                                                   cgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgac                                                                                                    >  1:2153301/1‑150 (MQ=255)
                                                      ggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGaa                                                                                                 <  2:1659259/150‑1 (MQ=255)
                                                        cGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                               >  1:3027370/1‑150 (MQ=255)
                                                         gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                              <  2:415656/150‑1 (MQ=255)
                                                          aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCg                                                                                             <  2:3685543/150‑1 (MQ=255)
                                                            tCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGc                                                                                           <  2:1272780/150‑1 (MQ=255)
                                                             cAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCa                                                                                          <  2:721435/150‑1 (MQ=255)
                                                              aGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaa                                                                                         <  2:4136692/150‑1 (MQ=255)
                                                              aGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaa                                                                                         >  1:4165920/1‑150 (MQ=255)
                                                               gTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCAGGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaag                                                                                        >  2:3913860/1‑150 (MQ=255)
                                                                  ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                                     <  2:1400921/150‑1 (MQ=255)
                                                                      atCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                                 <  2:4611052/150‑1 (MQ=255)
                                                                        cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATaa                                                                               >  1:3384775/1‑150 (MQ=255)
                                                                        cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATaa                                                                               >  2:3472730/1‑150 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                              >  2:324422/1‑150 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                              <  1:632761/150‑1 (MQ=255)
                                                                          aCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGt                                                                             <  1:401705/150‑1 (MQ=255)
                                                                              gCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACAt                                                                         <  1:4294457/150‑1 (MQ=255)
                                                                                tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                       >  1:5094397/1‑150 (MQ=255)
                                                                                tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                       <  1:1086385/150‑1 (MQ=255)
                                                                                   cacaGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTa                                                                    <  2:1842997/150‑1 (MQ=255)
                                                                                   cacaGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTa                                                                    >  2:1795265/1‑150 (MQ=255)
                                                                                    acaGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAg                                                                   <  2:3053151/150‑1 (MQ=255)
                                                                                      aGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATACGTACATGGGTAGtt                                                                 >  2:1542432/1‑150 (MQ=255)
                                                                                         tgttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtat                                                              <  1:1001092/150‑1 (MQ=255)
                                                                                         tgttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtat                                                              <  2:518563/150‑1 (MQ=255)
                                                                                           ttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTatat                                                            >  1:2337212/1‑150 (MQ=255)
                                                                                           ttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTGTatat                                                            >  1:2802091/1‑150 (MQ=255)
                                                                                             gtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATAttt                                                          >  1:2023557/1‑150 (MQ=255)
                                                                                             gtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATAttt                                                          <  1:1519215/150‑1 (MQ=255)
                                                                                             gtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATAttt                                                          <  2:3527319/150‑1 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                         >  2:3273361/1‑150 (MQ=255)
                                                                                               caTTTCAGAAATTTCATAATGCCCGGTACTCCTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGc                                                        <  1:4471865/149‑1 (MQ=255)
                                                                                                aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                       >  2:824770/1‑150 (MQ=255)
                                                                                                aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                       >  1:4733798/1‑150 (MQ=255)
                                                                                                   tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTc                                                    >  1:364712/1‑150 (MQ=255)
                                                                                                    cAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTcc                                                   <  2:749966/150‑1 (MQ=255)
                                                                                                     aGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCg                                                  <  2:2703798/150‑1 (MQ=255)
                                                                                                       aaaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGt                                                <  1:722621/150‑1 (MQ=255)
                                                                                                         aTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGttt                                              <  2:399457/150‑1 (MQ=255)
                                                                                                          tttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttg                                             >  1:1621618/1‑150 (MQ=255)
                                                                                                          tttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttg                                             >  1:380326/1‑150 (MQ=255)
                                                                                                               taatGCCCGGTACTCATGTTTTCGGTTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATACGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttg                                        >  1:969033/1‑150 (MQ=255)
                                                                                                                 atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                      <  2:4772375/150‑1 (MQ=255)
                                                                                                                 atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                      <  2:1911970/150‑1 (MQ=255)
                                                                                                                  tGCCCGTTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                     <  1:324240/150‑1 (MQ=255)
                                                                                                                  tGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCag                                     >  1:3596248/1‑149 (MQ=255)
                                                                                                                  tGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                     <  1:4844074/150‑1 (MQ=255)
                                                                                                                  tGCCCGGTACTCAGGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATAATCATAGGATTAGGACGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                     >  1:1069930/1‑150 (MQ=255)
                                                                                                                     ccGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAcc                                  <  2:590502/150‑1 (MQ=255)
                                                                                                                      cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                 >  2:2236856/1‑150 (MQ=255)
                                                                                                                          aCTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtt                             <  1:2993297/150‑1 (MQ=255)
                                                                                                                            tCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                           >  2:899086/1‑150 (MQ=255)
                                                                                                                            tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                           <  1:4098369/150‑1 (MQ=255)
                                                                                                                               tGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctt                        >  1:3055800/1‑150 (MQ=255)
                                                                                                                                 ttttCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttt                      <  2:3261216/150‑1 (MQ=255)
                                                                                                                                 ttttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttt                      <  1:3841374/150‑1 (MQ=255)
                                                                                                                                  tttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                     >  1:266254/1‑150 (MQ=255)
                                                                                                                                  tttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                     <  2:2503326/150‑1 (MQ=255)
                                                                                                                                  gttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                     <  1:2739018/149‑1 (MQ=255)
                                                                                                                                   ttCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                    >  2:4030926/1‑150 (MQ=255)
                                                                                                                                   gtCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                    <  2:4501844/149‑1 (MQ=255)
                                                                                                                                     cGGGGTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCACTAAGTACATGGTTAGTGTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttcttgttctt                  >  2:2860979/1‑150 (MQ=255)
                                                                                                                                      gggTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa                 >  1:2528900/1‑150 (MQ=255)
                                                                                                                                        gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc               <  2:2077090/150‑1 (MQ=255)
                                                                                                                                        gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc               <  2:4552443/150‑1 (MQ=255)
                                                                                                                                                tttCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGaa       <  2:4953988/150‑1 (MQ=255)
                                                                                                                                                  tCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt     <  1:1011018/150‑1 (MQ=255)
                                                                                                                                                  tCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt     <  2:4691465/150‑1 (MQ=255)
                                                                                                                                                  tCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt     <  2:3585410/150‑1 (MQ=255)
                                                                                                                                                  tCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt     <  1:1681151/150‑1 (MQ=255)
                                                                                                                                                    tAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGttt   >  1:4540365/1‑150 (MQ=255)
                                                                                                                                                     aaTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttt  <  2:70678/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTT  >  CP009273/4122099‑4122397

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: