Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 244,550 T→A 6.8% V75E (GTG→GAG)  yafJ → type 2 glutamine amidotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913244,5500TA 6.8% 209.6 / 4.2 70V75E (GTG→GAG) yafJtype 2 glutamine amidotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base A (3/2);  total (36/34)
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.84e-01

TGGAACCGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCG  >  NC_000913/244407‑244698
                                                                                                                                               |                                                                                                                                                    
tggAACCGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCt                                                                                                                                                  >  3:224141/1‑148 (MQ=255)
 ggAACCGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGACGGTAAAGGCTGTCGCACATTTAAAGATCCCCCACCCAGCTTGAATTCCCCCATCGCCAACCTTGTCCCGGACTACCCGATAAAATAGTGTTCGGTGCTGGCGCa                                                                                                                                                >  6:256074/1‑149 (MQ=255)
     ccGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGCTCCACAACCCAGCTTTAATTCCCCCATCGCCAACCTTGTCCAGGACTACCCGATAAAATCCGGTTCGGTGGTGGCTCATAtt                                                                                                                                            >  2:121907/1‑149 (MQ=255)
      cGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTc                                                                                                                                           <  1:160599/149‑1 (MQ=255)
      cGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTc                                                                                                                                           <  5:240025/149‑1 (MQ=255)
         gCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGc                                                                                                                                         >  4:255278/1‑148 (MQ=255)
          ccACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGc                                                                                                                                         >  4:101667/1‑147 (MQ=255)
            acaTAAAGATGGCTGGGGCATTACCTTTTACGACGGTTAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTACAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAgg                                                                                                                                     <  6:20010/149‑1 (MQ=255)
                 aaGATGGCTGGGGCATTACCTTTTAAGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAAt                                                                                                                                <  2:76773/149‑1 (MQ=255)
                  tgATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATc                                                                                                                               >  5:231557/2‑149 (MQ=255)
                   gATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCg                                                                                                                              <  3:311652/149‑1 (MQ=255)
                       gCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCgggg                                                                                                                           <  2:129667/148‑1 (MQ=255)
                             gCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTg                                                                                                                    <  3:101667/149‑1 (MQ=255)
                                ttACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGgcg                                                                                                                 >  1:209737/1‑149 (MQ=255)
                                ttACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGACGAGGTGgcg                                                                                                                 >  8:193638/1‑149 (MQ=255)
                                  aCCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGgcgc                                                                                                                <  8:93394/148‑1 (MQ=255)
                                   ccTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGGGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCt                                                                                                               <  4:293525/148‑1 (MQ=255)
                                    cTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTgg                                                                                                             <  5:259514/149‑1 (MQ=255)
                                     ttttACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTgg                                                                                                             >  3:180862/1‑148 (MQ=255)
                                        ctaaaacgTAAAGGCGGTCGCACAATTAAACATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACATGTCCCGGCCTACCCGATAAAATCCTGTTCGGTGGGGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGaaaa                                                                                                         <  3:166069/142‑1 (MQ=255)
                                          cGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACACAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGTTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATa                                                                                                       <  3:187163/149‑1 (MQ=255)
                                           gAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGCTGGAGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATAc                                                                                                      >  7:246644/1‑149 (MQ=255)
                                               gTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTTCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCAc                                                                                                  <  2:251842/149‑1 (MQ=255)
                                                tAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTACTTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCAcc                                                                                                 <  6:311055/149‑1 (MQ=255)
                                                 aaaGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTTCCCCATCGCCAAACTTGTACAGGACTACCAGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTTGCGCTGGAAAATACTCaac                                                                                                 >  5:206269/1‑146 (MQ=255)
                                                 aaaGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCAcc                                                                                                 >  1:267443/1‑148 (MQ=255)
                                                  aaGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCAccc                                                                                                <  5:196640/148‑1 (MQ=255)
                                                  aaGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCa                                                                                               <  5:241286/149‑1 (MQ=255)
                                                            cacaTTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAGCTTGTCCAGGACTACCCGATAAAATCCTGATCGGTGGTGCCTCATATTCGCCAGGCAAATCGGGGCGAGCTGGCGCTGGAAAATACTCACCCATTAAccccc                                                                                      >  6:257871/1‑146 (MQ=255)
                                                             acaTTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGa                                                                                    <  4:224141/149‑1 (MQ=255)
                                                                 ttAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGCGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTa                                                                                >  8:256353/1‑149 (MQ=255)
                                                                   aaaGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTGCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATg                                                                              >  4:203540/1‑149 (MQ=255)
                                                                    aaGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATg                                                                              <  3:230540/148‑1 (MQ=255)
                                                                    aaGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATg                                                                              <  2:209737/148‑1 (MQ=255)
                                                                    aaGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTAGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATTCTCACCCATTTACCCGCGAGTTATgg                                                                             >  1:299376/1‑149 (MQ=255)
                                                                     aGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATggg                                                                            >  6:146766/1‑149 (MQ=255)
                                                                      gATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           <  3:255278/149‑1 (MQ=255)
                                                                      gATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           <  5:155841/149‑1 (MQ=255)
                                                                       aTCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCAtt                                                                                             >  4:62900/1‑130 (MQ=255)
                                                                       aTCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCAtt                                                                                             <  3:62900/130‑1 (MQ=255)
                                                                          cacaACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCggg                                                                                                                            <  7:141797/96‑1 (MQ=255)
                                                                          cacaACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCggg                                                                                                                            >  8:141797/1‑96 (MQ=255)
                                                                           acaACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGt                                                                       >  4:131416/1‑148 (MQ=255)
                                                                                  aGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGAc                                                               <  6:313261/149‑1 (MQ=255)
                                                                                       tAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTTGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTAt                                                           <  6:322427/148‑1 (MQ=255)
                                                                                         aTTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGcc                                                        >  3:90604/1‑149 (MQ=255)
                                                                                            cccccATCGCCAACCTTGTCCAGGACTACCCGATAAAAACGTGTTGGGTGGTGGGTCTTATTCGCCGGGCTAATCGGGGCGGGGTGGCGCGGGAAAATACTCACCAATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCa                                                      >  3:280948/1‑148 (MQ=255)
                                                                                            cccccATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATATACCCGCGAGTTATGGGGGCGTAACTGCACTTATGCCCa                                                      >  6:306746/1‑148 (MQ=255)
                                                                                                    gCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACgggtg                                             >  7:141681/1‑146 (MQ=255)
                                                                                                     ccAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACaa                                            >  5:105719/1‑149 (MQ=255)
                                                                                                      cAAACTTGTCCAGGACTACCCGATAAACTCCTGTTCGGTGGTGGCGCATATTCGCCAGGCTCATCGGGGCGAGGTGGCGCTGGAAAATACTCCCCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAAc                                           <  1:121907/149‑1 (MQ=255)
                                                                                                       aaaCTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGc                                                         <  4:288745/134‑1 (MQ=255)
                                                                                                       aaaCTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGCGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGc                                                         >  3:288745/1‑134 (MQ=255)
                                                                                                              tCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACggg                                    <  2:267443/148‑1 (MQ=255)
                                                                                                                cAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTa                                 >  8:36345/1‑149 (MQ=255)
                                                                                                                 aGGACTACCCGATAAAATCCTGTTCGGTGGAGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           >  4:163445/1‑106 (MQ=255)
                                                                                                                 aGGACTACCCGATAAAATCCTGTTCGGTGGAGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           >  4:158259/1‑106 (MQ=255)
                                                                                                                 aGGACTACCCGATAAAATCCTGTTCGGTGGAGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           <  3:163445/106‑1 (MQ=255)
                                                                                                                 aGGACTACCCGATAAAATCCTGTTCGGTGGAGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATgggg                                                                           <  3:158259/106‑1 (MQ=255)
                                                                                                                   gACTAACCGATAAAAACCTGATAGGTTCTCGCTCATATTCGAAAGCCTAATCTGGGCGAGGTGGCGCTAGAAAATACTAACCCATTTACCCGCGAGTTATAGATGCGTAACTTGACTTATGCCCATAACGTACAACTGACGGGCAACa                               >  4:203689/1‑148 (MQ=255)
                                                                                                                      tACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCAAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTAATGCGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAAt                            >  8:305206/1‑148 (MQ=255)
                                                                                                                        cccGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCAc                         >  4:216612/1‑149 (MQ=255)
                                                                                                                        cccGATAAAATCCTGTTCGGTGGGGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCa                          <  6:231557/148‑1 (MQ=255)
                                                                                                                         ccGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACaa                                            <  2:132784/129‑1 (MQ=255)
                                                                                                                         ccGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACaa                                            >  1:132784/1‑129 (MQ=255)
                                                                                                                         ccGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACt                        <  3:131416/149‑1 (MQ=255)
                                                                                                                         ccGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACt                        >  1:195496/1‑149 (MQ=255)
                                                                                                                          cGATAAAATCCTGTTCGGTGGTGACTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTg                       >  8:215958/1‑149 (MQ=255)
                                                                                                                                     tGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCa             >  1:321721/1‑148 (MQ=255)
                                                                                                                                       ttCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAAc           <  2:321721/148‑1 (MQ=255)
                                                                                                                                          ggtggtGGCTCATATTCGCCAGGCTAATCGGGACGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTcc       <  7:193638/149‑1 (MQ=255)
                                                                                                                                          ggtggtGACTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTcc       <  7:215958/149‑1 (MQ=255)
                                                                                                                                               tggCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCg  <  4:280948/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGGAACCGGGCCACATAAAGATGGCTGGGGCATTACCTTTTACGAAGGTAAAGGCTGTCGCACATTTAAAGATCCACAACCCAGCTTTAATTCCCCCATCGCCAAACTTGTCCAGGACTACCCGATAAAATCCTGTTCGGTGGTGGCTCATATTCGCCAGGCTAATCGGGGCGAGGTGGCGCTGGAAAATACTCACCCATTTACCCGCGAGTTATGGGGGCGTAACTGGACTTATGCCCATAACGGACAACTGACGGGCTACAAATCACTGGAAACCGGCAACTTCCGCCCG  >  NC_000913/244407‑244698

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