Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,336,411 2 bp→AT coding (83‑84/1260 nt) murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,4110GA97.4% 109.8 / NA 40I28N (ATC→AATmurAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (13/25);  total (14/26)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*NC_0009133,336,4120AT94.7% 113.6 / ‑1.9 40I28N (ATC→AAT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (13/23);  total (14/26)
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.30e-01

CCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATAT  >  NC_000913/3336176‑3336639
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                     
ccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAgcagca                                                                                                                                                                                                                        <  1:197046/251‑1 (MQ=255)
 cTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCTCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCAt                                                                                                                                                                                                                       >  2:185264/1‑251 (MQ=255)
                           ccAGCGCCTAGATAGCAGCAAGCATGGTGTGAACCATATCGTATGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACATAACCAGTACGTTCTACTTTCTCGCCCAGCTGGCTTAGCAGCTTTATTGATGTATCGACGTCTTTCATTTTCGGGACTTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGGGCGTCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGTAAATTGTGCCTTCg                                                                                                                                                                                             <  2:80988/251‑1 (MQ=255)
                                  ccAGATAGAAGCCCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATTTATCGACGTCTTTCAGTTGCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGa                                                                                                                                                                                      <  1:320242/251‑1 (MQ=255)
                                                         tAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACg                                                                                                                                                                <  1:130734/250‑1 (MQ=255)
                                                          aaCCAGATCGAAAGGTGCTCAGAGTACATTAACGTCGAGGGCATCACTATGCACAGAACCATTACGTTCTACTATCGCACCCCGCTGGCTTAGCAGCTTCCTTGAGGTATCGACGTCTTTCAGTTTAGGGACGTTATGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCACAAAGGATCGGCATAGCAGCATTTTTAGAGCAGGAAATTATGCATTCGCACTAGAGCTTCGTTGCCCCATGAACACGa                                                                                                                                                               >  2:18002/1‑250 (MQ=255)
                                                              aGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAAttt                                                                                                                                                          >  1:135604/1‑251 (MQ=255)
                                                                  cGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATcc                                                                                                                                                      <  1:183884/251‑1 (MQ=255)
                                                                     aaGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAtt                                                                                                                                                   >  1:150240/1‑251 (MQ=255)
                                                                        gTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAttt                                                                                                                                                  >  1:242571/1‑249 (MQ=255)
                                                                          gcgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCa                                                                                                                                                     <  1:104980/244‑1 (MQ=255)
                                                                          gcgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCa                                                                                                                                                     >  2:104980/1‑244 (MQ=255)
                                                                           cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGttt                                                                                                                                             <  1:62147/251‑1 (MQ=255)
                                                                             cAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGt                                                                                                                                           <  1:370419/251‑1 (MQ=255)
                                                                                aaTACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtct                                                                                                                                       <  1:172844/251‑1 (MQ=255)
                                                                                        accgtcggtggCATCAATACGCCCAGAGCCATTACGTTCTTCTTTCGCACCAAGCGGGCTTAGCAGCTTCATTAAGGTATCGACGTCTTGCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCTCTGGAGCTTCGTTGGCACCTGTACACGAGATTTATCCATTTAGTTTGGTCTCAGTTaa                                                                                                                                <  2:79729/241‑1 (MQ=255)
                                                                                        aaCGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAAc                                                                                                                               >  2:371689/1‑251 (MQ=255)
                                                                                          cGTCGCGGGCATCAATATGCACAGAACCCTTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCtt                                                                                                                                                                                  >  2:176377/1‑199 (MQ=255)
                                                                                          cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCtt                                                                                                                                                                                  <  1:176377/199‑1 (MQ=255)
                                                                                          cGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACGTTCGCACCCAGCTGGCTTAGCAGCTTCATTTATGTATCGTCGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATGAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACaa                                                                                                                             <  1:208278/251‑1 (MQ=255)
                                                                                              gcgGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCTCCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCa                                                                                                                         <  1:185264/251‑1 (MQ=255)
                                                                                               cgGGCATCAATATGCACAGAACCATTACGTTCTACTTGCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCat                                                                                                                        <  1:269519/251‑1 (MQ=255)
                                                                                                   cataaaTATGCACAGAACCATTCTGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGGTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGATCTTCGTTGGCCCCTGTACACGAAATTTATCCATTTAGTTTTTTCTCAGTTAACAACTCATATcc                                                                                                                    <  2:337985/247‑1 (MQ=255)
                                                                                                       aaTATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGTCGTCTTTCAGTTTCGGGACGTTCTGGATCTCTCCCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCata                                                                                                                       <  1:257568/244‑1 (MQ=255)
                                                                                                       aaTATGCAAAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGCCGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTACGACAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCata                                                                                                                       >  2:257568/1‑244 (MQ=255)
                                                                                                        atatgccaaGAACCATTACGTTCTACTTTCCCCCCCAGCTGGCTTAGCATCTTCATTGATGTATCGACGGCTTCCAGTTTCGGGACGTTCTGGTTCTCTACCGGTTCGTCCGTCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCAGGAAATTGTTACTTCGCCCTGGAGCTTCTTAGGCCCCTGTACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATAGCCGCTAc                                                                                                               <  2:63600/243‑1 (MQ=255)
                                                                                                                         ccgtactcccttCGCCCCCAGCTGTCGTAGCAGCTTTATTGATTTTTCGTCGTCTTTCCTTTGCGGGCCGTTCTGGATCTCTACCGGTCCTTCCGCCAGTAGTGCGGCAAAAAGAGTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTTACTTCGCCCTGGAGCTTCGTTGGCCCCTTAACGCGTAACTTATCCATTTAGTTTGTTCTCAGTTAACAAATCATATCCGCTACCGGCgcagcgctcata                                                                                               <  2:329836/240‑12 (MQ=255)
                                                                                                                           gTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCt                                                                                            <  1:198369/251‑1 (MQ=255)
                                                                                                                              cTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCa                                                                                          >  1:148425/1‑250 (MQ=255)
                                                                                                                                aCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTTTTCTCAGTTAACTATTCATATCCGCTACCGGCGAATCGCCCATAGCTCaaaa                                                                                       <  1:416633/251‑1 (MQ=255)
                                                                                                                                   ttctctcccAGGTGGCGTAGCCGCTTAATTGATGTATCGATGTCTTTCAGTTTCGGGACGTTCTGTATCGCTACCGGTTTTTCCGCCAGTAGTGCGGCAAAAAGACTAGGCAGAGCTGCAGTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTATTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCGCAAAAg                                                                                      <  2:127887/243‑1 (MQ=255)
                                                                                                                                               gggCTTAGTAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCg                                                                        <  1:318788/250‑1 (MQ=255)
                                                                                                                                                 gCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGAt                                                                      >  2:4173/1‑251 (MQ=255)
                                                                                                                                                            tCATTGATGTATCGACGTCTTTCAGTTGCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCa                                                            <  1:316312/250‑1 (MQ=255)
                                                                                                                                                                tGATGTATCGATGTCTTCCAGTTCCGGGTCTATCTTGATCTCTCCCGGTTCTTCCGCCAGGAGTGCGGCAAAAAGATTAGGCAGAGCTGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGATCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTATTTGTTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGGTCAGTTTGCGATCGCGCGCCCACTcc                                                        <  2:285831/250‑1 (MQ=255)
                                                                                                                                                                     tATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCCggg                                                  >  1:76839/1‑251 (MQ=255)
                                                                                                                                                                           cGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGCACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCAGCTACCGGCGACTCGCCCCTAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGCATAc                                            >  2:185830/1‑251 (MQ=255)
                                                                                                                                                                                           gACGTTCTGGATCTCTACCGGTTCTGCCGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGACCACTCCGCAGGGGTATACGCTTTGATCGACACAg                            <  1:182959/251‑1 (MQ=255)
                                                                                                                                                                                                      tctttACCGGTTCTTCCTCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGTCTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTTCGCAGGGGGATACGCTTTGATCGACACAGCATGAATGCg                  <  2:264733/246‑1 (MQ=255)
                                                                                                                                                                                                                     ccGCCAGTAGTGCGGCAAAAAGATTAGGCAGAGCAGCATTTTTAGCGCCGGCAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAatatat  >  1:192027/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                     
CCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCATATCCGCTACCGGCGAATCGCCCATAGCTCAAAAGCCGTTCAGTTTGCGATCGCGCGCCCACTCCGCAGGGGTATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATAT  >  NC_000913/3336176‑3336639

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: