Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,261,471 A→C 8.8% I135M (ATT→ATG)  prs ← phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,261,4710AC 8.8% 227.9 / 13.3 80I135M (ATT→ATG) prsphosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base C (3/4);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

CAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAGGGGATAACAG  >  NC_000913/1261327‑1261616
                                                                                                                                                |                                                                                                                                                 
cagcTTAGCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGCATCt                                                                                                                                               <  4:191696/149‑1 (MQ=255)
       gCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGACACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGc                                                                                                                                        <  2:208924/149‑1 (MQ=255)
       gCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGc                                                                                                                                        >  1:300940/1‑149 (MQ=255)
                  ggCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCa                                                                                                                             <  7:304141/149‑1 (MQ=255)
                      cgcacAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGt                                                                                                                         >  8:56980/1‑149 (MQ=255)
                         acaACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAg                                                                                                                      >  5:353646/1‑149 (MQ=255)
                            aCGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGTTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGACTCTGTTCAGCGTGCAGATCCACTGTCAGcac                                                                                                                   >  6:356512/1‑149 (MQ=255)
                              gccgccGATGTCCGGAGAAACCACAATTGGGTTATCCAGATACAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCacac                                                                                                                 <  7:392787/149‑1 (MQ=255)
                               ccgccgATGTCCGGAGAAACCACAATTGGTTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCTGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACg                                                                                                                <  1:230297/149‑1 (MQ=255)
                               ccgccgATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACg                                                                                                                <  5:55290/149‑1 (MQ=255)
                                    gATGTCCGGAGAATCCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCaa                                                                                                           <  5:112914/149‑1 (MQ=255)
                                        tCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAacac                                                                                                        <  2:54409/148‑1 (MQ=255)
                                        tCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAacac                                                                                                        <  1:162729/148‑1 (MQ=255)
                                        tCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACAcc                                                                                                       <  8:322789/149‑1 (MQ=255)
                                         ccGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCg                                                                                                      >  5:160906/1‑149 (MQ=255)
                                          cGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGa                                                                                                     >  8:114612/1‑149 (MQ=255)
                                            gagaAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCATCCAGCAGGATCGCGATACCAAATACGTAATCAACCGGAACTTCGAAGAAAACCTGAATCTCTTCAGCGTGCAGATCCACTGTAAGCACAATGTCAAAACCGACg                                                                                                   >  2:344484/1‑149 (MQ=255)
                                                 aCCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTg                                                                                                >  2:107182/1‑147 (MQ=255)
                                                   cacaATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGga                                                                                              >  6:24394/1‑147 (MQ=255)
                                                   cacaATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGGCAACACCGACGCTGgaga                                                                                            <  8:285162/149‑1 (MQ=255)
                                                    acaATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGagag                                                                                           >  3:239734/1‑149 (MQ=255)
                                                     caATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGTTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGACTCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAgg                                                                                          <  5:356512/149‑1 (MQ=255)
                                                       aTTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGTGCTACCAAATACGTTATCAAACGGAACGTAGAATAAACCCTGAATCTGTTAATCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGaa                                                                                        <  1:344484/149‑1 (MQ=255)
                                                        ttGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAg                                                                                       >  2:279729/1‑149 (MQ=255)
                                                             ttATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCGACCAATTACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACTCGGTCAACACCGACGCTGGAGAGGAAGTCTGc                                                                                  <  6:388479/149‑1 (MQ=255)
                                                                     aTTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCGTGAATCTGTGCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACtt                                                                          <  2:268893/149‑1 (MQ=255)
                                                                         aGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCg                                                                       <  5:24394/148‑1 (MQ=255)
                                                                             gcagcaTGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACCCCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGt                                                                   >  4:279877/1‑148 (MQ=255)
                                                                                 caTGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTg                                                              >  5:372802/1‑149 (MQ=255)
                                                                                 caTGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTg                                                              >  6:66310/1‑149 (MQ=255)
                                                                                  aTGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAACCCCTGAATCTGTTCAGCGTGCAGATCCACAGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTgg                                                             <  4:200995/149‑1 (MQ=255)
                                                                                   tGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTgg                                                             >  8:363412/1‑148 (MQ=255)
                                                                                    gTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGCAGACACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGAGATTGGTa                                                           >  5:228174/1‑149 (MQ=255)
                                                                                     tCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAAGCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTac                                                          <  4:222510/149‑1 (MQ=255)
                                                                                     tCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTg                                                              >  4:294535/1‑145 (MQ=255)
                                                                                     tCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTg                                                              <  3:294535/145‑1 (MQ=255)
                                                                                     tCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTa                                                           <  3:374009/148‑1 (MQ=255)
                                                                                     tCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTa                                                           <  1:107182/148‑1 (MQ=255)
                                                                                       ttCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTacac                                                        >  7:345962/1‑149 (MQ=255)
                                                                                       ttCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTacac                                                        >  4:303971/1‑149 (MQ=255)
                                                                                       ttCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTacac                                                        >  6:109548/1‑149 (MQ=255)
                                                                                         ccAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTacac                                                        >  7:2188/1‑147 (MQ=255)
                                                                                          cagcagGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGa                                                      <  6:160906/148‑1 (MQ=255)
                                                                                            gcagGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCg                                                   <  6:260974/149‑1 (MQ=255)
                                                                                               ggATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGGCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGaa                                                <  8:405533/149‑1 (MQ=255)
                                                                                                gATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAAc                                               >  8:140396/1‑149 (MQ=255)
                                                                                                   cGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACgg                                             <  6:372802/148‑1 (MQ=255)
                                                                                                     ggCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAAACACTTTCGCAGTGATTGGTACACGAGCGGAACGGAcg                                          >  7:114683/1‑149 (MQ=255)
                                                                                                        tACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGa                                       >  8:5764/1‑149 (MQ=255)
                                                                                                          ccAAATACGTTATCAACCGGATAGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGAc                                      <  6:228174/148‑1 (MQ=255)
                                                                                                          ccAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGa                                       <  4:239734/147‑1 (MQ=255)
                                                                                                           cAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACCCGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACgg                                    >  6:76465/1‑149 (MQ=255)
                                                                                                            aaaTACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGAc                                           >  1:311253/1‑141 (MQ=255)
                                                                                                            aaaTACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGAc                                           <  2:311253/141‑1 (MQ=255)
                                                                                                             aaTACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTc                                  >  1:304538/1‑149 (MQ=255)
                                                                                                               tACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCt                                >  4:2699/1‑149 (MQ=255)
                                                                                                               tACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCt                                >  8:120532/1‑149 (MQ=255)
                                                                                                                 cGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTg                               <  7:363412/148‑1 (MQ=255)
                                                                                                                 cGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGGCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCCGAACGGACGCGACGGTCCTg                               <  2:300940/148‑1 (MQ=255)
                                                                                                                   ttATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCt                                >  2:201199/1‑145 (MQ=255)
                                                                                                                   ttATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCt                                <  1:201199/145‑1 (MQ=255)
                                                                                                                      tCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGgcgcg                          >  4:356562/1‑148 (MQ=255)
                                                                                                                         aCCGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCa                        >  3:202590/1‑147 (MQ=255)
                                                                                                                         aCCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCa                        >  4:320846/1‑147 (MQ=255)
                                                                                                                         aCCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATa                      <  5:66310/149‑1 (MQ=255)
                                                                                                                          ccGGAACGTCGAAGAAACCCTGCATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGCGCGGAACGGACGCGACGGTCCTGGCGCGCATAg                     >  3:298792/1‑149 (MQ=255)
                                                                                                                               aCGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCaaa                <  3:356562/149‑1 (MQ=255)
                                                                                                                                cGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAg               >  2:193707/1‑149 (MQ=255)
                                                                                                                                  tCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGaa                                                                                        >  6:91988/1‑74 (MQ=255)
                                                                                                                                  tCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGaa                                                                                        <  5:91988/74‑1 (MQ=255)
                                                                                                                                    gaagaaACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAaggagg           >  2:91391/1‑149 (MQ=255)
                                                                                                                                    gaagaaACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAgg           <  2:304538/149‑1 (MQ=255)
                                                                                                                                    gaagaaACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAg            <  7:56980/148‑1 (MQ=255)
                                                                                                                                      agaaACCCTGAATCTGTTCAGAGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGAAGTGCTTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAggg          >  8:300599/1‑148 (MQ=255)
                                                                                                                                       gaaACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAGGGGa        <  1:279729/149‑1 (MQ=255)
                                                                                                                                        aaaCCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAGGGGa        >  2:218086/1‑148 (MQ=255)
                                                                                                                                           cccTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGTAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGTAACGGACGCTACTGTACTGGCGCGCATAGCCAAATTAggggctca     >  4:306779/1‑144 (MQ=255)
                                                                                                                                            ccTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGTAAGTCTGCAACCACTTTCGCAGTGATTGGTACAAGAGCGGAACGGACGCGCCGGTCCTGGCGCGCATAGCAAAAGAAGGGGATAACa   >  6:64137/1‑149 (MQ=255)
                                                                                                                                            ccTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCGGGAGAGGAAGTCTGCAACCACTTTCGCAGTGCGTGGTACACGAGCCGACCGGACGCGACGGGCCTGGCGCGCATAGCCTAAGTAGGGGATaa     >  5:231311/1‑147 (MQ=255)
                                                                                                                                             cTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAGGGGATAACAg  >  4:365933/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAGCTTAGCGATAGCGCGGGCACGCACAACGCCGCCGATGTCCGGAGAAACCACAATTGGGTTATCCAGATTCAGCTGCAGCATGTCTTCCAGCAGGATCGGGCTACCAAATACGTTATCAACCGGAACGTCGAAGAAACCCTGAATCTGTTCAGCGTGCAGATCCACTGTCAGCACACGGTCAACACCGACGCTGGAGAGGAAGTCTGCAACCACTTTCGCAGTGATTGGTACACGAGCGGAACGGACGCGACGGTCCTGGCGCGCATAGCCAAAGTAGGGGATAACAG  >  NC_000913/1261327‑1261616

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