Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,2990G.11.1% 139.5 / 9.9 45coding (931/1485 nt)gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base . (4/1);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected: Not supported by required number of reads on both strands.

GATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGGAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAG  >  NC_000913/3963002‑3963598
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
gaTCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGGAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAACGTACGGCt                                                                                                                                                                                                                                                                                                           >  2:333007/1‑301 (MQ=255)
    gACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGGAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGa                                                                                                                                                                                                                                                                                                      >  2:487874/1‑301 (MQ=255)
      gtggAGTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGTAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgtaa                                                                                                                                                                                                                                                                                                     <  1:388246/299‑4 (MQ=255)
        ggATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGTAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTACGCTGCgtaaaa                                                                                                                                                                                                                                                                                                   >  2:388760/1‑296 (MQ=255)
                                              tttCTTTGCGGGGGTAAAACCGGGAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAgg                                                                                                                                                                                                                                                            <  1:525452/301‑1 (MQ=255)
                                                                               tttACGTACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                                                                                                                                                                                           <  2:43136/300‑1 (MQ=255)
                                                                               cgatctcACCAGATAAGCATCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                                                                                                                                                                                           <  2:47290/295‑1 (MQ=255)
                                                                                   cGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTCATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATCCCGTAAACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGACAAacacc                                                                                                                                                                                                                       >  1:308771/1‑300 (MQ=255)
                                                                                     cACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGACTTCACGCCGCGCACCGCCCGACAGGGGCATacaaacg                                                                                                                                                                                                                     >  1:47237/1‑297 (MQ=255)
                                                                                           ttAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTc                                                                                                                                                                                                               <  2:309175/300‑1 (MQ=255)
                                                                                           atagtgggggtttgggtgggtttGTTTGAAGTCAACGTTCAGGGCGTTTTCATGAAGCTGACAGGGGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTc                                                                                                                                                                                                               <  1:502399/280‑1 (MQ=255)
                                                                                                                       aaGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCg                                                                                                                                                                                   <  2:458778/301‑1 (MQ=255)
                                                                                                                            aaCGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGTGGCTTATTCCTTCAAGAGGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTCCGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt                                                                                                                                                                              <  2:463101/301‑3 (MQ=255)
                                                                                                                               ccTCAGGCCGTTTTCATGAAGCTGACAGGGGCTGATGAGAAAATGGTGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCGACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCattttt                                                                                                                                                                           <  1:332602/300‑6 (MQ=255)
                                                                                                                                   ttccgttcttCATGAAGTTGACTGGCGTTGATTAGCAAATCGCGGCTTATTGCTTCAATAGGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatttttaccg                                                                                                                                                                       <  2:154861/293‑10 (MQ=255)
                                                                                                                                             cATGAAGTTGACAGGCGCTGATGAGCAAATCTTGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGTTCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTaa                                                                                                                                                             <  2:452307/301‑1 (MQ=255)
                                                                                                                                             cATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGAACAGGATCGCCAGACCAATAGGGAACACTaa                                                                                                                                                             >  2:627886/1‑301 (MQ=255)
                                                                                                                                                          ggCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATCTTCAGTTCGGTAAAAATGGCGATCAGGATCCCCAGACCACTCGGGAACACTAACGCACGTTCCAGc                                                                                                                                                >  1:500022/1‑301 (MQ=255)
                                                                                                                                                                                cTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCCCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAACTGGCGATCAGGATCGCCAGACCACACGGGAACACAAAAGCAAGTTCCAGCACCACCCCATCAATCCACAGtt                                                                                                                          >  1:395922/1‑301 (MQ=255)
                                                                                                                                                                                  tATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATCTCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAAAACTAACGCACGTTCCAGCGTCAGCCCGGCCATCCCCAGTtcc                                                                                                                        >  1:360744/1‑301 (MQ=255)
                                                                                                                                                                                      gttttaaGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGGGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATAGTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCa                                                                                                                    <  1:487228/296‑1 (MQ=255)
                                                                                                                                                                                            aGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGAc                                                                                                              <  2:500669/301‑1 (MQ=255)
                                                                                                                                                                                                    ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCaatga                                                                                                      >  2:113595/1‑301 (MQ=255)
                                                                                                                                                                                                              tcttGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCt                                                                                            <  2:284481/298‑1 (MQ=255)
                                                                                                                                                                                                                 tGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGAACGCAAAGAAAGGCCCGccaaa                                                                                         >  1:284178/1‑297 (MQ=255)
                                                                                                                                                                                                                       aagaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACACGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAAtt                                                                                   <  2:396436/301‑1 (MQ=255)
                                                                                                                                                                                                                                  cGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACttt                                                                        >  2:55825/1‑301 (MQ=255)
                                                                                                                                                                                                                                    cAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTc                                                                      >  2:47095/1‑301 (MQ=255)
                                                                                                                                                                                                                                          tGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGt                                                                >  2:406202/1‑301 (MQ=255)
                                                                                                                                                                                                                                           gggtACGTGCTGTGCCTGTTCAATATCGATCATAAAGGGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTa                                                               <  1:55757/298‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTACAGCAACTTTTCCAGGGTAATTCGTTc                                                       >  1:520917/1‑301 (MQ=255)
                                                                                                                                                                                                                                                      ttgcctggttatgtctggtCATAAAGCGTCGTTGAATATTACTCGTCGTATGGCTGCGAATAGCCTGTTCGACGGCCAGAGGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTGCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatc                                                    <  1:627009/283‑1 (MQ=255)
                                                                                                                                                                                                                                                       tgCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatcc                                                   >  2:597790/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          cctacaggacgctttgccgtgttggctTTGAATATTCCGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACAGGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatccatc                                                <  1:19589/278‑1 (MQ=255)
                                                                                                                                                                                                                                                            gATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCCGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCCCGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTTTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATccc                                              >  1:349032/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              tCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCt                                            >  2:345363/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 atatCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAAGGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCtgtg                                         <  1:405639/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCAATCATGATTTc                            >  2:583473/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 ggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTcatccatcagatcggaagagcacacgtctga                         >  2:19640/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                                 ggcgTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTg                         <  2:361228/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    agcaggatattcGCAGCCGTACGCTGCGTAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTACAGGGTAATGCGTTCATCAAACCCCTGTGCCATCATGATTTCCTGTa                       >  2:128667/12‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     gcaggatattcGCAGCCGTACGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCACTTGCTGCCACTTTTCCAGGGTCATGCATTCATCCATCCCCCTTGCCATCATTATTTCCTTTAAc                     >  1:366093/11‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           tttACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGc               <  2:349495/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            attcGCAGCCGTACGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCCACTTTTCCAGGGTAATTCGTTCATCCATCCCCTGAGCAATCATGATTTCCTGaaaaccatgca              >  2:462773/4‑293 (MQ=255)
                                                                                                                                                                                                                                                                                                         cGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGaaa  >  2:378097/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                            
GATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTGTGCGGGGGTAAAACCGGGAAGATCCAGATTACGCACCAGATAAGCAGCGTGTTGCGGCGCTTGTTTGAAGTCAACGCTCAGGCCGATTTCATGAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAACACTAACGCACGTTCCAGCGTCAGCCCGTCAATCTCCAGTTCCTCCAGCCGACCGCAATGAATGGCTCGCTGTTTCAATTGCTGCAACTTTTCCAGGGTAATGCGTTCATCCATCCCCTGTGCCATCATGATTTCCTGTAACGCCTGCACGGTGCCGGAAG  >  NC_000913/3963002‑3963598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: