Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,006,513 G→A 5.5% L84L (CTG→CTA)  ygeW → putative carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,006,5130GA 5.5% 239.8 / 2.9 73L84L (CTG→CTA) ygeWputative carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/26);  new base A (2/2);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

AAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCTATA  >  NC_000913/3006371‑3006661
                                                                                                                                              |                                                                                                                                                    
aaGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGc                                                                                                                                                <  2:447754/149‑1 (MQ=255)
      tGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGca                                                                                                                                          >  1:22626/1‑149 (MQ=255)
      tGGACGTTGCCGCAGCATTAAAAGAACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGca                                                                                                                                          >  2:395592/1‑149 (MQ=255)
        gACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                                                                                                                        >  4:152314/1‑149 (MQ=255)
         aCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaa                                                                                                                                       <  1:122833/149‑1 (MQ=255)
           gTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaaca                                                                                                                                     <  2:270500/149‑1 (MQ=255)
                cGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                                                                                >  3:345349/1‑149 (MQ=255)
                cGCAGCATGAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                                                                                                                >  1:280230/1‑149 (MQ=255)
                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtctc                                                                                                                              <  2:154883/149‑1 (MQ=255)
                  cagcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtctc                                                                                                                              <  2:22626/149‑1 (MQ=255)
                      aTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacg                                                                                                                          >  2:304608/1‑149 (MQ=255)
                       ttAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                                                                                                                         >  1:282366/1‑149 (MQ=255)
                                 tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATca                                                                                                               >  3:400025/1‑149 (MQ=255)
                                 tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTACTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATca                                                                                                               >  3:417840/1‑149 (MQ=255)
                                 tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTACTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATca                                                                                                               >  3:39984/1‑149 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                                                                        <  2:515356/133‑1 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                                                                                                                        >  1:515356/1‑133 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                                                                        >  1:226630/1‑149 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                                                                        <  4:345349/149‑1 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                                                                        <  2:280230/149‑1 (MQ=255)
                                        tGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCg                                                                                                        >  3:384224/1‑149 (MQ=255)
                                          aaaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGATATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCt                                                                                                      >  3:122660/1‑149 (MQ=255)
                                           aaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTc                                                                                                     >  3:300767/1‑149 (MQ=255)
                                            aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                                                                                    <  3:67343/149‑1 (MQ=255)
                                            aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                                                                                    >  2:206945/1‑149 (MQ=255)
                                            aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa                                                                                                    >  3:337217/1‑149 (MQ=255)
                                                tctcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGc                                                                                                >  2:33125/1‑149 (MQ=255)
                                                  tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGa                                                                                              >  1:292352/1‑149 (MQ=255)
                                                    aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGaaa                                                                                            >  1:583076/1‑149 (MQ=255)
                                                    aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGaaa                                                                                            >  1:223207/1‑149 (MQ=255)
                                                    aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGaaa                                                                                            >  2:323440/1‑149 (MQ=255)
                                                        aaaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgt                                                                                        <  3:83013/149‑1 (MQ=255)
                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgc                                                                                      <  2:583076/149‑1 (MQ=255)
                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgc                                                                                      <  1:304608/149‑1 (MQ=255)
                                                          aGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgc                                                                                      <  4:283184/149‑1 (MQ=255)
                                                            tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                                                                    >  2:397144/1‑149 (MQ=255)
                                                            tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                                                                    >  1:136205/1‑149 (MQ=255)
                                                            tCTTTAATAGTGGATAAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCGGCTTCTCTTATACTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACcgtgcgt                                                                                    >  4:250864/1‑149 (MQ=255)
                                                              tttAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGa                                                                                  >  4:75138/1‑149 (MQ=255)
                                                              tttAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGa                                                                                  >  1:345282/1‑149 (MQ=255)
                                                                 aaTAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                                                                      >  4:257783/1‑110 (MQ=255)
                                                                 aaTAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                                                                                                                      <  3:257783/110‑1 (MQ=255)
                                                                    aGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGc                                                                            <  2:283671/149‑1 (MQ=255)
                                                                     gTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGcc                                                                           <  3:177036/149‑1 (MQ=255)
                                                                      tGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCa                                                                          >  4:112441/1‑149 (MQ=255)
                                                                      tGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCACGACGAAGGCAAATCACAAATCGCTCACCGCGAAACCGTGCGTGAAACCGCCa                                                                          >  1:182154/1‑149 (MQ=255)
                                                                            aGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGa                                                                    <  1:543660/149‑1 (MQ=255)
                                                                                   tCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCtt                                                             >  2:531267/1‑149 (MQ=255)
                                                                                     cGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCt                                                           >  1:224976/1‑149 (MQ=255)
                                                                                       tATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTgc                                                         <  2:388935/149‑1 (MQ=255)
                                                                                            cgcgACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGa                                                    <  2:78742/149‑1 (MQ=255)
                                                                                                   aCTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTAtt                                             <  4:246434/149‑1 (MQ=255)
                                                                                                    cTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTg                                            <  4:300767/149‑1 (MQ=255)
                                                                                                     tCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTgg                                           >  4:266165/1‑149 (MQ=255)
                                                                                                        aCCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTAt                                        >  2:5969/1‑149 (MQ=255)
                                                                                                          ccGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACACCTCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTc                                      <  2:68575/149‑1 (MQ=255)
                                                                                                              aCCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGcga                                  >  3:391640/1‑149 (MQ=255)
                                                                                                                 cGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGcgacga                               <  4:334398/149‑1 (MQ=255)
                                                                                                                    ttCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGatat                            <  2:182154/149‑1 (MQ=255)
                                                                                                                         ttATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATc                       >  2:106151/1‑149 (MQ=255)
                                                                                                                         ttATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATc                       >  1:276532/1‑149 (MQ=255)
                                                                                                                            tGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTgg                    >  1:25088/1‑149 (MQ=255)
                                                                                                                              cTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGgc                  >  4:373828/1‑149 (MQ=255)
                                                                                                                               ttCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGgcg                 >  4:146292/1‑149 (MQ=255)
                                                                                                                                       ttAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAAc         >  4:144642/1‑149 (MQ=255)
                                                                                                                                         aaCCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGc       <  4:400025/149‑1 (MQ=255)
                                                                                                                                         aaCCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGc       >  3:12269/1‑149 (MQ=255)
                                                                                                                                         aaCCTACTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGc       <  4:39984/149‑1 (MQ=255)
                                                                                                                                         aaCCTACTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCAGGGCGCAGGCAACGc       <  4:417840/149‑1 (MQ=255)
                                                                                                                                           ccTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCt     >  2:13251/1‑149 (MQ=255)
                                                                                                                                            cTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCta    <  1:323440/149‑1 (MQ=255)
                                                                                                                                             tGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCtat   >  3:239668/1‑149 (MQ=255)
                                                                                                                                             tGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCtat   >  1:509476/1‑149 (MQ=255)
                                                                                                                                              gCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCtata  <  1:531267/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGCGAAACCGTGCGTGAAACCGCCAATATGATCTCCTTCTGCGCCGACGCTATTGGTATTCGCGACGATATGTATCTGGGCGCAGGCAACGCCTATA  >  NC_000913/3006371‑3006661

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑