Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,954 Δ1 bp 8.4% coding (123/567 nt) yecD → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,9540G. 8.4% 229.1 / 14.9 71coding (123/567 nt)yecDputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base . (4/2);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

TTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTA  >  NC_000913/1950812‑1951101
                                                                                                                                              |                                                                                                                                                    
ttGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAt                                                                                                                                               <  1:587706/150‑1 (MQ=255)
    tttttCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg                                                                                                                                           <  1:564565/150‑1 (MQ=255)
      tttCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGGGGGTGGTGATCGATTTCCAAGAAGGCATCTTACGTTTTGCCGGTGGTCCCCATACTGCCGATGAGGTAGTTAATCCCGCCGGGAACCTGGCGGCGAAATTTcg‑cc                                                                                                                                        >  2:89648/1‑149 (MQ=255)
           ggAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGc                                                                                                                                   <  2:1368556/150‑1 (MQ=255)
           ggAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                                                                   <  1:1079891/150‑1 (MQ=255)
                 aTTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCATACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAg                                                                                                                             >  1:1381807/1‑150 (MQ=255)
                      gCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                                                        >  1:986544/1‑150 (MQ=255)
                      gCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgt                                                                                                                        >  2:1164771/1‑150 (MQ=255)
                       cTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgtg                                                                                                                       >  2:560902/1‑150 (MQ=255)
                       cTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGct                                                                                                                                  <  2:746596/139‑2 (MQ=255)
                       cTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGct                                                                                                                                  >  1:746596/1‑138 (MQ=255)
                          gAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGttt                                                                                                                    <  1:390182/150‑1 (MQ=255)
                           aaCTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                                                                   >  1:504933/1‑150 (MQ=255)
                           aaCTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                                                                   <  1:175671/150‑1 (MQ=255)
                              ttAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTgg                                                                                                                <  2:10594/150‑1 (MQ=255)
                               tAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGt                                                                                                               >  1:1121616/1‑150 (MQ=255)
                                aaTGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCTGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTg                                                                                                              <  1:558832/150‑1 (MQ=255)
                                 aTGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                                                                             >  2:215301/1‑150 (MQ=255)
                                  tGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                                                                            >  2:346225/1‑150 (MQ=255)
                                  tGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                                                                            >  2:1082598/1‑150 (MQ=255)
                                     tAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGt                                                                                                         >  2:592669/1‑150 (MQ=255)
                                                    ggtggtGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                                                                          >  2:1162651/1‑150 (MQ=255)
                                                     gtggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                                                                         <  1:722145/150‑1 (MQ=255)
                                                       ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                                                       <  1:560902/150‑1 (MQ=255)
                                                              gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGaa                                                                                >  1:173366/1‑150 (MQ=255)
                                                               aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                                                                               <  1:1164771/150‑1 (MQ=255)
                                                                 ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCa                                                                             <  1:1378121/150‑1 (MQ=255)
                                                                   aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                                                           <  1:1129352/150‑1 (MQ=255)
                                                                    cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                                                          >  2:607486/1‑150 (MQ=255)
                                                                     aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                                                                         <  1:967048/150‑1 (MQ=255)
                                                                       gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaaa                                                                       >  2:467529/1‑150 (MQ=255)
                                                                        aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                                                                                 >  1:734321/1‑91 (MQ=255)
                                                                        aagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCgg                                                                                                                                 <  2:734321/91‑1 (MQ=255)
                                                                           gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGc                                                                   >  2:464308/1‑150 (MQ=255)
                                                                            cATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc                                                                  <  2:584790/150‑1 (MQ=255)
                                                                                    cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                                                           >  1:479106/1‑149 (MQ=255)
                                                                                    cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                                                           <  2:479106/149‑1 (MQ=255)
                                                                                              tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAcc                                                >  1:1078696/1‑150 (MQ=255)
                                                                                                gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGTTTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg                                              >  1:1238139/1‑150 (MQ=255)
                                                                                                 tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc                                             <  1:333160/150‑1 (MQ=255)
                                                                                                  ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCt                                            >  1:830422/1‑150 (MQ=255)
                                                                                                      aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg                                        >  2:1215158/1‑150 (MQ=255)
                                                                                                        aCTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAgtg                                      >  1:1047546/1‑150 (MQ=255)
                                                                                                          tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGtt                                    <  2:986544/150‑1 (MQ=255)
                                                                                                           gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTg                                   >  2:830818/1‑150 (MQ=255)
                                                                                                           gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTg                                   >  2:324781/1‑150 (MQ=255)
                                                                                                           cccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGCTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCAGGTCGATGCGCCCTAACCCGCTAAAGTGTTg                                   >  2:11095/2‑150 (MQ=255)
                                                                                                            ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGc                                  <  2:18527/150‑1 (MQ=255)
                                                                                                              gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGccc                                <  1:215301/150‑1 (MQ=255)
                                                                                                              gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGccc                                >  1:565606/1‑150 (MQ=255)
                                                                                                                   ggTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGGCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGaaaa                           <  2:1158715/150‑1 (MQ=255)
                                                                                                                    gTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAAt                          >  1:675098/1‑150 (MQ=255)
                                                                                                                       gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                       <  1:464308/150‑1 (MQ=255)
                                                                                                                         tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     <  1:1082598/150‑1 (MQ=255)
                                                                                                                         tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     >  1:798726/1‑150 (MQ=255)
                                                                                                                         tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     <  1:830818/150‑1 (MQ=255)
                                                                                                                         tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     >  2:1196445/1‑150 (MQ=255)
                                                                                                                          aaTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtgg                    >  2:954758/1‑150 (MQ=255)
                                                                                                                           aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGc                   >  2:220741/1‑150 (MQ=255)
                                                                                                                           aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGc                   >  1:983670/1‑150 (MQ=255)
                                                                                                                             cgcgCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCaa                 >  1:878287/1‑150 (MQ=255)
                                                                                                                              gcgcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     >  2:146097/1‑145 (MQ=255)
                                                                                                                              gcgcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtggtg                     <  1:146097/145‑1 (MQ=255)
                                                                                                                              gcgcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAAc                >  1:1335999/1‑150 (MQ=255)
                                                                                                                                      ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgct        >  1:702987/3‑150 (MQ=255)
                                                                                                                                      ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgct        >  1:576773/3‑150 (MQ=255)
                                                                                                                                       aGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgctg       <  2:1121616/150‑1 (MQ=255)
                                                                                                                                        gCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCtgctgc      <  2:173366/150‑1 (MQ=255)
                                                                                                                                        gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCtgctgc      >  1:1131530/1‑150 (MQ=255)
                                                                                                                                         cTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCa     <  2:1335999/150‑1 (MQ=255)
                                                                                                                                          tGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCAt    >  1:1158110/1‑150 (MQ=255)
                                                                                                                                            gcggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTa  <  1:607486/150‑1 (MQ=255)
                                                                                                                                            gcggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTa  >  1:1105500/1‑150 (MQ=255)
                                                                                                                                            gcggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTa  <  2:1381807/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TTGTTTTTTCAGGAGTCATTATGCTTGAACTTAATGCTAAAACCACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTA  >  NC_000913/1950812‑1951101

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