Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,288,332 A→G 5.1% L98P (CTT→CCT)  purU ← formyltetrahydrofolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,288,3320AG 5.1% 252.0 / 2.9 79L98P (CTT→CCT) purUformyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base G (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TTGCGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATA  >  NC_000913/1288188‑1288478
                                                                                                                                                |                                                                                                                                                  
ttGCGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGc                                                                                                                                                 <  2:316266/148‑1 (MQ=255)
  gCGGGTTAACCCTTCATGGCTTACCAGATCAATCGGAATATCAAAACGCTCAAGCAGAGAAGGTAAAGTATCGTGGTTACCAATAACTGGCGCGATTTCGACTTCCAGGCCGCAGTAATTGGCTTTCATCAACAAATCGCCAAGGCa                                                                                                                                                <  7:118405/147‑1 (MQ=255)
   cGGGTTAACCCTTCATGGCTTACCATCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCATGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCaa                                                                                                                                               <  5:375996/147‑1 (MQ=255)
   cGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAAt                                                                                                                                              <  8:202503/148‑1 (MQ=255)
           cccTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGcttct                                                                                                                                      >  3:350607/1‑148 (MQ=255)
           cccTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGcttct                                                                                                                                      <  8:187851/148‑1 (MQ=255)
               tCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATgcg                                                                                                                                           >  4:284881/1‑139 (MQ=255)
               tCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATgcg                                                                                                                                           <  3:284881/139‑1 (MQ=255)
               tCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGcttctt                                                                                                                                     >  2:322070/1‑145 (MQ=255)
               tCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATAGTGGTTACCAATAACTGCCGCGATTTAGACATCCAGGCCGCCGTAATTGGCTTTCATAAACAAATCGCCAAGGCAATGCGCTTCTTTAg                                                                                                                                  >  3:388466/1‑148 (MQ=255)
                gaTGGCTTACGAGCTCATACGGAATATCATAACGCTCAACTAGTGATCGTATAGTATCGTGGTTACCAATAACTGCAGCGATTTCGACATCCAGGCCGCCGTATTTGGCTTTGATCAACAAATCGCCAAGGCAATGCGCTTCTTTAg                                                                                                                                  <  2:383198/146‑1 (MQ=255)
                cATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCTCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTa                                                                                                                                   <  4:148716/146‑1 (MQ=255)
                cATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCttt                                                                                                                                    <  1:322070/145‑1 (MQ=255)
                cATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAg                                                                                                                                  <  8:188587/147‑1 (MQ=255)
                cATGGCTTACCAGCTCAAACGGAATATCAAAACCCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTa                                                                                                                                   <  2:198867/146‑1 (MQ=255)
                 aTGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGt                                                                                                                                 <  4:350607/147‑1 (MQ=255)
                    gCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGAc                                                                                                                              >  2:201788/1‑147 (MQ=255)
                                  aCGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCa                                                                                                                      <  1:191118/141‑1 (MQ=255)
                                  aCGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCa                                                                                                                      >  2:191118/1‑141 (MQ=255)
                                     gAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGCCATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGcg                                                                                                            <  4:287668/148‑1 (MQ=255)
                                     gAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGc                                                                                                             >  4:198362/1‑147 (MQ=255)
                                         atCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGcgacga                                                                                                        >  2:104783/1‑148 (MQ=255)
                                            aaaaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAGGGCAATGCGCTTCTTTAGTGACCAGAATCAc                                                                                                                     <  4:411376/132‑1 (MQ=255)
                                            aaaaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAGGGCAATGCGCTTCTTTAGTGACCAGAATCAc                                                                                                                     <  4:404934/132‑1 (MQ=255)
                                            aaaaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAGGGCAATGCGCTTCTTTAGTGACCAGAATCAc                                                                                                                     >  3:411376/1‑132 (MQ=255)
                                            aaaaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAGGGCAATGCGCTTCTTTAGTGACCAGAATCAc                                                                                                                     >  3:404934/1‑132 (MQ=255)
                                             aaaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGTGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCgg                                                                                                    >  3:65141/1‑148 (MQ=255)
                                              aaCGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGc                                                                                                   <  6:4994/148‑1 (MQ=255)
                                                     aaCCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTc                                                                                            <  5:71463/148‑1 (MQ=255)
                                                      aCCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGAtt                                                                                             >  1:358462/1‑146 (MQ=255)
                                                      aCCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGAGACGACCGGCAGGAtt                                                                                             >  6:391664/1‑146 (MQ=255)
                                                        cAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATAACTATCCGGCGCAGACCGGCAGGATTCAg                                                                                          >  1:352165/1‑147 (MQ=255)
                                                           agaACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTca                                                                                      <  3:198362/148‑1 (MQ=255)
                                                             aaCGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTcacg                                                                                    <  4:113615/148‑1 (MQ=255)
                                                                gTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCGACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCAGCAGGATTCAGCTcacgcac                                                                                 >  1:9680/1‑148 (MQ=255)
                                                                   aaGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACgg                                                                               <  2:358462/147‑1 (MQ=255)
                                                                      tATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCATAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGcc                                                                           <  1:304444/148‑1 (MQ=255)
                                                                       aTCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCttt                                                                                                                                    >  8:110690/1‑90 (MQ=255)
                                                                       aTCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCttt                                                                                                                                    <  7:110690/90‑1 (MQ=255)
                                                                         cGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCtt                                                                         <  1:104783/147‑1 (MQ=255)
                                                                         cGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACAATACGCCCACGACCGCCAGGATACAGCTCACGCACGGAGCCTTc                                                                        >  3:138559/1‑148 (MQ=255)
                                                                            ggTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTa                                                                                                                   <  2:207520/102‑1 (MQ=255)
                                                                            ggTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTa                                                                                                                   >  1:207520/1‑102 (MQ=255)
                                                                            ggTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAACCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTgg                                                                     >  2:374908/1‑148 (MQ=255)
                                                                                aCCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCaag                                                                 >  5:366471/1‑147 (MQ=255)
                                                                                aCCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCaa                                                                  <  6:366471/147‑1 (MQ=255)
                                                                                 ccAATAACTGCCGCGATTTCGACATCCAGGCCGCCCTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATg                                                                >  5:406577/1‑148 (MQ=255)
                                                                                   aaTAACTGCCTCGATTTCGACATCCAGGCCGCCGTAATTGGATTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTATTGACCTGAATCACTATCCTGAGACGACCGGCAGGATTCAGCTCACGCACGGATCCTTCTGGCAATgcg                                                              <  5:391664/148‑1 (MQ=255)
                                                                                       aCTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTAt                                                          <  2:425562/148‑1 (MQ=255)
                                                                                       aCTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAACCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTAt                                                          <  1:374908/148‑1 (MQ=255)
                                                                                            cgcgATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCCCTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATc                                                         >  7:332527/1‑144 (MQ=255)
                                                                                            cgcgATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATc                                                         >  5:352188/1‑144 (MQ=255)
                                                                                            cgcgATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCgag                                                      >  4:38518/1‑147 (MQ=255)
                                                                                                    cGACATCCAGGCCGCCGTAATTGGCTTTCATCGACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCAGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCca                                              <  2:9680/147‑1 (MQ=255)
                                                                                                    cGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATg                                                                >  8:295316/1‑129 (MQ=255)
                                                                                                    cGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATg                                                                <  7:295316/129‑1 (MQ=255)
                                                                                                      aCATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCcagca                                           <  4:65141/148‑1 (MQ=255)
                                                                                                           cAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCa                                                                                                                                                >  4:145716/1‑42 (MQ=255)
                                                                                                           cAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCa                                                                                                                                                <  3:145716/42‑1 (MQ=255)
                                                                                                           cAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGt                                       <  3:38518/147‑1 (MQ=255)
                                                                                                                  ccgTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCAt                               >  6:84577/1‑148 (MQ=255)
                                                                                                                       aTTGGCTTTCATCATCAAATCGCCAAGGGAATCCTCTTCTTTATTGACCATAATCAATATCCGGCGACGCCCGTCAGTATTCAGCTCACGCACGTATACTTCTGGAAATGCGATATCGAGATCAGACAGCAGGGTGGAATCATTaaaa                          <  1:70900/148‑1 (MQ=255)
                                                                                                                       aTTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTaaaa                          <  1:201788/148‑1 (MQ=255)
                                                                                                                       aTTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTaaa                           <  4:169533/147‑1 (MQ=255)
                                                                                                                       aTTGGCTTTCATAAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTaaa                           <  4:388466/147‑1 (MQ=255)
                                                                                                                              ttCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCtt                   >  6:239129/1‑148 (MQ=255)
                                                                                                                              ttCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCtt                   >  8:365664/1‑148 (MQ=255)
                                                                                                                               tcatcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCtt                   <  7:365664/147‑1 (MQ=255)
                                                                                                                               tcatcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCt                    <  1:262816/146‑1 (MQ=255)
                                                                                                                                 atcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCa                >  5:251842/1‑148 (MQ=255)
                                                                                                                                  tcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCgg                                                                                                    >  8:51218/1‑63 (MQ=255)
                                                                                                                                  tcaACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCgg                                                                                                    <  7:51218/63‑1 (MQ=255)
                                                                                                                                   caacaaATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGt              >  2:59297/1‑148 (MQ=255)
                                                                                                                                     acaaATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTc            >  4:467/1‑148 (MQ=255)
                                                                                                                                      caaATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTcc           >  3:148639/1‑148 (MQ=255)
                                                                                                                                       aaaTCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCg          >  1:398259/1‑148 (MQ=255)
                                                                                                                                        aaTCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGt         <  6:406577/148‑1 (MQ=255)
                                                                                                                                         aTCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCa                <  2:374199/140‑1 (MQ=255)
                                                                                                                                         aTCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCa                >  1:374199/1‑140 (MQ=255)
                                                                                                                                               aaGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCAt   >  7:23912/1‑147 (MQ=255)
                                                                                                                                                aGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATa  >  6:5437/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TTGCGGGTTAACCCTTCATGGCTTACCAGCTCAAACGGAATATCAAAACGCTCAACCAGAGAACGTAAAGTATCGTGGTTACCAATAACTGCCGCGATTTCGACATCCAGGCCGCCGTAATTGGCTTTCATCAACAAATCGCCAAGGCAATGCGCTTCTTTAGTGACCAGAATCACTATCCGGCGACGACCGGCAGGATTCAGCTCACGCACGGAGCCTTCTGGCAATGCGCTATCGAGATCCGCCAGCAGGGTGGAATCATTAAAAATCCCTTCCAGTTCCGTGCGCATA  >  NC_000913/1288188‑1288478

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