Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,999 T→G 10.5% L217V (TTG→GTG)  yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,9990TG10.5% 958.8 / 27.2 345L217V (TTG→GTG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (117/184);  new base G (42/2);  total (159/186)
Fisher's exact test for biased strand distribution p-value = 9.47e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

ATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACG  >  NC_000913/71850‑72148
                                                                                                                                                     |                                                                                                                                                     
 tatTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGttt                                                                                                                                                      <  1:4931616/150‑1 (MQ=255)
     ccTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTttgttg                                                                                                                                                  >  2:3874294/1‑150 (MQ=255)
       tGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTtgtg                                                                                                                                                >  2:3362418/1‑150 (MQ=255)
       tGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTtgtg                                                                                                                                                <  2:3624566/150‑1 (MQ=255)
       tGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTtgtg                                                                                                                                                >  1:4219288/1‑150 (MQ=255)
       tGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTtgtg                                                                                                                                                >  2:1352222/1‑150 (MQ=255)
       tGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTtgtg                                                                                                                                                >  1:2223532/1‑150 (MQ=255)
          cGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCt                                                                                                                                             <  2:1122450/150‑1 (MQ=255)
          cGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCt                                                                                                                                             <  2:2847843/150‑1 (MQ=255)
          cGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCt                                                                                                                                             <  1:738923/150‑1 (MQ=255)
            gAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGGTGTGGCGGa                                                                                                                                           >  1:1932748/1‑150 (MQ=255)
               tGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGc                                                                                                                                        <  2:2507433/150‑1 (MQ=255)
               tGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGc                                                                                                                                        <  2:601958/150‑1 (MQ=255)
                   gagCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGtt                                                                                                                                    <  2:2144973/150‑1 (MQ=255)
                        gTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGAttt                                                                                                                               >  1:4955223/1‑150 (MQ=255)
                          gAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCt                                                                                                                             <  1:4795775/150‑1 (MQ=255)
                           aGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTg                                                                                                                            <  2:2956067/150‑1 (MQ=255)
                            gTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGc                                                                                                                           <  1:3957999/150‑1 (MQ=255)
                                  aTGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGGTGTGGCTGACGCCGTTGATTTCTGCCATCgg                                                                                                                     >  1:4138391/1‑150 (MQ=255)
                                     gTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgt                                                                                                                  >  2:400732/1‑150 (MQ=255)
                                       tgctgcTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATGATTTGTCCCGCGGTCGTGTGTTGTGGCGGTCGCCGTTGGTTTATGCCCTCGGCgtgg                                                                                                                >  1:2663352/1‑150 (MQ=255)
                                        gctgctGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggt                                                                                                               <  2:564833/150‑1 (MQ=255)
                                         ctgctgGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCgtggtg                                                                                                              <  2:3778033/150‑1 (MQ=255)
                                          tgctgGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGAAGCCGTTTATTTCTGCCATCGGCGtggtgg                                                                                                             >  2:1251607/1‑150 (MQ=255)
                                            gtgGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGTGTGATTATTTGTGCCGCGGTCGTTTGTTTTGGCTTTCGGCGTTGTTTTATTGCATGGGGGTGGTTgcg                                                                                                           >  2:956858/2‑150 (MQ=255)
                                             tgGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcgc                                                                                                          >  2:2942326/1‑150 (MQ=255)
                                             tgGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcgc                                                                                                          <  1:2128293/150‑1 (MQ=255)
                                              ggCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCt                                                                                                         >  2:495944/1‑150 (MQ=255)
                                              ggCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCt                                                                                                         <  1:532864/150‑1 (MQ=255)
                                                  aCAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtg                                                                                                     <  1:839492/150‑1 (MQ=255)
                                                  aCAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtg                                                                                                     <  1:3265197/150‑1 (MQ=255)
                                                  aCAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtg                                                                                                     <  2:4219288/150‑1 (MQ=255)
                                                   cAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtgg                                                                                                    <  1:5700468/150‑1 (MQ=255)
                                                   cAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggtgg                                                                                                    >  2:4128601/1‑150 (MQ=255)
                                                    aGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                                                                   <  2:2056952/150‑1 (MQ=255)
                                                    aGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                                                                   <  2:2399453/150‑1 (MQ=255)
                                                    aGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                                                                   <  2:1967924/150‑1 (MQ=255)
                                                      cGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGgtgt                                                                                                 <  2:3006532/150‑1 (MQ=255)
                                                        gtgtTTTAGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCTACTCACCGTTTGAGTCATTATTTGTCCCGCGGTAGTTTTTTTTGGCTGAAACCCTTGAAATCTGCCAAAGGCCAGGTTGCGCAGGTGGTGTTa                                                                                               >  1:4649458/1‑150 (MQ=255)
                                                            ttttGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTc                                                                                           <  1:3196359/150‑1 (MQ=255)
                                                             tttGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTTCCCTAGGCGTGGTGGCGCTGGTGGTGTTAATTCg                                                                                          >  1:4184917/1‑150 (MQ=255)
                                                                 tgGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCAc                                                                                      <  2:5475380/150‑1 (MQ=255)
                                                                   ggTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTTCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCAccc                                                                                    >  2:1023412/1‑150 (MQ=255)
                                                                   ggTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGGGTGAGTCCTTATTTGTCCCGCGGTCGTGTGTTGTGGGTGGCGCCGGTGTTTTCTTCCATCGGCGTGGGGGGGCTGGGGGTGTTAATTCGcccccc                                                                                    >  1:1124241/1‑150 (MQ=255)
                                                                    gTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCAAAGGCGTGGTGGCGCAGGTGGGGTTAATTCGCCAcccc                                                                                   >  1:912410/1‑149 (MQ=255)
                                                                      tggtggCTGTCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGtt                                                                                 <  1:1469711/150‑1 (MQ=255)
                                                                            ctggctgTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCTGACGCCGGTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTAGCCACCCCTTGATGcc                                                                           >  2:1532954/1‑150 (MQ=255)
                                                                             tggctgTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCGGACGCCGTTGATTTGATCCCTCGGCGTGGGGGCCCAGGTGGGGGTAATTCGCCCCCCCGTGATGCCg                                                                          >  1:2766725/1‑150 (MQ=255)
                                                                                     cTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTatat                                                                  <  1:554594/150‑1 (MQ=255)
                                                                                      tGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTTGCTGACGCCGGTGGTTTCTGCCCACGGCGTGGTGGGGCGGGTGGTGGTAATTCGCCCCCCCGTGATGCCCGGGTATATc                                                                 >  1:3439472/1‑150 (MQ=255)
                                                                                      tGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCTGACGCCGTTGATTTCTGCCCTCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCCTTGATGCCGGTGTATATc                                                                 >  2:4058983/1‑150 (MQ=255)
                                                                                        gCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGa                                                               <  2:2613646/150‑1 (MQ=255)
                                                                                          ggTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCTGTCGCCGGTGATTTCTGCCAACGGCGGGGTGGCGCTGGTGGGGTTATTTAGCCACCCGTTGATGCCGGGGTATATCGata                                                             >  1:666607/1‑150 (MQ=255)
                                                                                           gTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGatat                                                            >  2:5202865/1‑150 (MQ=255)
                                                                                            ttATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAtt                                                           >  2:4791620/1‑150 (MQ=255)
                                                                                             tATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGGGTTGTGGGGGACGCCGGTGCTTTGTGGCAACGGCGGGGTGGCGCTGGTGGGGTTGATTCGCCACCCCGTGATTGCGGGGTATATAGATAttt                                                          >  1:1131769/1‑150 (MQ=255)
                                                                                             tATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCAATATTTGTCCCGCGGTCGTTTGTTGTTGCTGACGCCGTTGATTTATGCCAACGGCGTGGTGGCGCTGGTGGTGTTAATTAGCCACCCCTTGTTGGCGGTGTATATAGAGAttt                                                          >  1:3736638/1‑150 (MQ=255)
                                                                                               tGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTg                                                        <  1:3878530/150‑1 (MQ=255)
                                                                                                ggCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGGTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGc                                                       >  2:5342226/1‑150 (MQ=255)
                                                                                                     aGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaaa                                                  <  2:4471243/150‑1 (MQ=255)
                                                                                                       cGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTAACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGt                                                <  1:4070848/150‑1 (MQ=255)
                                                                                                        ggTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTg                                               <  1:5760239/150‑1 (MQ=255)
                                                                                                         gTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTgg                                              >  2:2453583/1‑150 (MQ=255)
                                                                                                           aaaGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTTGCAGACGCCCGTGATTTATGTCATCGGAGTGGTGGCGCTGGTGGTGTTAATTCACCACCCCTTGATGCCCGTGTATATCGATAATTTGAGTAAAGTGGtt                                            >  2:5171812/1‑150 (MQ=255)
                                                                                                                 aCTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTGTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTAGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGtt                                      >  2:2340332/1‑150 (MQ=255)
                                                                                                                   tGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTa                                    <  2:1429244/150‑1 (MQ=255)
                                                                                                                   tGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTa                                    <  2:2663278/150‑1 (MQ=255)
                                                                                                                    gACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGGGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAg                                   <  2:4849851/150‑1 (MQ=255)
                                                                                                                     cccGTTTGAGTCAGTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAgg                                  <  2:5763939/149‑1 (MQ=255)
                                                                                                                     aCCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAgg                                  <  1:2551902/150‑1 (MQ=255)
                                                                                                                     aCCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAgg                                  >  2:3886177/1‑150 (MQ=255)
                                                                                                                     aCCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAAGTCGCCACCCGTTGATGCCGGTGTATATCGATATTGTGCGTAAAGTGGTTGGGGTTTAgg                                  <  2:79047/150‑1 (MQ=255)
                                                                                                                        gTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGGTGGGGTTTAGgaga                               >  2:3748498/1‑150 (MQ=255)
                                                                                                                         tttGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGAt                              <  1:457367/150‑1 (MQ=255)
                                                                                                                          ttGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATa                             <  1:5422814/150‑1 (MQ=255)
                                                                                                                          ttGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATa                             <  2:429396/150‑1 (MQ=255)
                                                                                                                             aGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAATTGGTTTGGGTTTATGAGATAGta                          >  2:3212142/1‑149 (MQ=255)
                                                                                                                             aGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTc                          <  1:2894593/150‑1 (MQ=255)
                                                                                                                               tCATTATTTGTCCCGCGGTCGTTTGTTGTTGCTGAAGCCGTTGATTTATTCCAACGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATTCCGGTGTATATAGATATTTTTCGTAAATTGGTTTGGGTTTAGGAGATAGTAtt                        >  2:1695196/1‑150 (MQ=255)
                                                                                                                               tCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCtt                        <  1:4457398/150‑1 (MQ=255)
                                                                                                                                 attattTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCCTTTATGCCGGTGTATATCGATATTTTTCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTtgt                      >  2:2770572/1‑150 (MQ=255)
                                                                                                                                  ttattTGTCCCGCGGTCGTTTTTTGTGGCTGACGCCGTTGATTTATGCCATAGGCGTGGTGGCGCTGGTGGTGTTAAATAGCCACCCCTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTTGGGTTTAGGAGATAGTCTtgtg                     >  2:553974/1‑150 (MQ=255)
                                                                                                                                     tttGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCgg                  <  2:5026315/150‑1 (MQ=255)
                                                                                                                                      ttGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCggg                 <  2:4366756/150‑1 (MQ=255)
                                                                                                                                         tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTg              <  1:5096763/150‑1 (MQ=255)
                                                                                                                                         tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTg              <  1:1440500/150‑1 (MQ=255)
                                                                                                                                           ccGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGcc            <  2:3302182/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATg       <  1:5574326/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATg       <  2:4649458/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATg       <  1:2383265/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATg       <  1:199678/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATg       <  1:174065/150‑1 (MQ=255)
                                                                                                                                                  cGTGTGTTGTGGCGGACGCCGTTGCTTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTCAGTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTGGGGGTTGAGGAGATAGTCTTGTGCGGGTTGCCTGATGCg     <  1:1222561/150‑1 (MQ=255)
                                                                                                                                                   gTTTGTTGTGGCTGACGCCGATGATTTCTGCCAACGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGa    >  1:1545866/1‑150 (MQ=255)
                                                                                                                                                    tttGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGAc   >  2:311373/1‑150 (MQ=255)
                                                                                                                                                    tttGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGAc   >  2:3386211/1‑150 (MQ=255)
                                                                                                                                                    tttGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGAc   >  1:1331371/1‑150 (MQ=255)
                                                                                                                                                    tttGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGAc   >  1:3001592/1‑150 (MQ=255)
                                                                                                                                                    tttGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGAc   >  1:559788/1‑150 (MQ=255)
                                                                                                                                                     ttgttgTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACg  <  2:1882508/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACG  >  NC_000913/71850‑72148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: