Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,086 A→G 6.4% intergenic (+30/+168) 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,0860AG 6.4% 187.7 / 3.3 63intergenic (+30/+168)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (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 = 8.97e-01

GAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGA  >  NC_000913/2683946‑2684234
                                                                                                                                            |                                                                                                                                                    
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)
                                                                                                                                            aggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa  <  6:78121/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGA  >  NC_000913/2683946‑2684234

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: