Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,158,653 T→G 18.6% V55G (GTG→GGG)  oxyR → oxidative and nitrosative stress transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,158,6530TG18.6% 143.3 / 3.2 81V55G (GTG→GGG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/38);  new base G (14/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 2.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATAATCGTTCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAC  >  NC_000913/4158429‑4158875
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
aTAATCGTTCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGCCCCGCCGTAAA‑Gggtt                                                                                                                                                                                                                              >  1:281099/1‑228 (MQ=255)
aTAATCGTTCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑gtgct                                                                                                                                                                                                                              >  1:315384/1‑226 (MQ=255)
   aTCGTTCATTGCTATTCTACCTATCGCAATGAACTATCGTGGCGATGGAGGATGGATACTGAATATTCGTGTTCTTGAGTACCTGGTGGCATTGGCTGCACACCGCCATTTTCGGCGTGCGGCAGAGTCGTGCCCCGTTAGCCAGCCGACGCTTAGCGGGC‑AATTAGTAAACATGTAAGCTGAGCTGGGCGTGATGTTGGGGGAGCGGCACAGGCGGAAA‑GTgtggtt                                                                                                                                                                                                                           >  2:481022/1‑228 (MQ=255)
        ttaTTGCTATTCTAGCTATCGCCGTGAGCTATCGTGGCGATGGTGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAAAACCTCCATTTTCGGCGTGCGGCGGATTCCTGCAACGTTAGCCAGCCGACGCTTAGCGGGTAAATTCGTAAGC‑AGGAAGTAGAGCTGCGCGAAAGGTTGGGGGAGCGGACCAGCCGTCAA‑TGGTGGTTccac                                                                                                                                                                                                                       >  2:481389/3‑224 (MQ=255)
        tCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTGGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCAccc                                                                                                                                                                                                                      >  1:230743/1‑228 (MQ=255)
        tCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATATTCGGCGTCCGGCAGATTCCTGCCACGTTAGCCAGCCGCCGCTTATCGGGGACATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTGGCGGGAGCGGAACCGCGCTAAATGTTTTGTTcccc                                                                                                                                                                                                                       >  2:510971/1‑228 (MQ=255)
         cATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCGTATCGGGCAAATTCGTAAGC‑TGGCAGATGAGCTGGGCGTGCTGTTGCTGGAGCGGACCACCCGTAAA‑GGGTTGTTACCCCa                                                                                                                                                                                                                     >  2:205467/1‑228 (MQ=255)
                 ttCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGCCCCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTAACCCAGGCGGGaa                                                                                                                                                                                                             >  1:495669/1‑228 (MQ=255)
                      ccTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGCTGAGCGTGGCGTGATGTTGCTGGGGCGGACCACCCGTAAA‑GGGTAGTTCACCCCGGGGGGATtggtg                                                                                                                                                                                                        >  2:559386/1‑228 (MQ=255)
                      ccTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GGGTTGTTAACCCAGGGGGGAAtgctg                                                                                                                                                                                                        >  1:521038/1‑228 (MQ=255)
                      ccTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAACC‑TGGAAGATTAGCTGGGCGTGATGTTGCTGGAACGGCCCAGCCGGGAA‑GGGTTGTT‑ACCCCCGCGGGaaaagatg                                                                                                                                                                                                       >  2:457974/1‑222 (MQ=255)
                        tATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTTGCTGGAGCGGACCAACCGTAAA‑GTGTTGTTCACCCAGGCGGGAATgctgct                                                                                                                                                                                                      >  1:97644/1‑228 (MQ=255)
                         aTCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGGGATGGGCGGGATGTTGCTGGAGCGTACCAGCCGTAAA‑GTGTTGTTCACCCAGGGGGGAAAGcggctg                                                                                                                                                                                                     >  2:160495/1‑228 (MQ=255)
                         aTCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAGTGttgctg                                                                                                                                                                                                     >  1:588696/1‑228 (MQ=255)
                         aTCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACCCTTAGCGGGCCAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGCTGTTGCTGTGTCGGGCCAGCCCTAAA‑GTGGTTTCCACACAGGCGGGAATGctgctg                                                                                                                                                                                                     >  2:562265/1‑228 (MQ=255)
                           cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGGGCTGGGCGGGATGTTGCTGGAGCGGACCAGCCGTAAA‑GGGTTGTTCAACCAGGCGGGAAAGGtggtggg                                                                                                                                                                                                   >  2:439939/1‑227 (MQ=255)
                           cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTAACCCAGGCGGGACTGCTGCtggt                                                                                                                                                                                                   >  1:274445/1‑228 (MQ=255)
                           cGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCGGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGCTGAGCTGGGCGGGATGTTGCTGGAGCGGACCGGCCGTAAA‑GCGTGTTTCACCAGGGCGGGGATGTTGCtggg                                                                                                                                                                                                   >  2:407064/1‑227 (MQ=255)
                              cATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTACA‑GGGTTGTTCACCCAGGCGGGAATGCtggtggtgg                                                                                                                                                                                                 >  2:573263/1‑227 (MQ=255)
                                 gaccttgcGTGGGGATGGAGGATGGATAATGAAATTTCGTGATCTGGAGTACCGGGTGGCATTGGCTGACCACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGAt                                                                                                                                                                                               <  1:371592/219‑1 (MQ=255)
                                    ctttCGGGGCGATGGAGGATGATTAATGAATATTCGTGCTTTTGGTTCCTTGGTGGCATTGGCTGACCACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGgc                                                                                                                                                                                          <  1:439939/225‑1 (MQ=255)
                                       tCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTTGCTGGAGCGGGCCACCCGTAAA‑GGGGTGTTCACCCAGGCGGGAAGTGTGCTGGGGGACAAGgcggg                                                                                                                                                                                       >  2:371592/1‑226 (MQ=255)
                                            gCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTACTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGGGTGATGTTGCTGGAGCGGACCAGCCGCAAA‑GGGGAGTTACCCCAGGCGGGAATCTTGGTGTTGGGTCAGTCGGGTACCg                                                                                                                                                                                  >  1:256206/1‑228 (MQ=255)
                                               aTGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGTTAAGCTGGCCGTGATGTTGCTGGAGCGGACCACCCGTAAA‑GTGTTTTTCACCAAGGCGGGAATGCTCCTGGGGGGTCAGGCGCTAACCggtc                                                                                                                                                                               >  1:121824/1‑225 (MQ=255)
                                                tGGAGGATGGATAAGGAATATTCGTGATCTTGAGTACCGGGTGGCATGGGCTGAACCCCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGtgct                                                                                                                                                                              <  1:443496/228‑1 (MQ=255)
                                                 ggaggaTGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGGCTACCctgctg                                                                                                                                                                             >  1:390439/1‑228 (MQ=255)
                                                       gggatAATGAATGTTGTGGACCTTGAGTACCTGGTGGCATTGGCTAACCACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGa                                                                                                                                                                       <  1:513950/226‑1 (MQ=255)
                                                          atAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGGAAA‑GGGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGGGCTT‑CGGGAgggg                                                                                                                                                                   >  1:319361/1‑226 (MQ=255)
                                                            agtGAATATGCGGGATCGTGAGTACCGGGGGGCATTGGCTAAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGaa                                                                                                                                                                 <  1:59090/226‑1 (MQ=255)
                                                             aTGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGGATGCGGCTGGTGGATCAGGCGGGTCCCGGGCTG‑Cgggggggggaa                                                                                                                                                                >  2:436357/1‑218 (MQ=255)
                                                              tgaatttgggTGATGTTGAGTACCTGGTGCCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGaaa                                                                                                                                                                <  1:391970/218‑1 (MQ=255)
                                                              agAATTTTCGTTATCTTGAGTACTGGGGGGAATGGGCTGAACCCCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGaa                                                                                                                                                                 <  1:436357/225‑1 (MQ=255)
                                                               gAATATTCGTGATCTTGAGTCCCTGGTGGCAGTGGATGACTACCGTCATTTTCGGCGTGCGGCAGATTCCTGCCGAGTTAGCCAGCCGACGCGTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAg                                                                                                                                                               <  1:347282/227‑1 (MQ=255)
                                                                aaTATTCGTGATTTTGAGTACCTGGTGTCATGGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGGTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAg                                                                                                                                                               <  1:611712/226‑1 (MQ=255)
                                                                    ttCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATTTTGCTGGAGCGGACCAGCCGTAAA‑GGGTGGTTCACCCAGGGGGGAATGCTCCTGGTGGATCAGGCGCGTACCGGGCGG‑CGTGGGGGGAAAGACCTTa                                                                                                                                                         >  2:611712/1‑228 (MQ=255)
                                                                          aTCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGCTGTTCCTGGGGCGGACCAGCCGTAAA‑GGGTTGTTCACCCAGGCGGGAATGCTGCCGGTGAGTCAGGCGCGGACCCTGTCG‑GGGGAGGGGAAAGGCCTTAAAgaga                                                                                                                                                   >  2:310709/1‑228 (MQ=255)
                                                                                    ccTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTTCCTGGAGCGGACCAGCCGGAAA‑GTGTTGTTCACCCAGGCGGGAACGGTGCTGGGGGTGCAGGCGCGTACCGTGCTG‑CGGGGGGGGAAGAGCCTTAAagagcagggaaaccc                                                                                                                                         >  2:525658/1‑217 (MQ=255)
                                                                                     cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGGATGCTGCTGGTGGACCAGGCGCGTACCGTGCCG‑CGTGGGGTGACAGTCCTTCAAGAGATGGCAAGcccg                                                                                                                                        >  1:484983/1‑226 (MQ=255)
                                                                                     cTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGGGGGTCAGGGGCGGACCGTGCTGCCGTG‑GGGGAAAGTCCTTAAAGGGAGGGAAAGCcag                                                                                                                                        >  2:245131/1‑228 (MQ=255)
                                                                                       ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGGGAAAGGCCTTAAAGAGATGGCAAGCcagca                                                                                                                                      >  1:336317/1‑228 (MQ=255)
                                                                                       ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTAACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGCGGTGAAAGGCCTTAAAGAGAGGGCAAGccacca                                                                                                                                      >  1:200054/1‑228 (MQ=255)
                                                                                       ggtggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGCCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GGGTAGTTCACCCAGGCGGGAATGCTGCTGGGGGATCAGGCGCGTACCATGCTG‑CGAGAGGTGGAAGTCCTTAAAACGATGGAAAGCccgca                                                                                                                                      >  2:322926/1‑228 (MQ=255)
                                                                                          ggCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGGTGCCGGGGGGTCAGGGCCGTACCGTGCGG‑GCGGAGGGGAAAGCCCTTAAAAAAAGGGCAAGCCAGCAggg                                                                                                                                   >  2:250903/1‑228 (MQ=255)
                                                                                                ggcgaacccccGCCTTTTTGGGGGGGCGGCAGTTTCCTGCCAGGTTAGGCAGCCGGCGCTTGGGGGGCAAATTCGGAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGAc                                                                                                                             <  2:336317/219‑1 (MQ=255)
                                                                                                 gtgaaaacccGCCATTGTGGGGGTCGGGCAGATTCCTGCAACGTTGGCCGGCCGACGCTTAGGGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGAGGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGCGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACg                                                                                                                            <  2:258177/220‑1 (MQ=255)
                                                                                                          cccctttttGGGGGTGCGGCAGATCCCGGCCAGGTTAGCCGGGCGACGTTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGGGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGa                                                                                                                   <  2:95075/223‑1 (MQ=255)
                                                                                                          cGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGGTGTCCGGa                                                                                                                   >  1:534144/1‑228 (MQ=255)
                                                                                                                    ggCGTGCGGTAGCTTCCTGCTAGGTTTGCCGGCCGACGCTTAGGGGGCAAATTAGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGcaca                                                                                                         <  1:205467/228‑1 (MQ=255)
                                                                                                                         gCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCCTAAA‑GGGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCATGCGCGTACCTAGCTT‑CGTGGGGGAAAAGGCCTTCAAGAGATGGCAAGCCAGCAGGGCGGAACGAGGTCCCGACCGCTGCACCTTGGt                                                                                                    >  2:175398/1‑228 (MQ=255)
                                                                                                                          cGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGGGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCAACTTGGtt                                                                                                   >  1:12015/1‑228 (MQ=255)
                                                                                                                                  tCCTGGCAGGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTcc                                                                                            <  1:457974/227‑1 (MQ=255)
                                                                                                                                    ctgcagggTTAGCCAGGCGACGCTTGGCGGGCAAATTCGGAAGC‑TGGAAGGTGGGCTGGGCGTGATGTTGCTGGAGGGGGCCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCca                                                                                          <  1:322926/220‑1 (MQ=255)
                                                                                                                                      gccgtgttcggcaggcGGCGGTTAGCGGGCAAGTGCGTAAGC‑TGGAAGATGGGCTGGGCGGGGTGGTGGTGGTGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAg                                                                                       <  2:495669/217‑1 (MQ=255)
                                                                                                                                          cGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTg                                                                                    <  1:245131/227‑1 (MQ=255)
                                                                                                                                           tttAGCCAGCCGGAGTTTGGGGGGAAAATTCGGAAGC‑GGGAAGTAGGGCTGGGCGGGAGTGGGAGGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGa                                                                                  <  2:274445/227‑1 (MQ=255)
                                                                                                                                             tggccagccGAGGTTTGGCGGGAAAATTCGGAAGC‑TGGAAGATGGGCTGGGGGGGTTGTTGCTGGGGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGAcc                                                                                <  2:534144/226‑1 (MQ=255)
                                                                                                                                                   gccGACGGTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGGGCGGACCAGCCGTAAA‑GGGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTAcc                                                                           <  1:328450/227‑1 (MQ=255)
                                                                                                                                                          gttAGCGGGCAAATTGGTAAGC‑TGGAAGATGAGCGGGGGGGGAGGTTGCGGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCg                                                                   <  1:460736/227‑1 (MQ=255)
                                                                                                                                                              gCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCat                                                                >  1:8052/1‑227 (MQ=255)
                                                                                                                                                               ggggAAAATTCGTAGGC‑TGGAAGATGAGGGGGGTGTTGTTGTGCTGGGGGGGCCCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                              <  2:588696/227‑1 (MQ=255)
                                                                                                                                                                gggCAAATTCGTAAGG‑TGGAAGATTAGCTGGGCGTGATGTTGCTGGATCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                              <  1:160495/227‑1 (MQ=255)
                                                                                                                                                                gggCAAATTCGTAAGC‑TGGAAGATAAGCGGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCatat                                                              <  1:472698/227‑1 (MQ=255)
                                                                                                                                                                    aattttgTAAGG‑TGGAAGGTTGGGTGGGGGTGATTTTGGGGGGGGGGACCAGCCGTAAA‑GGGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGGGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                         <  2:97644/222‑1 (MQ=255)
                                                                                                                                                                    aaaTTGGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                         <  1:573263/228‑1 (MQ=255)
                                                                                                                                                                    aaaTTCGTAAGC‑TGAAAGATGGGCGGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCCCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATcc                                                         <  1:250903/228‑1 (MQ=255)
                                                                                                                                                                       tgcggaaGC‑TGGAAGGTGGGCGTGGCGGGATGTTGCTGGAGCGGACCTGCCGTAAA‑GTGTTGTTCGCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTa                                                      <  1:407996/223‑1 (MQ=255)
                                                                                                                                                                         tgTAAGG‑TGGAAGGTGGGGTGGGGGTTGTGTTGCGGGGGGGGACCGGCCGTAAA‑GTGTTGTTCCCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtg                                                    <  2:12015/227‑1 (MQ=255)
                                                                                                                                                                           tAAGC‑TGGAAGAGGAGCTGGGCGTGATGGTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtg                                                    <  1:407064/226‑1 (MQ=255)
                                                                                                                                                                           tAAGC‑TGGAAGACGAGCTGGGCGTGATGTTGGTGGAGGGGACGAGCCGTAAA‑GGGTTGTGCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTAtgc                                                   <  1:510971/227‑1 (MQ=255)
                                                                                                                                                                                   gAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGGGACGATGTCCGGACCGCTGCCCATTGGGTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAg                                           >  1:271417/1‑228 (MQ=255)
                                                                                                                                                                                   gAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCt                                                                          >  1:623500/1‑197 (MQ=255)
                                                                                                                                                                                           gCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTcc                                   >  1:463532/1‑228 (MQ=255)
                                                                                                                                                                                                   tggtgtgGCTGGCGCGGACCAGCCGTAAA‑GTGTTGTTCCCCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCt                             <  1:525658/223‑1 (MQ=255)
                                                                                                                                                                                                      tGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACg                                                                                                                            <  2:306891/128‑1 (MQ=255)
                                                                                                                                                                                                            tGGGGCGGCCCAGCCGTAAA‑GTTTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTAt                   <  1:562265/227‑1 (MQ=255)
                                                                                                                                                                                                            ggggagggACCCGCCGGTAA‑GTGTTGGTCGACCAGGCGGGAATGCGGCTGGTGGAACAGGCGCGATCCGGGCTG‑CGTGGGGTGAAAGTGCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCATACCTTTCCAAAGCTGGAAATGTAt                   <  1:310709/221‑1 (MQ=255)
                                                                                                                                                                                                                 cGGGCCAGCCGTAAA‑GTGTTTTTCACCCGGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTTCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCa              <  2:307769/227‑1 (MQ=255)
                                                                                                                                                                                                                  ggACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGGCGATGTCCGGACCGGTGCACATTGGTTTGATTCCCACAGTGGGAACGTACCTGCTCTCGCATATTATCCCTATTCTGAACCAGACCTTTCCCAAGCTGGAAATGAACCTGCAt             >  1:168679/1‑227 (MQ=255)
                                                                                                                                                                                                                       aGCCGTAAA‑GTGTGTTTCACCCAGGCGGGAATGCTGCTGGTTGTTCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGGCCTTAAAGAGATGGCAAGCCAGCATGGCGATACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGTCCTTTCCAAAGCTGGAAATGTATCTGCATGAAGca       <  1:481022/228‑1 (MQ=255)
                                                                                                                                                                                                                        gcgtaaaa‑GTGTTGGTCCCCCAGGCGGGAATTCTGCTGGTGGATGAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGcac      <  2:230743/223‑1 (MQ=255)
                                                                                                                                                                                                                         ccGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGcac      <  1:559386/227‑1 (MQ=255)
                                                                                                                                                                                                                             aaa‑GTGTTGTTCACCCAGGGGGGAATGCTGCTGGTGGATCAGGGCGGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAc  <  1:490018/227‑1 (MQ=255)
                                                                                                                                                                                                                             aaa‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAc  <  1:116499/227‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
ATAATCGTTCATTGCTATTCTACCTATCGCCATGAACTATCGTGGCGATGGAGGATGGATAATGAATATTCGTGATCTTGAGTACCTGGTGGCATTGGCTGAACACCGCCATTTTCGGCGTGCGGCAGATTCCTGCCACGTTAGCCAGCCGACGCTTAGCGGGCAAATTCGTAAGC‑TGGAAGATGAGCTGGGCGTGATGTTGCTGGAGCGGACCAGCCGTAAA‑GTGTTGTTCACCCAGGCGGGAATGCTGCTGGTGGATCAGGCGCGTACCGTGCTG‑CGTGAGGTGAAAGTCCTTAAAGAGATGGCAAGCCAGCAGGGCGAGACGATGTCCGGACCGCTGCACATTGGTTTGATTCCCACAGTTGGACCGTACCTGCTACCGCATATTATCCCTATGCTGCACCAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGAC  >  NC_000913/4158429‑4158875

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