Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,248 A→T 100% 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,2480AT99.2% 829.7 / ‑2.7 250intergenic (+53/‑296)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (111/137);  total (112/138)
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.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 250 total aligned reads displayed.

ACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGA  >  CP009273/4122104‑4122391
                                                                                                                                                |                                                                                                                                               
 cccTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGa                                                                                                                                             <  1:1667295/148‑1 (MQ=255)
    tctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAat                                                                                                                                          >  2:1239048/1‑148 (MQ=255)
      tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatat                                                                                                                                        <  1:2819334/148‑1 (MQ=255)
      tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatat                                                                                                                                        <  1:1862243/148‑1 (MQ=255)
       ggTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatata                                                                                                                                       <  1:2642960/148‑1 (MQ=255)
       ggTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatata                                                                                                                                       >  1:161917/1‑148 (MQ=255)
        gTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatatat                                                                                                                                      >  1:912805/1‑148 (MQ=255)
               gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTa                                                                                                                               <  2:1151810/148‑1 (MQ=255)
                cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTAt                                                                                                                              <  2:191532/148‑1 (MQ=255)
                 ggACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATc                                                                                                                             <  1:1334130/148‑1 (MQ=255)
                   aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCAt                                                                                                                           <  1:3011064/148‑1 (MQ=255)
                   aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTATTATGAAATATATTGAATTATCAt                                                                                                                           >  1:2111438/1‑148 (MQ=255)
                      tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAgg                                                                                                                        >  1:496395/1‑148 (MQ=255)
                      tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAgg                                                                                                                        <  1:447608/148‑1 (MQ=255)
                      tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAgg                                                                                                                        >  1:3396117/1‑148 (MQ=255)
                       ggCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGa                                                                                                                       <  2:958354/148‑1 (MQ=255)
                       ggCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGa                                                                                                                       <  2:1073347/148‑1 (MQ=255)
                        gcgcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGAt                                                                                                                      >  1:2515466/1‑148 (MQ=255)
                         cgcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTGTCTTATGAAATATATTGAATTATCATAGGAtt                                                                                                                     <  2:1964819/148‑1 (MQ=255)
                          gcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTa                                                                                                                    <  2:1469127/148‑1 (MQ=255)
                             gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGc                                                                                                                 >  2:156071/1‑148 (MQ=255)
                                cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                                              <  1:2847962/148‑1 (MQ=255)
                                 gcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTTTAATTTCAGAAATTTAATAATGCCCGGTACTAATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGa                                                                                                             <  2:2951409/148‑1 (MQ=255)
                                        cACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCg                                                                                                      <  1:2014741/148‑1 (MQ=255)
                                         aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGt                                                                                                     >  2:2793002/1‑148 (MQ=255)
                                         aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGt                                                                                                     >  1:272930/1‑148 (MQ=255)
                                          ccgccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGtt                                                                                                    <  1:130060/148‑1 (MQ=255)
                                             ccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTac                                                                                                 >  2:2451552/1‑148 (MQ=255)
                                                gggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTTGGCCGGATTAAGCGTTTacgac                                                                                              >  2:2991738/1‑148 (MQ=255)
                                                  gCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTAcgacga                                                                                            <  1:3206284/148‑1 (MQ=255)
                                                    gATTCAGTCGGATATCAACCGGCTGTCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAAt                                                                                          <  1:3066979/148‑1 (MQ=255)
                                                    gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAAt                                                                                          <  1:315420/148‑1 (MQ=255)
                                                    gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAAt                                                                                          <  2:3073965/148‑1 (MQ=255)
                                                     aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                         <  1:1437630/148‑1 (MQ=255)
                                                     aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                         <  2:514357/148‑1 (MQ=255)
                                                     aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                         <  2:935418/148‑1 (MQ=255)
                                                     aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                         <  2:2978215/148‑1 (MQ=255)
                                                      ttCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                        >  2:1353806/1‑148 (MQ=255)
                                                       tCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCg                                                                                       >  2:1504945/1‑148 (MQ=255)
                                                           tCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaa                                                                                   <  2:2155022/148‑1 (MQ=255)
                                                            cGGATATCAACCGGCTGGCACAGTTGTTGTAATTTTAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaag                                                                                  <  2:3389779/148‑1 (MQ=255)
                                                             ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATTCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaaga                                                                                 >  1:1227419/1‑148 (MQ=255)
                                                             ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATGGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaaga                                                                                 <  2:2329393/148‑1 (MQ=255)
                                                             ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaaga                                                                                 <  2:2605268/148‑1 (MQ=255)
                                                               atatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                               <  2:2123268/148‑1 (MQ=255)
                                                               atatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                               <  2:2123116/148‑1 (MQ=255)
                                                                tatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGc                                                                              <  2:1204754/148‑1 (MQ=255)
                                                                   cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                           >  1:194553/1‑148 (MQ=255)
                                                                   cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                           >  2:1541680/1‑148 (MQ=255)
                                                                    aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATa                                                                          <  2:2330136/148‑1 (MQ=255)
                                                                    aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATa                                                                          <  2:1933594/148‑1 (MQ=255)
                                                                       cGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGt                                                                       >  2:680418/1‑148 (MQ=255)
                                                                                caGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTa                                                              >  2:1401343/1‑148 (MQ=255)
                                                                                       tgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtata                                                       <  1:205523/148‑1 (MQ=255)
                                                                                         tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATAtt                                                     >  1:1561466/1‑148 (MQ=255)
                                                                                           aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                   >  1:1876070/1‑148 (MQ=255)
                                                                                           aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                   >  1:1688505/1‑148 (MQ=255)
                                                                                            tttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGc                                                  >  1:758433/1‑148 (MQ=255)
                                                                                             ttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                 <  1:1030864/148‑1 (MQ=255)
                                                                                             ttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                 >  1:2238704/1‑148 (MQ=255)
                                                                                             ttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                 <  1:2783707/148‑1 (MQ=255)
                                                                                              tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGAGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAg                                                >  2:1862314/1‑148 (MQ=255)
                                                                                                  aaaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCg                                            >  2:3349757/1‑148 (MQ=255)
                                                                                                      ttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGttt                                        >  1:3451686/1‑148 (MQ=255)
                                                                                                        aataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTACGGtttgc                                      >  2:3241472/1‑148 (MQ=255)
                                                                                                         ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgct                                     <  1:1449832/148‑1 (MQ=255)
                                                                                                         ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgct                                     <  2:1225432/148‑1 (MQ=255)
                                                                                                         ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgct                                     <  1:2434848/148‑1 (MQ=255)
                                                                                                            atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttg                                  >  1:656391/1‑148 (MQ=255)
                                                                                                                 cGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCGTTGCATAc                             >  1:1837082/1‑148 (MQ=255)
                                                                                                                      cTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAt                        >  2:2057188/1‑148 (MQ=255)
                                                                                                                      cTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAt                        <  1:156528/148‑1 (MQ=255)
                                                                                                                       tCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtt                       <  2:89966/148‑1 (MQ=255)
                                                                                                                        cATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCTGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                      >  2:2507973/1‑148 (MQ=255)
                                                                                                                        cATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                      >  1:2659408/1‑148 (MQ=255)
                                                                                                                         aTGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                     <  1:789780/148‑1 (MQ=255)
                                                                                                                         aTGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                     >  1:404884/1‑148 (MQ=255)
                                                                                                                          tGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                    <  2:2839422/148‑1 (MQ=255)
                                                                                                                          tGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                    <  2:3162296/148‑1 (MQ=255)
                                                                                                                          tGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGAAGTCCGGTTTGCTTTGCATACCGGAttttc                    >  2:3208118/1‑148 (MQ=255)
                                                                                                                           gTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                   <  2:2390745/148‑1 (MQ=255)
                                                                                                                           gTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                   >  2:3123769/1‑148 (MQ=255)
                                                                                                                             tttCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                 <  2:2243996/148‑1 (MQ=255)
                                                                                                                             tttCGGGTTTATGGTTTCTTATGAAATATATTGAAGTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                 <  2:2243867/148‑1 (MQ=255)
                                                                                                                               tCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt               <  1:1239048/148‑1 (MQ=255)
                                                                                                                                cGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc              <  1:1504945/148‑1 (MQ=255)
                                                                                                                                  ggTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt            <  2:1163085/148‑1 (MQ=255)
                                                                                                                                   gTTTATGGTTTCTTATGACATATATTGAATTATCATAGGAATAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa           <  2:1809123/148‑1 (MQ=255)
                                                                                                                                   gTTTATGGTTTCTTATGAAATATCTTGACTTATCAGAGGCTGAGGCCGGATTAAGCGTTTACGACGACTCCGGCAAGAAGCAATAAGTACATGGTTAGATTCTATTTGCAGTCCGGTTTGCTTTGCAGACCGGATTTTCTTTTTCTTa           <  2:784438/148‑1 (MQ=255)
                                                                                                                                   gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa           <  2:215271/148‑1 (MQ=255)
                                                                                                                                    tttATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc          >  2:1741496/1‑148 (MQ=255)
                                                                                                                                    tttATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc          >  2:1741535/1‑148 (MQ=255)
                                                                                                                                     ttATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc         <  1:2451552/148‑1 (MQ=255)
                                                                                                                                        tGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc      >  2:3476962/1‑148 (MQ=255)
                                                                                                                                        tGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc      >  2:57905/1‑148 (MQ=255)
                                                                                                                                        tGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATc      >  2:1177567/1‑148 (MQ=255)
                                                                                                                                         ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtat                                                        <  1:1948641/97‑1 (MQ=255)
                                                                                                                                           tttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg   >  2:420103/1‑148 (MQ=255)
                                                                                                                                           tttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg   >  2:636899/1‑148 (MQ=255)
                                                                                                                                            ttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGa  >  1:2386604/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGA  >  CP009273/4122104‑4122391

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: