Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,894,839 C→T intergenic (+68/+80) ygcP → / ← ygcQ putative anti‑terminator regulatory protein/putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,894,8390CT98.0% 149.9 / NA 49intergenic (+68/+80)ygcP/ygcQputative anti‑terminator regulatory protein/putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/22);  minor base A (1/0);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGC  >  NC_000913/2894593‑2895074
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
cccTGCCAGGCTGTATGCCCAAAGGGCGGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGa                                                                                                                                                                                                                                         <  2:373404/250‑1 (MQ=255)
  cTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCAGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGGGTCGTGGCGCTTTCCACCACGAATACCGGGGCCTGGACGTTAGCGAAACAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACTCCACATCCGACATTCAGCTACTGATg                                                                                                                                                                                                                                       >  1:257296/1‑251 (MQ=255)
   tGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTTATGc                                                                                                                                                                                                                                      >  1:459842/1‑251 (MQ=255)
                    aaaGGGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGGGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTGGTAGGCCGCACGCCACATCCGACATTCAGCGTCTTATGCGACGCTTGACGCGTCtt                                                                                                                                                                                                                     <  2:459842/251‑1 (MQ=255)
                           tgggctgggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCCGCGTCTGATGCTACGCTTGACGCGTCTTATCAGGc                                                                                                                                                                                                              >  2:413363/1‑251 (MQ=255)
                                                       aCCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTa                                                                                                                                                                                                                    <  1:167580/217‑1 (MQ=255)
                                                       aCCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTa                                                                                                                                                                                                                    >  2:167580/1‑217 (MQ=255)
                                                             gATTGCTGGTGGGCTGGTGTGCGATGAAGGAGATGCTCATAATGAGATTAACGCGGGTTTCGGGTCGCTTTCCACAACGAATACCGGGGTCGGGATGTTAGCGAAAAAATTACTTTGACGGGAGAAGTGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTTATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                                                                                                                                                                             <  1:485341/250‑1 (MQ=255)
                                                              aTTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatc                                                                                                                                                                           >  1:388252/1‑251 (MQ=255)
                                                                tGCCGGGGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatca                                                                                                                                                                          <  1:218176/250‑1 (MQ=255)
                                                                           tGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGa                                                                                                                                                                                          <  1:213628/223‑1 (MQ=255)
                                                                           tGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGa                                                                                                                                                                                          >  2:213628/1‑223 (MQ=255)
                                                                               ggactaTGAAAAATATGCGCGTAATGCGGTGAACGAGGGGTTCGTGGCGCTTTCCCTCACGAATACTGGGGACTGGTCGTGTGCGAAAAAATTTCTTTGACGAGATAATCGTAACCAATGGAATTTTGTTTGATTGGGAGGCCGATCGCCACACCCGACATTCAGCGTCTGATGAGACGCTTGACGCGTCTGATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCAGCCGGCAATACGTTa                                                                                                                                                           <  2:9693/245‑1 (MQ=255)
                                                                                  cGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTgg                                                                                                                                                       >  2:349210/1‑251 (MQ=255)
                                                                                         gaagaTGCGCGTAATTCGATTAACGCTGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTg                                                                                                                                                <  1:414952/251‑1 (MQ=255)
                                                                                           agaTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCACTCTGCg                                                                                                                                              >  2:154189/1‑251 (MQ=255)
                                                                                             aTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTACCACCACGAATCCCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                                                                                            <  1:139073/251‑1 (MQ=255)
                                                                                                  cgTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGt                                                                                                                                       <  1:402969/251‑1 (MQ=255)
                                                                                                    tAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCa                                                                                                                                                                                                 >  2:18133/1‑191 (MQ=255)
                                                                                                    tAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCa                                                                                                                                                                                                 <  1:18133/191‑1 (MQ=255)
                                                                                                        gCGATTAACGCGGGTGTCGTGGCGTTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCa                                                                                                                                 <  2:388252/251‑1 (MQ=255)
                                                                                                                     gtgtCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTcc                                                                                                                      >  1:362223/1‑249 (MQ=255)
                                                                                                                                cTTTCCACCACGAATCCCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAAGCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGt                                                                                                                                             <  1:482485/215‑1 (MQ=255)
                                                                                                                                            aatacagaggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTCGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAAcc                                                                                              >  2:441840/9‑250 (MQ=255)
                                                                                                                                                 cGGGGTCTGGACGTTAGCGAAAAAATTCCTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCCGCACCACTTTCCACGCATCCACCACGccc                                                                                        >  2:72063/1‑250 (MQ=255)
                                                                                                                                                    ggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGGTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGAGGCGACGCTGGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGAGCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCTCTTTCCAGTCATCAACCACGCCAACa                                                                                     <  1:441840/251‑1 (MQ=255)
                                                                                                                                                                  cGAAAAAATAACTTTGACGGGCTGATAGTCACCCTTTGAATTTGGTTTGATTTGTACTCCGCCCGACACATCCGGCATTCAACGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCACAGGCCAGACCCGTAGGCAGGTACAACAGGCAACACGAAATTGGCAATCTGCGTGAGTGTTTGTCACCCGTGCCCAAAGCACAACTTTCCCGTCGTCAACCACGcaaagctcagcccgtggca                                                                        >  2:67077/1‑232 (MQ=255)
                                                                                                                                                                      aaaaTTTCTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATAcc                                                                   <  2:362223/251‑1 (MQ=255)
                                                                                                                                                                       acgttACTTTGACGGGATAATAGTAACCAATTGCATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGCACCGTAGGCAGGCTCATCCGGCAACACGTTATTGGCAATCTGCGTGAACGTTTGTCCCCAGAGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGCAAAtcccc                                                                  >  2:72116/4‑247 (MQ=255)
                                                                                                                                                                        aaTTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACACGTCCCCAGGCCAGACCCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGc                                                                 >  2:111319/1‑251 (MQ=255)
                                                                                                                                                                         aTTACTTTGACGGGATAATCTAAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACAGCTGCCTGTGAAAATACCGCa                                                                <  2:43973/251‑1 (MQ=255)
                                                                                                                                                                                    gggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTc                                                     >  1:410350/1‑251 (MQ=255)
                                                                                                                                                                                      gATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCCGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGCCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGAACCACTTTCCAGTCATCAACCACGCCCACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTcct                                                   >  1:458677/1‑249 (MQ=255)
                                                                                                                                                                                             gTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAAt                                            <  1:99031/251‑1 (MQ=255)
                                                                                                                                                                                               aaCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATc                                           <  1:349210/250‑1 (MQ=255)
                                                                                                                                                                                                   aatggAATTGGGTTTAATTTGTAGGCCGCACGCCAAACCCGACATTCAGCGTCTGAGGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGccac                                      <  2:230517/247‑1 (MQ=255)
                                                                                                                                                                                                    aTTGAATTTAATTTGATTTGTAGGCCGCACGCCACATCCGACATTAAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGACCCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCGCCCCTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGccacc                                     >  2:284109/1‑251 (MQ=255)
                                                                                                                                                                                                      tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCcaccac                                   >  2:40770/1‑251 (MQ=255)
                                                                                                                                                                                                        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCACCGCTGGCGTCATGATTAATCTCCACcaca                                  >  1:250401/1‑250 (MQ=255)
                                                                                                                                                                                                              gTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTACAGCACCACTTTCCAGTCATCAACCACGCCAACATCTCCCTGTAACCATACCTCCGCGCTGACGTCATGATTAATCCCCACCACAAACTTGc                           >  2:124433/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTg                      >  2:412645/1‑251 (MQ=255)
                                                                                                                                                                                                                     ttGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCACCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCg                    >  2:442506/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgc               >  1:43766/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggCCGCACGCCACAGCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGACATGATGAATCGCCACCACAAACTTGCTATTGCGTAcgc               <  1:284109/251‑1 (MQ=255)
                                                                                                                                                                                                                             cccACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATGAATCGCCACCACAAACTTGCTATTGCGTAcgcccg            <  2:275309/249‑1 (MQ=255)
                                                                                                                                                                                                                             cGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGACCCGTAGGCAGGCTCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCACCCCCACAACCTTGCTATTGCGTAcgcccc            >  2:451062/1‑250 (MQ=255)
                                                                                                                                                                                                                               cACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTc                                                                                                                                  >  1:318079/1‑131 (MQ=255)
                                                                                                                                                                                                                               cACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTc                                                                                                                                  <  2:318079/131‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCCTTAACgc  >  1:180324/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
TCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGC  >  NC_000913/2894593‑2895074

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