Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,896,477 C→A 6.6% P103H (CCC→CAC)  nudL → putative NUDIX hydrolase with low 3‑phosphohydroxypyruvate phosphatase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,896,4770CA 6.6% 510.0 / 14.2 151P103H (CCC→CAC) nudLputative NUDIX hydrolase with low 3‑phosphohydroxypyruvate phosphatase activity
Reads supporting (aligned to +/- strand):  ref base C (69/72);  new base A (7/3);  total (76/75)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 170 total aligned reads displayed.

TCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTAG  >  NC_000913/1896328‑1896626
                                                                                                                                                     |                                                                                                                                                     
tCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTgccgcc                                                                                                                                                       <  1:2325249/150‑1 (MQ=255)
tCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTgccgcc                                                                                                                                                       <  2:1575524/150‑1 (MQ=255)
 cATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGccc                                                                                                                                                      >  2:1713559/1‑150 (MQ=255)
 cATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAAGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGccc                                                                                                                                                      >  2:1930351/1‑150 (MQ=255)
     tGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCg                                                                                                                                                  <  1:1203619/150‑1 (MQ=255)
        gTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATa                                                                                                                                               >  1:2031951/1‑150 (MQ=255)
         tAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAg                                                                                                                                              >  2:528177/1‑150 (MQ=255)
          aaaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGc                                                                                                                                             <  2:1299404/150‑1 (MQ=255)
           aaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCg                                                                                                                                            <  1:193513/150‑1 (MQ=255)
             cacGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGACGGGCGGCAGCACGCCGATaacttc                                                                                                                                          <  2:515895/150‑6 (MQ=255)
              acGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGACGGGCGGCAGCACGCCGATAACTTCa                                                                                                                                         <  1:1002974/150‑1 (MQ=255)
                    gACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCt                                                                                                                                   >  1:609692/1‑150 (MQ=255)
                     aCAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTa                                                                                                                                  <  1:207894/150‑1 (MQ=255)
                      cAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTAc                                                                                                                                 >  2:858765/1‑150 (MQ=255)
                      cAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTAc                                                                                                                                 >  2:2199399/1‑150 (MQ=255)
                      cAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTAc                                                                                                                                 >  1:1294224/1‑150 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               >  1:1523653/1‑150 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               >  2:430895/1‑150 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               <  2:1412138/150‑1 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               <  1:1952088/150‑1 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               >  2:1062063/1‑150 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGGCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGGTACCa                                                                                                                               <  2:626210/150‑1 (MQ=255)
                        aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGGCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                               <  1:1894185/150‑1 (MQ=255)
                           ggCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGt                                                                                                                            >  2:1000223/1‑150 (MQ=255)
                              aTTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGCAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAAc                                                                                                                         >  2:141997/1‑150 (MQ=255)
                              aTTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAAc                                                                                                                         >  1:188634/1‑150 (MQ=255)
                                 cccTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAAcccc                                                                                                                      <  2:2316576/150‑1 (MQ=255)
                                 cccTGGAGGTGCAGTCGAAGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGCTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACAGGCTACCAGGTAAcccc                                                                                                                      >  2:1657238/1‑150 (MQ=255)
                                  ccTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCa                                                                                                                     <  2:2067139/150‑1 (MQ=255)
                                  ccTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCa                                                                                                                     <  2:1999389/150‑1 (MQ=255)
                                    tGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGt                                                                                                                   <  2:2487840/150‑1 (MQ=255)
                                     ggAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTg                                                                                                                  >  2:1565044/1‑150 (MQ=255)
                                        ggTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGAGGTc                                                                                                               >  1:783614/1‑150 (MQ=255)
                                          tGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCgg                                                                                                             >  2:312649/1‑150 (MQ=255)
                                            cAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCa                                                                                                           >  1:91281/1‑150 (MQ=255)
                                             aGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCAt                                                                                                          <  1:103363/150‑1 (MQ=255)
                                               tCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTa                                                                                                        <  1:480273/150‑1 (MQ=255)
                                                cGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTAt                                                                                                       <  1:1359020/150‑1 (MQ=255)
                                                cGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTAt                                                                                                       >  2:1177170/1‑150 (MQ=255)
                                                 gATGACACGGACGCATCAGCTATCGCCGCCGCGCTGTGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATc                                                                                                      >  1:1378284/1‑150 (MQ=255)
                                                  aTGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTAAcc                                                                                                     >  1:855269/1‑150 (MQ=255)
                                                   tGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCTGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATccc                                                                                                    >  2:1252958/1‑150 (MQ=255)
                                                   tGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTAGGCATTATccc                                                                                                    >  1:3951/1‑150 (MQ=255)
                                                        cGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATcccgcccg                                                                                               >  1:370021/1‑150 (MQ=255)
                                                            cGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCt                                                                                           >  2:2227715/1‑150 (MQ=255)
                                                             gCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTg                                                                                          <  1:813528/150‑1 (MQ=255)
                                                                 cAGCTATCGCCGCCGCGTTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGt                                                                                      <  1:1857188/150‑1 (MQ=255)
                                                                 cAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGt                                                                                      <  1:1713559/150‑1 (MQ=255)
                                                                    cTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATc                                                                                   >  2:145337/1‑150 (MQ=255)
                                                                          ccgccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCa                                                                             >  2:1725359/1‑150 (MQ=255)
                                                                           cgccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAg                                                                            <  2:2031951/150‑1 (MQ=255)
                                                                             ccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTg                                                                          >  1:1853521/1‑150 (MQ=255)
                                                                             ccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACAGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTg                                                                          >  1:1767682/1‑150 (MQ=255)
                                                                                   tGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATg                                                                    <  1:430895/150‑1 (MQ=255)
                                                                                    gcgcGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGa                                                                   >  2:20243/1‑150 (MQ=255)
                                                                                    gcgcGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGa                                                                   >  1:1481506/1‑150 (MQ=255)
                                                                                     cgcgAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGaa                                                                  <  2:1249489/150‑1 (MQ=255)
                                                                                      gcgAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAg                                                                 <  1:2479841/150‑1 (MQ=255)
                                                                                             tGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTcggc                                                          <  2:2390591/150‑1 (MQ=255)
                                                                                                    gagGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGttt                                                   <  1:822237/150‑1 (MQ=255)
                                                                                                     agGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGGCTCGGCGGTGTTTg                                                  <  1:1215058/150‑1 (MQ=255)
                                                                                                      ggTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGa                                                 >  2:1895809/1‑150 (MQ=255)
                                                                                                       gTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGaa                                                <  1:2488439/150‑1 (MQ=255)
                                                                                                       gTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGaa                                                <  1:1757966/150‑1 (MQ=255)
                                                                                                        tCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGCGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGaaa                                               >  2:719434/1‑150 (MQ=255)
                                                                                                         cGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAAt                                              <  1:1562073/150‑1 (MQ=255)
                                                                                                              taCCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGc                                         <  1:858765/150‑1 (MQ=255)
                                                                                                               aCCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGCCGGTGTTTGAAATGCCGCt                                        >  2:2379127/1‑150 (MQ=255)
                                                                                                                ccGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTc                                       >  1:1959554/1‑150 (MQ=255)
                                                                                                                ccGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTc                                       >  2:1614919/1‑150 (MQ=255)
                                                                                                                    cTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGccc                                   <  1:2227715/150‑1 (MQ=255)
                                                                                                                        cgccgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGc                               >  2:1856519/1‑150 (MQ=255)
                                                                                                                        cgccgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGc                               <  2:1496800/150‑1 (MQ=255)
                                                                                                                        cgccgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGc                               <  1:493737/150‑1 (MQ=255)
                                                                                                                            gTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTa                           >  1:2167379/1‑150 (MQ=255)
                                                                                                                                acgctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATc                       >  1:65904/5‑150 (MQ=255)
                                                                                                                                 cgctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCt                      >  2:352903/4‑150 (MQ=255)
                                                                                                                                  gctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAACATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTg                     >  2:472977/3‑150 (MQ=255)
                                                                                                                                  gTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTg                     <  1:939002/150‑1 (MQ=255)
                                                                                                                                   ttATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTgg                    <  1:725413/150‑1 (MQ=255)
                                                                                                                                   ctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTgg                    >  2:844530/2‑150 (MQ=255)
                                                                                                                                      tCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTc                 >  1:1122394/1‑150 (MQ=255)
                                                                                                                                       cGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCg                <  1:312649/150‑1 (MQ=255)
                                                                                                                                          cGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTa             >  2:1788704/1‑150 (MQ=255)
                                                                                                                                          cGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTa             <  1:1171426/150‑1 (MQ=255)
                                                                                                                                          cGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTa             <  2:2252527/150‑1 (MQ=255)
                                                                                                                                          cGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTa             <  1:1735630/150‑1 (MQ=255)
                                                                                                                                           gTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTAt            >  1:1790923/1‑150 (MQ=255)
                                                                                                                                             gctgcCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCACGGCGGTGTTTGAAATACCACTCGCCCAGCAATAACATCAGGGTccaattaa          >  1:1047632/1‑143 (MQ=255)
                                                                                                                                                gccgccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCAccc       >  2:1009242/1‑150 (MQ=255)
                                                                                                                                                gccgccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCAccc       >  2:653382/1‑150 (MQ=255)
                                                                                                                                                 ccgccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCTTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCTTTACATCTGGTTCTTTATCACCCt      <  2:547882/150‑1 (MQ=255)
                                                                                                                                                  cgccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATc                       >  1:461772/1‑132 (MQ=255)
                                                                                                                                                  cgccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATc                       <  2:461772/132‑1 (MQ=255)
                                                                                                                                                   gccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTTCCGTATTGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTAGCCCAGGCATTACATCTGGGTCGTTATCACCCttt    <  2:1294224/150‑1 (MQ=255)
                                                                                                                                                   gccCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCttt    <  2:1481506/150‑1 (MQ=255)
                                                                                                                                                    cccGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTa   >  2:482551/1‑150 (MQ=255)
                                                                                                                                                    cccGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTa   <  1:1895809/150‑1 (MQ=255)
                                                                                                                                                    cccGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTa   <  1:1565044/150‑1 (MQ=255)
                                                                                                                                                     ccGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTAg  <  1:1257654/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACCCTTTAG  >  NC_000913/1896328‑1896626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: