Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913505,4400CA20.0% 143.5 / 12.6 62A176E (GCA→GAA) ushA5'‑nucleotidase/UDP‑sugar hydrolase
Reads supporting (aligned to +/- strand):  ref base C (40/8);  new base A (0/12);  total (41/21)
Fisher's exact test for biased strand distribution p-value = 9.00e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATG  >  NC_000913/505199‑505690
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                             
gATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTgg                                                                                                                                                                                                                                                      >  1:302638/1‑251 (MQ=255)
    ccGAAACTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAAc                                                                                                                                                                                                                                                 >  1:69900/1‑251 (MQ=255)
     cGAACCTGCTTTTCGCGGCATGAATCTGGTGGGCTATGACGCGATGGCGATAGGTAATCATGCATTTGATAATCCGCTCACCGTATTACGCCCGCAGGAAAAGTGGGCCAAGATCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAACCCGTGGGCGCTGGTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGCAACAGCAAAAATTGGTaccc                                                                                                                                                                                                                                                >  2:487723/1‑248 (MQ=255)
      gAACCTGCCTTTCGCGGTATGAATCTGGAGGGCTATGACGCGACCGCCATCCGTCATAATGACTTTGATCATCCGCCCACCGTATTACGCCAGCAGGAAAACCGGGCAAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTACACCGTGGGCGCTGTTTAACCGTCACGATCTGAAACTTGCCGGTATTGGGCTGCCAACCGATGACCCCGCACAACTTGGTAcccc                                                                                                                                                                                                                                               >  2:179440/1‑251 (MQ=255)
      gAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAAccc                                                                                                                                                                                                                                               <  1:453832/251‑1 (MQ=255)
        aCCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCgg                                                                                                                                                                                                                                            >  2:312321/1‑251 (MQ=255)
         ccTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTATCCCGGa                                                                                                                                                                                                                                           <  1:340387/251‑1 (MQ=255)
            gATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATa                                                                                                                                                                                                                                        <  1:321405/251‑1 (MQ=255)
                   gcgGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACt                                                                                                                                                                                                                                <  1:168539/251‑1 (MQ=255)
                          gAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGAGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCACCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGAAATACTTCACTGCTATCg                                                                                                                                                                                                                         >  2:505778/1‑251 (MQ=255)
                           aaTCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCg                                                                                                                                                                                                                         >  1:463420/1‑250 (MQ=255)
                              cTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAAtt                                                                                                                                                                                                                     >  2:307410/1‑251 (MQ=255)
                                        aTGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCCCCGTATTACGACAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATCTCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGCCACAGCAAAAATTGGTAACCCGGAATACTTCCCTGATATCGAATTTCGTAAGccc                                                                                                                                                                                                           <  2:193302/251‑1 (MQ=255)
                                        aTGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGccc                                                                                                                                                                                                           <  2:69900/251‑1 (MQ=255)
                                           ccgcgATGGCGATCGGTAATCATGAATTTTATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCccgcc                                                                                                                                                                                                        <  2:523011/250‑1 (MQ=255)
                                                    cGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATgaagcg                                                                                                                                                                                               >  2:500675/1‑251 (MQ=255)
                                                         ggTAATCATGAATTTGATAATCCGCTCACCGTATTCCGCCGGAAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAGGTACTGGCGAGCGCCTGTTAAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTAAAAATTGCCGTTATTGGGCTGACACCCGATGACACAGAAAAAATTGGTACCCCGAAATACTTCACTGAGATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCt                                                                                                                                                                                          <  1:98564/251‑1 (MQ=255)
                                                                    aTTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAgg                                                                                                                                                                               >  2:537135/1‑251 (MQ=255)
                                                                       tGATAATCCCCTCACCGTATTCCCCCAGCAGAAAATGGGGGCCAAGTCCCCTTTGCTTTCCGCGAATATCTCCAAAAAAAGTACTGGCGACCGCTTTTTAAAACCGTGGCCGCTGTTTAAGCGTCAGGATCTAAAAATTGCCTTTATGGGGCTGACACCCGATGACCCAAAAAAAATTGGAACCCCGAAATACTTCCCTGATATCGATTTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAg                                                                                                                                                                             <  1:332837/250‑1 (MQ=255)
                                                                                  tCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGCAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATAATTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGAGATTAAGGAGCTGCAacagac                                                                                                                                                                  >  1:37869/1‑250 (MQ=255)
                                                                                      ggtattccccaagcagAAAAAGGGGGCAAGGTCCCCGTTGCTTCCCGCGAATATCTACAAAAAAAGTACGGGGGAGCCCTTGTTAAAGCGGGGGGCGTTGTTAAAGCGACAGGATAAGAAAAGTGCGGTGATGGGGCGGAAACCGGAGGAAAAAAAAAAAATGGGAAACAAGAAAGACTCCAGGGATATGAAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGaaa                                                                                                                                                             <  2:127713/240‑1 (MQ=255)
                                                                                         ttacccacggaagaaagggTGGGCAAAGTTCGCTTTGTTTTCGGCGAATATCTGCAAAAAAAGGAGTGGCGGTCGCCTTTTAAAACGGTGGTCGTTGTTTAACCGACAGGATATGAAAATTGCCGTTTTGGGGCTAAACCCCGATGACACCAAAAAATTGGGAAACCCGAATTCCTTCTCGGATATCGAATTCCGAAAACCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAg                                                                                                                                                           <  2:379650/232‑1 (MQ=255)
                                                                                         attccGCCAGCAGAAAAAGTGGGCCAAGTCCCCGTTGCTTTCCGCGAATATCTACAAAAAAAGTACTGGCGAGCGCCTGTTAAACCCGTGGGCGTTTTTAAAGCGTAAGGATCTGAAAATTGCCGTTATTGGGCGGACACCCGATGACACAGAAAAAATGGGTACCCCGGAATCCTTCTCTAATATCGATTTTGGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGc                                                                                                                                                          <  1:63054/247‑1 (MQ=255)
                                                                                         aTTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCCGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAACAATTGGTACCCCGGCATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGTTGTTTCATGAGCTGCAACAGACAGTAAAGc                                                                                                                                                          >  2:109502/1‑251 (MQ=255)
                                                                                              gCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGAc                                                                                                                                                     >  1:309352/1‑251 (MQ=255)
                                                                                               ccAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGCAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGAc                                                                                                                                                     >  1:486077/1‑250 (MQ=255)
                                                                                                      aaaaaggggcccaggtccccGTTGTTTCCCGGGAATATCTACAAAAAAAGTACTGGCGAGCCCCTGTTTAACCCGGGGGCGTTTTTAAAGCGTCAGGATCTAAAAATTGCCGTTATGGGGCTGACAACCGATGACACAGAAAAAATTGGTACCCCGGAATCCTTCACTGATATGGATTTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATtattat                                                                                                                                             <  1:29709/234‑1 (MQ=255)
                                                                                                       aaaaGGGGGCCAAGTTCCCGTTGCTTCCCGCGAATATCTACCAAAAAAGTACGGGGGGGCGCTTTTTTAACCCGGGGGCGTTTTTAAGGCGTCAGGATCTAAAATTTGCCGTTATGGGGCTAACACCCGAGGACACAGAAAAAATTGGTACCCCGGAATCCTCCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATc                                                                                                                                            <  1:19158/251‑1 (MQ=255)
                                                                                                           gTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcgg                                                                                                                                        >  2:228496/1‑251 (MQ=255)
                                                                                                            ggggCCAAGTCCCCGTTGTTTCCCGCGAATATCCACCAAAAAAGTACGGGCGGGCCCTTGTTAAACCCGGGGGCGCTGTTTAGGCGGAAGGATCTAAAAATTGCCGTTATGGGGCTGACACCCGATGACACAGCAAAAATTGGTACCCCGAAATACTCCACTGATATCGAATTTCGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcggc                                                                                                                                       <  1:114535/250‑1 (MQ=255)
                                                                                                             gggCCAATTTCCCGTTGCTTCCCGCAAATATCAACCAAAAAAGTCCTGGCGGGCGCTTGTTTAAACCGGGGGCGTTGTTTAGGCGTCAGTATCTAAAAATGGCCGTTTTGGGGCTGAAAACCAATGACACAGAAAAAATTGGTACCCGGAATTTCTCCACTAATAGCGATTTTCTAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCgcggcg                                                                                                                                      <  1:178860/251‑1 (MQ=255)
                                                                                                                ccAAGTGCCCATTGCTTTCAGCCAATATCTACCCGCAACGTACTGGCGAGAGCCTGTTTAAACCGTGGCCGCTGCTTAAGCGTCAGGACCTGAAAATTGCCGTTATTGGGCTGCCAACCGATGAAACAGCCAAAATTGGTAACCCGGAATACTTCACTGATATCGCATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATAGCGGCGAcc                                                                                                                                   >  2:443817/1‑251 (MQ=255)
                                                                                                                 caattgccctttGCTTCCCGAGAATATCTACCAAAAAAGTCCTGGGGAGCGCTTGTTTAACCCGGGGGCGCTGTTTAAGCGTCAGTATCAAAAAATTCCCGTTATTGGGCTGACAACCGATGACACAAAAAAAATTGGTACCCCGGAATCCTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGAccc                                                                                                                                  <  1:300187/241‑1 (MQ=255)
                                                                                                                  aaGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGAccc                                                                                                                                  >  1:66209/1‑250 (MQ=255)
                                                                                                                        ccGTTGCTTCCCGCAAATTTTTCCCAAAAAAGTACTGGGGGGCGCTGTTTTAAACGGGGGGCTCTTTTTAAGCGTCAGGATCTAAAATTGGCCGTTATTGGGATGACAACCGATGACCCACAAAAAATTGGTACCCCGAAATACTTCCCTGAAAACGAATTTGGAAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATggg                                                                                                                           <  1:352387/251‑1 (MQ=255)
                                                                                                                              cTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGTCAtt                                                                                                                      >  1:453123/1‑250 (MQ=255)
                                                                                                                                        atatCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCAATATGGGGCATTACGATAATggg                                                                                                           >  2:292303/1‑250 (MQ=255)
                                                                                                                                           tCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTg                                                                                                          >  1:31801/1‑249 (MQ=255)
                                                                                                                                                    aaaGAACTGGCGAGCGACTATTTAAAACGTGGGCGCTGATTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAACATTGGGAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGAAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCTTTACGATAATGGTGAGCACGGctct                                                                                               >  2:409281/1‑251 (MQ=255)
                                                                                                                                                        tACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACg                                                                                           >  1:272311/1‑251 (MQ=255)
                                                                                                                                                         aCTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGc                                                                                          >  1:47047/1‑251 (MQ=255)
                                                                                                                                                             gCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCg                                                                                      >  1:289660/1‑251 (MQ=255)
                                                                                                                                                             gCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCg                                                                                      >  1:31718/1‑251 (MQ=255)
                                                                                                                                                                   gcCTATTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAt                                                                                >  2:58152/1‑251 (MQ=255)
                                                                                                                                                                   ccctttttAAACCGTGGGCGCTGTTAAAGCGTGAGTATCAAAAATTTCCCGTTATTGGCCTGACACCCGATGACACAGAAAAATTTGCAACCCCGAAATACTCCACTGATATCAATTTTCGTAGGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAt                                                                                <  1:312321/246‑1 (MQ=255)
                                                                                                                                                                              aCCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGCTGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCCGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTTAGCACTGCTCTAACGCACCGGACGCTGTGgagaggg                                                                      >  2:306892/1‑248 (MQ=255)
                                                                                                                                                                                        cTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGCTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                                            >  1:424094/1‑251 (MQ=255)
                                                                                                                                                                                           tttAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                                                         >  2:83316/1‑251 (MQ=255)
                                                                                                                                                                                                  cGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCgg                                                   >  2:198514/1‑251 (MQ=255)
                                                                                                                                                                                                      aGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCAAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGTTGGCGCGCGCGCTGCCTGCCGGATCg                                               >  1:273434/1‑251 (MQ=255)
                                                                                                                                                                                                              aaaaTTGCTTTTATGGGGCTGACACCCGATGACACAGAAAAAATTGTAACCCCGAAATTCTTAACGGATATCGTATTTGGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat                                       <  1:307410/251‑1 (MQ=255)
                                                                                                                                                                                                                  atggacGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTAACTGATATCGAATTTCGTAAGCCCGCCGATGAAGAGAAGCTGGCGATTCAGGAGCTGCAACAGACAGAAACGCCAGACATTATTATCGCGGCGACACATATGGGGCATTACGATAATGGTGAGCACTGCTCTAACGCACCGGGCGATGTGGAGATGGAACGCGCTCTTCCTGCCGGAGTGCTGGCgatgat                                    >  2:255252/6‑251 (MQ=255)
                                                                                                                                                                                                                    gCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt                                 >  1:297622/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgtc                                >  1:305792/1‑251 (MQ=255)
                                                                                                                                                                                                                        gTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCgg                              >  1:467623/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtgg                           >  2:161791/1‑251 (MQ=255)
                                                                                                                                                                                                                              gggCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTc                         >  1:422846/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gaaacccGATGACACAGAAAAAATTGATACCCCGAATTACTTCTCTGATATAGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTc                     <  1:79348/244‑1 (MQ=255)
                                                                                                                                                                                                                                   gACAACCGATGACACAGGAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCAAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGCGGTCCCTCGc                   >  1:100505/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGc                   >  1:34826/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aacaccgaatgcaccacAAAAATTGGGTACCCCGAAATCCTTCATGGATATCTAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTc                     <  1:114130/233‑1 (MQ=255)
                                                                                                                                                                                                                                       aCCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGCATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCa                  >  1:435287/1‑248 (MQ=255)
                                                                                                                                                                                                                                               acacAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTc        >  1:439326/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cAAAAATTGGTAACCCGGAATACTTCACTGCTATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATg  >  2:204071/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCAGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATg  >  2:126724/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                             
GATGCCGAACCTGATTTTCGCGGTATGAATCTGGTGGGCTATGACGCGATGGCGATCGGTAATCATGAATTTGATAATCCGCTCACCGTATTACGCCAGCAGGAAAAGTGGGCCAAGTTCCCGTTGCTTTCCGCGAATATCTACCAGAAAAGTACTGGCGAGCGCCTGTTTAAACCGTGGGCGCTGTTTAAGCGTCAGGATCTGAAAATTGCCGTTATTGGGCTGACAACCGATGACACAGCAAAAATTGGTAACCCGGAATACTTCACTGATATCGAATTTCGTAAGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATG  >  NC_000913/505199‑505690

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: