Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,085 A→C 6.4% intergenic (+29/+169) yphH → / ← glyA putative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,0850AC 6.4% 184.3 / 4.7 63intergenic (+29/+169)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base C (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATT  >  NC_000913/2683946‑2684232
                                                                                                                                           |                                                                                                                                                   
gaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCt                                                                                                                                            <  4:198790/149‑1 (MQ=255)
gaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGAGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCt                                                                                                                                            <  1:154513/149‑1 (MQ=255)
gaAGGACCGAGCGTGAAAGCGACGCAAATTCGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTGGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCt                                                                                                                                            <  5:106197/149‑1 (MQ=255)
        gAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATc                                                                                                                                    <  8:244543/149‑1 (MQ=255)
             tGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTAGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATaaa                                                                                                                               <  1:12493/149‑1 (MQ=255)
              gAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                              >  5:162159/1‑149 (MQ=255)
              gAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                              >  3:54442/1‑149 (MQ=255)
              gAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                              >  6:232558/1‑149 (MQ=255)
                     aCGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCTCAGTTACGACAGATTTGATCGCGCGATAACGTAGAAAGGCTTCCCGACGGAAGCCTTGATGATCATAAACGAAAAAt                                                                                                                       <  3:180978/149‑1 (MQ=255)
                        cAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGc                                                                                                                    <  2:31991/149‑1 (MQ=255)
                         aaaTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGcc                                                                                                                   >  7:114068/1‑149 (MQ=255)
                          aaTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCATGATGATCATAAACGAAAAATTGCCt                                                                                                                  <  7:189151/149‑1 (MQ=255)
                           aTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAATCGAAAAATTGCCt                                                                                                                  >  7:128858/1‑148 (MQ=255)
                           aTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGAGAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCt                                                                                                                  >  4:32774/1‑148 (MQ=255)
                              ggCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGAt                                                                                                               <  5:232558/148‑1 (MQ=255)
                               gCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgc                                                                                                             <  5:160850/149‑1 (MQ=255)
                                 tttGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATTGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGAt                                                                                                               >  5:157881/1‑145 (MQ=255)
                                 tttGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATATCGTAGAAAGGCTTGCCGAAGGATGCCTGGATGGTCATAAACGATATATTGCCTGATgcgg                                                                                                           >  4:143220/1‑148 (MQ=255)
                                               ccGCGCACAAGATGTGCTGGGAATTGGCTTTTTGTATGTTGAGCCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                             <  7:183920/149‑1 (MQ=255)
                                               ccGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCtttc                                                                                                                                         <  3:49770/105‑3 (MQ=255)
                                               ccGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCtttc                                                                                                                                         >  4:49770/1‑103 (MQ=255)
                                               ccGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                             >  7:12591/1‑149 (MQ=255)
                                                 gcgcACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                            >  3:78683/1‑148 (MQ=255)
                                                          aGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                                      >  8:99818/1‑145 (MQ=255)
                                                          aGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                                      >  7:68412/1‑145 (MQ=255)
                                                               cTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGc                                                                              >  7:83284/1‑148 (MQ=255)
                                                                      ttGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATAtt                                                                      >  4:204182/1‑149 (MQ=255)
                                                                         gCTAGTTGTATGTTGACGCGCAGTTACGACACATATGATGGCGCGAGAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATAAAAAACGAAAAATTGCCTGATGCTCTACGCTTATAAGGCCTACACGGAGATTGCAATATATTGaa                                                                   >  4:36963/1‑149 (MQ=255)
                                                                             tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAttt                                                                <  7:92073/148‑1 (MQ=255)
                                                                             tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAtt                                                                 <  8:72178/147‑1 (MQ=255)
                                                                             tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                                               <  7:99818/149‑1 (MQ=255)
                                                                             tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                                               >  8:16964/1‑149 (MQ=255)
                                                                              ttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                                               >  1:101052/1‑148 (MQ=255)
                                                                                gTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCa                                                             >  3:80386/1‑148 (MQ=255)
                                                                                gTATGTTGAGGAGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTAAGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCa                                                             >  4:13046/1‑148 (MQ=255)
                                                                                       gAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                      <  8:68412/148‑1 (MQ=255)
                                                                                       gAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                                      <  3:143220/148‑1 (MQ=255)
                                                                                       gAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                                     >  6:48939/1‑149 (MQ=255)
                                                                                        aGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                                     <  8:83284/148‑1 (MQ=255)
                                                                                            gcAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGc                                                <  5:48939/149‑1 (MQ=255)
                                                                                                   cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                         >  6:138954/1‑149 (MQ=255)
                                                                                                    gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                                              >  4:193308/1‑143 (MQ=255)
                                                                                                    gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                                              <  3:193308/143‑1 (MQ=255)
                                                                                                    gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                        <  5:138954/149‑1 (MQ=255)
                                                                                                    gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                        >  8:35771/1‑149 (MQ=255)
                                                                                                    gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                        >  1:73699/1‑149 (MQ=255)
                                                                                                     aCAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                       >  5:70284/1‑149 (MQ=255)
                                                                                                     aCAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                       >  7:77924/1‑149 (MQ=255)
                                                                                                      cAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                      >  5:78121/1‑149 (MQ=255)
                                                                                                            tGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                                 <  2:101052/148‑1 (MQ=255)
                                                                                                              aTGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgc                              >  4:23090/1‑149 (MQ=255)
                                                                                                               tGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgcc                             >  1:206809/1‑149 (MQ=255)
                                                                                                                         aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                    <  4:78683/148‑1 (MQ=255)
                                                                                                                          cGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                  >  2:110179/1‑149 (MQ=255)
                                                                                                                          agTAGAAAGGCTTCGCGTAGTAAGCCTTGATTATCATAAAGGAAAAATTGCCTGATGCTCTAGGCTTATCATGCCTACATGGAGATTGCAATATATTGAATTTGCAAATTTTTGTAGGCCTGATAAGGCGTTAACGCCGAATCCGGCa                   <  4:90143/147‑1 (MQ=21)
                                                                                                                              gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa               <  2:73699/148‑1 (MQ=255)
                                                                                                                               caaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca             >  4:204528/2‑149 (MQ=255)
                                                                                                                                    cTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg        <  2:206809/149‑1 (MQ=255)
                                                                                                                                     ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg        >  8:113860/1‑148 (MQ=255)
                                                                                                                                     gtCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGc       <  8:195517/148‑1 (MQ=255)
                                                                                                                                      tCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca      >  8:161538/1‑149 (MQ=255)
                                                                                                                                       cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca      >  3:28030/1‑148 (MQ=255)
                                                                                                                                        ccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca    <  3:32774/149‑1 (MQ=255)
                                                                                                                                          gaaggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAtt  >  3:110751/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATT  >  NC_000913/2683946‑2684232

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