Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,167 C→T 6.5% S473S (TCG→TCA yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,1670CT 6.5% 188.9 / 3.4 61S473S (TCG→TCAyoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (2/2);  total (32/29)
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.93e-01

GTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGT  >  NC_000913/1900021‑1900315
                                                                                                                                                  |                                                                                                                                                    
gttaACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCCCACGGGGAATGTGACCATTTGcc                                                                                                                                                      >  5:116040/1‑147 (MQ=255)
 ttaACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCgag                                                                                                                                                   <  8:19449/149‑1 (MQ=255)
   aaCGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGcc                                                                                                                                                      >  5:120443/1‑144 (MQ=255)
     cGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCa                                                                                                                                               >  7:385192/1‑149 (MQ=255)
              aTCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATAACAACGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGccct                                                                                                                                      >  2:367448/1‑148 (MQ=255)
                cTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCcgc                                                                                                                                     >  3:89208/1‑148 (MQ=255)
                       gTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTcg                                                                                                                             <  8:385192/149‑1 (MQ=255)
                        tGCGGCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTcgc                                                                                                                            <  1:367448/149‑1 (MQ=255)
                           gcCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGa                                                                                                                          >  8:310616/1‑148 (MQ=255)
                            ccGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGata                                                                                                                        <  7:315814/149‑1 (MQ=255)
                             cGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGatat                                                                                                                       >  6:55004/1‑149 (MQ=255)
                                ttGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATatcat                                                                                                                    <  4:89208/149‑1 (MQ=255)
                                  gCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATatcatc                                                                                                                   <  7:310616/148‑1 (MQ=255)
                                    tctTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgt                                                                                                                 >  2:171241/1‑148 (MQ=255)
                                    tctTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgt                                                                                                                 >  3:306508/1‑148 (MQ=255)
                                      tttAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATTATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCTGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcat                                                                                                              >  1:151287/1‑149 (MQ=255)
                                      tttAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCTGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcat                                                                                                              >  1:145949/1‑149 (MQ=255)
                                      tttAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcat                                                                                                              >  3:8043/1‑149 (MQ=255)
                                      tttAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcat                                                                                                              >  6:350838/1‑149 (MQ=255)
                                      tttAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTAAATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcat                                                                                                              >  5:145158/1‑149 (MQ=255)
                                        tAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATAtcatcgtcatcg                                                                                                            >  7:284258/1‑149 (MQ=255)
                                            aaTGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTCGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCaa                                                                                                        <  4:602/149‑1 (MQ=255)
                                             aTGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAg                                                                                                       <  5:350838/149‑1 (MQ=255)
                                              ttctaaCGAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGACCTACATCAATCACATCGCCCACACTGGGAATGAGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAgg                                                                                                      <  3:345972/145‑1 (MQ=255)
                                                      aGATCACCACGATATTCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACTCGGGGTATGTGACCATTTGCCGATATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAAc                                                                                              <  8:276354/149‑1 (MQ=255)
                                                      aGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAAc                                                                                              <  5:186665/149‑1 (MQ=255)
                                                            acacGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAg                                                                                        >  7:206989/1‑149 (MQ=255)
                                                                gATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCaa                                                                                         >  7:73888/1‑144 (MQ=255)
                                                                gATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCaa                                                                                         <  8:73888/144‑1 (MQ=255)
                                                                  tAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGAGACCATTTGCCGAGCTCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCctg                                                                                  >  1:266409/1‑149 (MQ=255)
                                                                       aTTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCtgctgc                                                                              >  7:61049/1‑148 (MQ=255)
                                                                           gCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGc                                                                         <  1:171241/149‑1 (MQ=255)
                                                                           gCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGc                                                                         <  8:61049/149‑1 (MQ=255)
                                                                               cAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCAt                                                                       <  6:145158/147‑1 (MQ=255)
                                                                                aaTGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACCCGGGGAATGTGACCATTTGCCGCGATCACGAGGCCCGCGACCGTCGCGATATCATCTTCATCGGCAAGGTGCTCACCATCAAGCGCCTGCTGCACGGCACGc                                                                     >  2:295928/1‑148 (MQ=255)
                                                                                    atgGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAAt                                                                 <  4:8043/148‑1 (MQ=255)
                                                                                    atgGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATc                                                                >  3:340423/1‑149 (MQ=255)
                                                                                     tgGTGATATGCAGTGGACCTACTTCAATCACATAGCCAACACGGGGAATGTGACCATTTTCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCGTGCTGCAAGGCATGCAAATCt                                                               <  4:340423/149‑1 (MQ=255)
                                                                                        tGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTgca                                                            >  4:34116/1‑147 (MQ=255)
                                                                                        tGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTa                                                            <  5:164293/149‑1 (MQ=255)
                                                                                            atGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTAc                                                           >  8:233263/1‑146 (MQ=255)
                                                                                            atGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTAc                                                           <  7:233263/146‑1 (MQ=255)
                                                                                             tGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCTGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGc                                                        <  2:145949/148‑1 (MQ=255)
                                                                                             tGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCTGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGc                                                        <  2:151287/148‑1 (MQ=255)
                                                                                             tGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCAGCTGCAAGGCATGCAAATCTGTACCGcc                                                       >  3:139205/1‑149 (MQ=255)
                                                                                              gCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTg                                                              <  3:115100/141‑1 (MQ=255)
                                                                                              gCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTg                                                              >  4:115100/1‑141 (MQ=255)
                                                                                              gCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCt                                                      >  5:124534/1‑149 (MQ=255)
                                                                                                    gCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTAc                                                 >  7:341515/1‑148 (MQ=255)
                                                                                                     cccTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTAcca                                               <  5:64912/149‑1 (MQ=255)
                                                                                                        tACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTAccagcc                                            >  6:190943/1‑149 (MQ=255)
                                                                                                         aCATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCaa                                                                                                        >  6:25252/1‑88 (MQ=255)
                                                                                                         aCATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCaa                                                                                                        <  5:25252/88‑1 (MQ=255)
                                                                                                           atcaatcaCATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACcagccagc                                         <  6:120443/149‑1 (MQ=255)
                                                                                                                tcaCATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGCCCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCa                                    >  7:305285/1‑149 (MQ=255)
                                                                                                                tcaCATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCa                                    <  5:190943/149‑1 (MQ=255)
                                                                                                                  acaTCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGACTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCAc                                   >  5:9386/1‑148 (MQ=255)
                                                                                                                    aTCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCAt                                <  8:284258/149‑1 (MQ=255)
                                                                                                                          cacaCGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTaa                          <  5:55004/149‑1 (MQ=255)
                                                                                                                           acacGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAAt                         <  8:206989/149‑1 (MQ=255)
                                                                                                                                        gACCATTTGCCGAGCTCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGtt            >  5:98299/1‑149 (MQ=255)
                                                                                                                                        gACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATCATTTACGGCgat            >  6:226687/1‑147 (MQ=255)
                                                                                                                                            aTTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGTTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGt        >  3:107884/1‑149 (MQ=255)
                                                                                                                                                 ccGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGa                       >  2:277868/1‑129 (MQ=255)
                                                                                                                                                 ccGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGa                       <  1:277868/129‑1 (MQ=255)
                                                                                                                                                  cGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGt  >  4:343336/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTAACGCTTACTCATCTTCATCGTGCGCCGGTTGCTCTTTAACAATGCGAACCAGATCAACACGATAATCATTGGCTTCAATGATGGTGATATGCAGTGGCCCTACATCAATCACATCGCCCACACGGGGAATGTGACCATTTGCCGAGATCACGAGGCCCGCGACCGTCGCGATATCATCGTCATCGGCAAGGTGCTCAACATCAAGCGCCTGCTGCAAGGCATGCAAATCTGTACCGCCTTTTACCAGCCAGCCGTCACCATCAGTAATGATTTCCGGCGTTTCGTCAGCGT  >  NC_000913/1900021‑1900315

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