Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,044 A→G Q50Q (CAA→CAG serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,0440AG91.4% 96.1 / ‑1.1 37Q50Q (CAA→CAGserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (16/16);  total (17/20)
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.79e-01

TTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCC  >  NC_000913/4624808‑4625287
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
ttGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAAt                                                                                                                                                                                                                                        >  1:69716/1‑250 (MQ=255)
         aaTGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCa                                                                                                                                                                                                                               >  1:250709/1‑250 (MQ=255)
          aTGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCaaa                                                                                                                                                                                                                             <  1:373387/251‑1 (MQ=255)
               gCCAGAAGCAATGGATACAAGGGAGCCGCATGCGTTATTTTCCCTGCTTCGAACGATTTTTCAGGAGCCTTAATGCCTAACATTTCCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTgg                                                                                                                                                                                                                         <  2:123986/250‑1 (MQ=255)
                  aGAAGCAATGGATACAAGGTAGCATCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATTCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGATGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTgg                                                                                                                                                                                                                     >  2:376953/1‑250 (MQ=255)
                                  aGGTAGCCTCATTCGCTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGTCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGAttgtt                                                                                                                                                                                                     <  1:4914/251‑1 (MQ=255)
                                  aGGTAGCCTCATGCGTTATTTTCCCTGACTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACAAGAGCAAACTGGGTGCGGCGATGGTGAttgtt                                                                                                                                                                                                     >  1:309461/1‑251 (MQ=255)
                                           tcagCGTTAGTTTCCCTGTTTCTGACGAGTTTACAGTAGCCTTAATGCCTACCAGTACCCGTTGCTACCTGCCTGAAGATGACTCTTTTTGGCCGGGTCTGCCTGTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTTGTGGGCTGGATAAACAGTGTCTTACCCAATACCATAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCt                                                                                                                                                                                              <  2:203831/246‑1 (MQ=255)
                                           cATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTg                                                                                                                                                                                             <  1:364485/250‑1 (MQ=255)
                                                                    cGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGTCTATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTc                                                                                                                                                                   <  1:338019/251‑1 (MQ=255)
                                                                    cGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTc                                                                                                                                                                   >  1:274414/1‑251 (MQ=255)
                                                                    cGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTc                                                                                                                                                                   <  1:209647/251‑1 (MQ=255)
                                                                                 gCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGATGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTtc                                                                                                                                                      >  1:142139/1‑251 (MQ=255)
                                                                                   cTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTTGTGCCGCCTGGTGCGAGGGAGAATATCAGGGGGTTCGTCTGGCAGGtttac                                                                                                                                                    >  2:273034/1‑248 (MQ=255)
                                                                                              aCATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACgg                                                                                                                                         >  1:176776/1‑251 (MQ=255)
                                                                                               cATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACggg                                                                                                                                        >  1:210915/1‑251 (MQ=255)
                                                                                                aTTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGc                                                                                                                                       >  1:314099/1‑251 (MQ=255)
                                                                                                      tGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTacgcg                                                                                                                                 >  1:124079/1‑248 (MQ=255)
                                                                                                        gTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTACTTAAGTGGTGATGAAGAGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGATGGATAAACAGCGTCTGACCCAATACCAGAGCAAACAGGGTGAGGCGATGGTGATTGTTGCCACCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTacgcgac                                                                                                                               >  1:321568/1‑246 (MQ=255)
                                                                                                                                 ttATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGa                                                                                                         <  1:24855/248‑1 (MQ=255)
                                                                                                                                  tATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGt                                                                                                      <  1:30053/250‑1 (MQ=255)
                                                                                                                                   aTGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCg                                                                                                    >  1:52430/1‑251 (MQ=255)
                                                                                                                                   aTGGCCGGGTCTGCCTCTTTCATTAAGGGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACAGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCg                                                                                                    >  1:52472/1‑251 (MQ=255)
                                                                                                                                     ggCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTGCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGTCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGcc                                                                                                  <  2:210915/251‑1 (MQ=255)
                                                                                                                                          ggTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCt                                                                                             <  1:117698/251‑1 (MQ=255)
                                                                                                                                                                  gatgaAGTGATGCCCCTGGAGTGCCACGCAGGTCGTATCGGCCGGCTGCTGAATGTGCGGGGGCTGGATAATCAGCGTCTGACCCAATTCCAGAGCAATCTGGGTGCGGCGAAGGTGATTGTTGCCTCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTTACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGAAGTCGCCCCGCTGTGGAAATTCCCGCACCGGCGCAc                                                                     <  2:363087/251‑1 (MQ=255)
                                                                                                                                                                                    gATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGTGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCTCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGAt                                                   <  1:236778/251‑1 (MQ=255)
                                                                                                                                                                                       tACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGGGATgg                                                 >  1:177047/1‑250 (MQ=255)
                                                                                                                                                                                        aCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGTCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAGCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGa                                                <  1:351315/250‑1 (MQ=255)
                                                                                                                                                                                        aCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGac                                               <  1:376953/251‑2 (MQ=255)
                                                                                                                                                                                           ccgacggtctttgCGGCTGGCTGCTGTATGGTAGTGGGATGGATAAACATCGGCTGACCCAATACCAGAGCAAACTGGGGGCGGCGAGGGTGATGGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTGCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                            <  2:142139/239‑1 (MQ=255)
                                                                                                                                                                                            cGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATgg                                           >  1:10222/1‑251 (MQ=255)
                                                                                                                                                                                                  ccGTAGCGGCTGGCTGCTGAATGGTCGTGGGCTGGATAAACATCGTCTGACCCAATCCCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCGGGTGCGTGAAAGATTAGCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGACACGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTTCTGGTGATGGACATGGATTCCa                                     <  2:76472/250‑1 (MQ=255)
                                                                                                                                                                                                    gTAGCGGCTGTCTGCTGTATGGTCGTGGGCGGGATAAACAGCGGTTGCCCCCATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                   <  2:314099/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGTCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGAGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGCTGGAAATGGACTCCACCGCCATCCAGATTGAATGTATtggtgg         >  2:163712/1‑247 (MQ=255)
                                                                                                                                                                                                                              gggCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGAGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtgatg          <  1:163712/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aTAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc  <  1:284465/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
TTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCC  >  NC_000913/4624808‑4625287

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