Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,618,946 G→C 5.4% I405M (ATC→ATG) ‡ uraA ← uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,9460GC 5.4% 642.4 / 22.0 179I405M (ATC→ATG) ‡uraAuracil permease
Reads supporting (aligned to +/- strand):  ref base G (76/92);  new base C (2/8);  total (79/100)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 184 total aligned reads displayed.

CACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTCGA  >  NC_000913/2618797‑2619095
                                                                                                                                                     |                                                                                                                                                      
 acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  1:583867/1‑150 (MQ=255)
 acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  1:2037588/1‑150 (MQ=255)
 acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                       >  1:2037617/1‑150 (MQ=255)
  caAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGAt                                                                                                                                                      <  2:3740463/150‑1 (MQ=255)
  caAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGAt                                                                                                                                                      <  2:3863649/150‑1 (MQ=255)
    aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                    <  2:1762176/150‑1 (MQ=255)
      aTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCaaa                                                                                                                                                  <  2:1971233/150‑1 (MQ=255)
       tCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAAc                                                                                                                                                 <  1:1347661/150‑1 (MQ=255)
           gagaCGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  2:4129815/150‑1 (MQ=255)
           gagaCGTGTACCACAGAATCCCATGCTGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  2:2954360/150‑1 (MQ=255)
           gagaAGTTTACCTCAGTATCCCATTATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCTCGTCTTCTGCGTCCATCTACACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                             <  1:489410/150‑4 (MQ=255)
           gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                             <  1:2781960/150‑1 (MQ=255)
           gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                             <  1:1578509/150‑4 (MQ=255)
             gaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGc                                                                                                                                           <  1:3510538/150‑1 (MQ=255)
             gaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGc                                                                                                                                           <  2:3347265/150‑1 (MQ=255)
               aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                         >  1:2138720/1‑150 (MQ=255)
               aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                         >  1:3792658/1‑150 (MQ=255)
                 tttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGAt                                                                                                                                       >  1:1006623/1‑150 (MQ=255)
                  ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                      <  1:3496340/150‑1 (MQ=255)
                      cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                  <  2:3930244/150‑1 (MQ=255)
                      cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                  <  2:1066054/150‑1 (MQ=255)
                        caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                                >  2:1222232/1‑150 (MQ=255)
                                cATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGc                                                                                                                        >  1:2197550/1‑150 (MQ=255)
                                   gatgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTTTCTGTTATATAAGCGTCTTATGCGTACAGCACCACATCTTAAGGGCGCAGCACCAAGATAATTTTGAAGATAAAAATTAAGCCGATACCGAAAATGTTCGCCac                                                                                                                     >  1:354020/1‑149 (MQ=255)
                                     tgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcg                                                                                                                   <  2:1314312/150‑1 (MQ=255)
                                      gTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcgc                                                                                                                  >  1:2153312/1‑150 (MQ=255)
                                        ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                                <  1:3691841/150‑1 (MQ=255)
                                           aaCCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATAc                                                                                                             >  1:964276/1‑150 (MQ=255)
                                            aCCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATAcc                                                                                                            >  1:3018067/1‑150 (MQ=255)
                                               gggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑ttt                                                                                                        <  1:2451514/150‑1 (MQ=255)
                                                ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑tttt                                                                                                       >  1:3872975/1‑150 (MQ=255)
                                                ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑tttt                                                                                                       <  1:1025133/150‑1 (MQ=255)
                                                ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑tttt                                                                                                       <  2:703318/150‑1 (MQ=255)
                                                 gcagcaACACTGCCCGGTCGGGACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTa                                                                                                      >  1:3584189/1‑150 (MQ=255)
                                                 gcagcaACACTGCCAGGTCGGTACCATATTTGTGTGTTATATCCGCGTATTCAGCGTCCAGCACCACTTCTTCAGGGCGCCGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCCACGATGGTAGCCAGCGCCAAACC‑Ttata                                                                                                      >  1:3584225/1‑150 (MQ=255)
                                                             cccGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCcgccgcgc                                                                                          <  2:2697723/150‑1 (MQ=255)
                                                                       aCATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTAc                                                                                <  2:1246957/150‑1 (MQ=255)
                                                                           tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                            >  1:1912946/1‑150 (MQ=255)
                                                                            attTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTcg                                                                           <  2:180612/150‑1 (MQ=255)
                                                                               tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                        <  1:158594/150‑1 (MQ=255)
                                                                               tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCACCGCGCCGATGTTTACCTTCgcgc                                                                        <  1:814412/150‑1 (MQ=255)
                                                                                 tCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCa                                                                      <  1:3344394/150‑1 (MQ=255)
                                                                                  cTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCAc                                                                     <  1:1133470/150‑1 (MQ=255)
                                                                                     ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                  <  2:141700/150‑1 (MQ=255)
                                                                                     ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                  >  2:1463810/1‑150 (MQ=255)
                                                                                      tataTCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGAc                                                                 <  2:2731129/150‑1 (MQ=255)
                                                                                         atCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGcc                                                              <  2:1850287/150‑1 (MQ=255)
                                                                                                tcttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                       <  2:2439147/150‑1 (MQ=255)
                                                                                                tcgtctGCGTCCAGCACCACGTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                       <  2:3856144/147‑1 (MQ=255)
                                                                                                 cttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCa                                                      >  1:930698/1‑150 (MQ=255)
                                                                                                  ttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGAGGTTTACCTTCGCGCCACTGACGCCGATGATCaa                                                     >  2:1908626/1‑150 (MQ=255)
                                                                                                    ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                   >  1:2741053/1‑150 (MQ=255)
                                                                                                     tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                  <  2:3062936/150‑1 (MQ=255)
                                                                                                     tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                  >  1:3768965/1‑150 (MQ=255)
                                                                                                     tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                  >  2:2612543/1‑150 (MQ=255)
                                                                                                       cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                                <  2:3291589/150‑1 (MQ=255)
                                                                                                       cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                                <  2:964276/150‑1 (MQ=255)
                                                                                                       cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                                <  2:1021049/150‑1 (MQ=255)
                                                                                                            aGCACCACTTCTTCAGGGCGCAGCACTCTGATCTGTTTGTAGATCAAAGTTATGCCGTTAGCGACGATGGTCGCCAGCGCCATAGC‑TTTTATCTCCGCCGCGCCGTTGTTTTCCTTCGCGCCTCTGACGGCGATGTTCTTATTCTCGGaa                                           <  2:86254/150‑1 (MQ=255)
                                                                                                            aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCTCCAGCGCCATACCATTTAAACACCCCCGCCCCAATTTTAACGTACGCCCCACTGACCCCAAGAAGCAAAATCACGGa                                            >  1:2484300/1‑150 (MQ=255)
                                                                                                            aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                           >  1:2574449/1‑150 (MQ=255)
                                                                                                            aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                           >  1:2688114/1‑150 (MQ=255)
                                                                                                             gCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAg                                          >  2:3697553/1‑150 (MQ=255)
                                                                                                               accacTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTc                                        <  1:2443773/150‑1 (MQ=255)
                                                                                                                ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                       <  2:2822365/150‑1 (MQ=255)
                                                                                                                ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                       <  2:2385587/150‑1 (MQ=255)
                                                                                                                ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                       <  2:2257140/150‑1 (MQ=255)
                                                                                                                  acTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAgg                                     >  1:1207140/1‑150 (MQ=255)
                                                                                                                   cttcttCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGa                                    >  1:445222/1‑150 (MQ=255)
                                                                                                                    ttcttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGGTTACCGTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGAt                                   >  2:2516750/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACTATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                <  2:16894/150‑1 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCAGGGAAGTCAGGATCAg                                >  2:730719/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                >  2:106336/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                >  2:40732/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                >  1:2512902/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                >  1:3188374/1‑150 (MQ=255)
                                                                                                                       ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                                >  1:3528342/1‑150 (MQ=255)
                                                                                                                            ggCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCt                           <  2:1208388/150‑1 (MQ=255)
                                                                                                                              cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                         <  1:1694359/150‑1 (MQ=255)
                                                                                                                              cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                         <  1:1694634/150‑1 (MQ=255)
                                                                                                                              cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                         <  1:2369493/150‑1 (MQ=255)
                                                                                                                              cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                         <  2:152869/150‑1 (MQ=255)
                                                                                                                               gcagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgtg                        <  2:3284587/150‑1 (MQ=255)
                                                                                                                                cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                       <  1:2780412/150‑1 (MQ=255)
                                                                                                                                cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                       <  2:4009579/150‑1 (MQ=255)
                                                                                                                                  gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                     <  1:3967786/150‑1 (MQ=255)
                                                                                                                                  gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                     >  1:1251180/1‑150 (MQ=255)
                                                                                                                                  gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                     <  2:3784247/150‑1 (MQ=255)
                                                                                                                                   cacaCTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCttt                    <  1:1501792/150‑1 (MQ=255)
                                                                                                                                      acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt                 >  2:217563/1‑150 (MQ=255)
                                                                                                                                         cttaaGTTTGATCTTCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTa              >  2:2898769/5‑150 (MQ=255)
                                                                                                                                           tCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTACCTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAAt            <  2:3405140/150‑1 (MQ=255)
                                                                                                                                            cAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATc           <  1:90942/150‑1 (MQ=255)
                                                                                                                                              gTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCa         <  2:3538021/150‑1 (MQ=255)
                                                                                                                                                ttGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACt       >  1:2966780/1‑150 (MQ=255)
                                                                                                                                                 tGAAGATCAATCTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGTTCAAATTCACGGTGGTCAGGATCAGGTTCTGTGCTTTGTTGTATTCCTCtt      >  2:3097849/1‑150 (MQ=255)
                                                                                                                                                 tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACtt      >  2:2127383/1‑150 (MQ=255)
                                                                                                                                                 tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTATGTTGTAATCCACtt      >  2:1427117/1‑150 (MQ=255)
                                                                                                                                                     gATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttcga  <  2:424085/150‑1 (MQ=255)
                                                                                                                                                     gATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttcga  >  2:2274775/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTCGA  >  NC_000913/2618797‑2619095

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: ‑