Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,602 2 bp→TT intergenic (‑80/‑130) ettA ← / → slt energy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,6020CT98.2% 159.4 / NA 54intergenic (‑80/‑131)ettA/sltenergy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/25);  total (29/25)
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.66e-01
*NC_0009134,630,6030CT98.2% 168.0 / ‑6.3 54intergenic (‑81/‑130)ettA/sltenergy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/24);  minor base G (0/1);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTT  >  NC_000913/4630362‑4630841
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                              
cccGCCATAATGCGCAGAAGGGTTGACTGACCCGCGCCGTTCAGACACAGGACACCAATTTTTGCCCCAGGGAAGAAACTTAGATAGATGTTTTTCAAAATTTGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTAGAAACGAATTGAGCCACGTTGGACTTCGCCTTTATGTTTATCGGGATTAGGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCATTTGGTCGATCa                                                                                                                                                                                                                                       <  2:245877/251‑1 (MQ=255)
           gcgcAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCTAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagt                                                                                                                                                                                                                             <  1:171732/250‑1 (MQ=255)
                  agGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGCAGTGTAGCCTTTTCCCCGCCTAATCCCAGTTGGTCGTTCACCTTCGCagtaaaagtaa                                                                                                                                                                                                                     >  1:32281/1‑250 (MQ=255)
                     gtagACTTACCCGCGCCGTTCAGACCCAGGTCTCCAATTTTTGCCCCAGGTAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCTACACGATGCATGGTATAAACGAATTTAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTaaaaa                                                                                                                                                                                                                  <  1:213717/248‑1 (MQ=255)
                            tACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTc                                                                                                                                                                                                            <  1:20169/250‑1 (MQ=255)
                                cgcgCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGACACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGCACAACTTTGCCGACACGATGCCTGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAACAAGTGTCCGTa                                                                                                                                                                                                        >  1:29461/1‑250 (MQ=255)
                                            aCCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAAc                                                                                                                                                                                           <  1:19305/251‑1 (MQ=255)
                                              ccAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACgg                                                                                                                                                                                         >  1:229106/1‑251 (MQ=255)
                                                       cAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGGTTCGGCGGAACAACTTTGCCGACACGATGCATGTTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCATTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTg                                                                                                                                                                                <  1:171618/251‑1 (MQ=255)
                                                        aaTTTTTTCCCCATGGAAGAAACGTAGAGAGATTGTTTTCAAAATATGACGTTTCGGCGGAACAACTGGGCCGTCGCGAGGCATGGTATAAACGAATTGAGCCTCGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGTCGAAGTGTAGACTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCGTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTAGCGTAAACAGCAATGaaggc                                                                                                                                                                               <  2:167533/251‑5 (MQ=255)
                                                        aaTTTTTGCCCCAGTGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGTATTGAGTCACGTTGGACTTCGCCTCTATGTTTAGTGTGATGATGAGTTTGCAAAGGCGGTGTGTAGCCTTTTTCTCCGCCTGATCTCAGTTGGTCGATCACCTACGCAGTAGAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTgg                                                                                                                                                                               <  1:126025/251‑1 (MQ=255)
                                                        aaTTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCCGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTgg                                                                                                                                                                               >  2:205972/1‑251 (MQ=255)
                                                            tttGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCTACACGATGCATGGTATAAACGAATTGAGCCACGTTGGTCTTCGCCTCTATGTTGATCGTGATGATGATTTTTCAAAGGCTAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGGCCGTAACGTGGCGTAAACGGCAATGACTGGTTa                                                                                                                                                                            <  1:243772/250‑1 (MQ=255)
                                                                    aGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATc                                                                                                                                                                   >  1:140203/1‑251 (MQ=255)
                                                                      ggAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGCGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCt                                                                                                                                                                  >  2:16623/1‑250 (MQ=255)
                                                                           aaaCTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAAGACTGGTTAGCATAAATCTATTAcgc                                                                                                                                                           >  1:76810/1‑251 (MQ=255)
                                                                                 agagagATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGACTTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCAt                                                                                                                                                      >  1:11688/1‑251 (MQ=255)
                                                                                   agagaTGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTCTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATg                                                                                                                                                     >  2:71070/1‑250 (MQ=255)
                                                                                        tGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGTTGCATGGTATAAACGCATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGCGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTg                                                                                                                                               >  2:181798/1‑251 (MQ=255)
                                                                                        tGTTTTTAAAAATATGAAGTTTCGGCTGAACAACTTTGCCGACATGATGCATAGTATAAACGAATTGAGCAACGTTGGACTTCGCCTCTATGTTTATCGAGATGATGGGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCATTTGGTCGATAAACTTCGCAGTAAACGTAAAAAATTGCCCGTAACGAGGCGTAAACGCCAAGGACTGGTAAACATAAATCTATTACGGGGCATGACgatg                                                                                                                                               >  2:12318/1‑248 (MQ=255)
                                                                                             ttCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAAGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTg                                                                                                                                          <  2:76810/250‑1 (MQ=255)
                                                                                               cAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGAt                                                                                                                                        <  1:71070/251‑1 (MQ=255)
                                                                                                       gACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTAGGTTTATCGTGATGATGAGTTTTAAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAACTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTac                                                                                                                                >  1:109958/1‑251 (MQ=255)
                                                                                                            ttCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACtt                                                                                                                            >  1:215121/1‑250 (MQ=255)
                                                                                                                             gCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTgg                                                                                                           >  1:51203/1‑250 (MQ=255)
                                                                                                                               cGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGaa                                                                                                         >  1:146931/1‑250 (MQ=255)
                                                                                                                               cGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAACATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGCTGTATTTACACTTATAGGATGCGCTTGTGGaa                                                                                                         >  2:153345/1‑250 (MQ=255)
                                                                                                                                   acGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATCAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGc                                                                                                    >  1:27585/1‑251 (MQ=255)
                                                                                                                                            ggTATAAACGAATTGAGACACGTTGGACTAGGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGGGTCCGTAACGTGGCGTAAACGGCAATGACTGGATAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCAAAACAAGt                                                                                           >  1:207464/1‑251 (MQ=255)
                                                                                                                                                 aaaCGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTAcc                                                                                       <  2:82361/250‑1 (MQ=255)
                                                                                                                                                               aCGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTGGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCGATTACGCGGCATGACGCTGCATTGATGTATTTACGCTTAGAGGATGCGCTTGTGGAAAAAGCCAAAGAAGTTACCTGTCGGCTGTTGGCt                                                                        <  2:29461/251‑1 (MQ=255)
                                                                                                                                                                cGTTGGACTTCGCCGCTATGTTTATTGTGTGGATGATTGTTCAAAGGCGAAGTGTAGCCTTTTACCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGTCATGACGCTGCATTTATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTg                                                                       <  2:201209/251‑1 (MQ=255)
                                                                                                                                                                 gTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGc                                                                      >  1:6373/1‑251 (MQ=255)
                                                                                                                                                                    ggACTTCGCCTCTATGTTTATCGTGATGATGATTTTTCAAAGGCGAAGAGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGACGCTCACCTTCGCAGTAACAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTAAGCGGCATGACGATGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGACgg                                                                   >  1:139391/1‑251 (MQ=255)
                                                                                                                                                                     gACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCACCGCCTAATCCCAGTTGGTAGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCAGTTGGCTGCCGgt                                                                  >  1:93742/1‑251 (MQ=255)
                                                                                                                                                                      aCTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGgtg                                                                 <  1:12318/251‑1 (MQ=255)
                                                                                                                                                                       cTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCTGGTCGATCACCTTCGCAGTACAGGTATAAAAGTGTCCGTCACGTGGCGTAACCCGCAATAACTGGTTCGCATCAATCTATTACGCGGCATGACGCTGCATTGCTGTATTTACACTTAGAGGATGCGCTTGTGGTAAAAGCCAAACCAGTTACCTGGCGGCTGCTCGCTGCCgtggt                                                                >  2:87023/1‑247 (MQ=255)
                                                                                                                                                                                 aTGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTctgct                                                       <  1:16623/250‑1 (MQ=255)
                                                                                                                                                                                   gTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTc                                                                                                                                                                                                            <  1:156831/99‑1 (MQ=255)
                                                                                                                                                                                   gTTTATCGTGATGATGAGTTTTCAAAGGAGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTc                                                                                                                                                                                                            >  2:156831/1‑99 (MQ=255)
                                                                                                                                                                                    tttATCGTTATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTACCGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCAGGTGTCTGTCTTCTGAc                                                    <  2:193749/250‑1 (MQ=255)
                                                                                                                                                                                    tttATCGTGATGGTGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCCGTTGGTCGATCACCTTCGCAGTAAAAGTCAAAAAGTGTCCGTAACGTGGCGTAAACGGCGATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTTCTTTGATGTATTTACCCTTAGAGCATGCGCTTGTGGAAAAATCCAAACAAGTTCCCTGGCGACTGTTGTCTGCCGTTGTCTGTCTGCTGACg                                                   >  2:130861/1‑251 (MQ=255)
                                                                                                                                                                                       aTCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTc                                                >  1:45717/1‑251 (MQ=255)
                                                                                                                                                                                             atgaTGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCagcagc                                          <  2:229106/248‑1 (MQ=255)
                                                                                                                                                                                                     ttttCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCCGTAATAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGTTGTCTGTCTGCTGACGGTCAGCAGCGTGGcgcg                                  >  2:72570/1‑251 (MQ=255)
                                                                                                                                                                                                        tCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCGGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAg                                <  2:140203/250‑1 (MQ=255)
                                                                                                                                                                                                                 aaGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAAGCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCGTGTGGAAAAAGCCAAACTAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTcagca                                            <  1:231458/229‑1 (MQ=255)
                                                                                                                                                                                                                 aaGTGTAGCCTTTTTCCCCGCCTAATCACAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGAGTAAACAGCAATGACTGGCTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTcagca                                            >  2:231458/1‑229 (MQ=255)
                                                                                                                                                                                                                      tAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGAt                  <  2:6373/250‑1 (MQ=255)
                                                                                                                                                                                                                             tttCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGagcagc           <  1:235056/250‑1 (MQ=255)
                                                                                                                                                                                                                             tttCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCg          >  1:196678/1‑251 (MQ=255)
                                                                                                                                                                                                                               tCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGt         <  1:72333/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ccGCCTAATCCCAGTTGGTCGAGCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTGAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGGTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCTCTGGATGAGCAGCGTAgcc     <  2:170305/251‑3 (MQ=255)
                                                                                                                                                                                                                                      cTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAACAAGAGTCCGTAACGTGGCGTAAAAGGCAATGACTGGCTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGACGtt  >  2:141140/1‑250 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                              
CCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTT  >  NC_000913/4630362‑4630841

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