Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,047,524 C→G 5.2% M506I (ATG→ATC)  gcvP ← glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,047,5240CG 5.2% 249.6 / 2.4 74M506I (ATG→ATC) gcvPglycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base G (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTG  >  NC_000913/3047378‑3047670
                                                                                                                                                  |                                                                                                                                                  
cagTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGatca                                                                                                                                                   <  5:260617/148‑1 (MQ=255)
       ccagcTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCt                                                                                                                                             >  6:54014/1‑147 (MQ=255)
           cTGCGCAATCATCTGCTGATAACCTTCGGACTGCTCCGGCGGGCAGAAAGGGTGCATTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGAtt                                                                                                                                         >  3:135867/1‑147 (MQ=255)
              cgcAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAg                                                                                                                                      <  6:324656/147‑1 (MQ=255)
               gcAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAgcg                                                                                                                                    >  3:90115/1‑148 (MQ=255)
                  atcatcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAgcgc                                                                                                                                   >  4:90758/1‑146 (MQ=255)
                   tcatcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAgcgc                                                                                                                                   >  8:220042/1‑145 (MQ=255)
                    catcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAgcgc                                                                                                                                   >  2:242443/1‑144 (MQ=255)
                     ctcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCGGCGCCAGa                                                                                                                               <  5:54014/146‑1 (MQ=255)
                     atcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGa                                                                                                                               <  4:135867/147‑1 (MQ=255)
                      tcTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATc                                                                                                                             >  7:220302/1‑148 (MQ=255)
                         gctgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCtt                                                                                                                           >  2:308681/1‑147 (MQ=255)
                          ctgATAACCTTCTGCCTTCTCCGGCGTGAAGAACGGGTGCAGTTCGGCAAATTACGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCttt                                                                                                                          <  4:260128/147‑1 (MQ=255)
                          ctgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCttt                                                                                                                          <  7:220042/147‑1 (MQ=255)
                          ctgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCt                                                                                                                            <  4:90115/145‑1 (MQ=255)
                          ctgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGGGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCt                                                                                                                            <  6:77812/145‑1 (MQ=255)
                           tgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCTGCGCCAGATCTTTAc                                                                                                                        <  5:250941/148‑1 (MQ=255)
                           tgATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATAATCGCCTGATTCAGCGCCAGATCTTTAc                                                                                                                        >  6:250941/1‑148 (MQ=255)
                             aTAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGc                                                                                                                      >  7:323681/1‑148 (MQ=255)
                                        ccTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTg                                                                                                           <  3:39895/148‑1 (MQ=255)
                                         cTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGc                                                                                                          <  5:438222/148‑1 (MQ=255)
                                            cTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCat                                                                                                        <  5:66190/147‑1 (MQ=255)
                                            cTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCat                                                                                                        <  1:308681/147‑1 (MQ=255)
                                              ccGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCata                                                                                                       >  3:54453/1‑146 (MQ=255)
                                               cggcggGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCatat                                                                                                      <  1:211659/146‑1 (MQ=255)
                                                ggcggGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAgc                                                                                                   >  3:39198/1‑148 (MQ=255)
                                                       aGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCTcatcat                                                                                            <  4:39198/148‑1 (MQ=255)
                                                        gAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCAtt                                                                                           <  3:90758/148‑1 (MQ=255)
                                                          aCGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCAttt                                                                                          <  7:371578/147‑1 (MQ=255)
                                                          aCGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTc                                                                                         >  3:86007/1‑148 (MQ=255)
                                                              gTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGtt                                                                                     <  1:242443/148‑1 (MQ=255)
                                                              gTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGtt                                                                                     >  1:66594/1‑148 (MQ=255)
                                                               tGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGttt                                                                                    >  4:202109/1‑148 (MQ=255)
                                                               tGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGttt                                                                                    <  2:66594/148‑1 (MQ=255)
                                                                 cAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCg                                                                                  <  8:323681/148‑1 (MQ=255)
                                                                    ttCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCtg                                                                               <  3:380009/148‑1 (MQ=255)
                                                                         cAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTgg                                                                            >  8:352312/1‑146 (MQ=255)
                                                                         cAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTa                                                                          <  6:416900/148‑1 (MQ=255)
                                                                               ccGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATCGCGCATCATTTCGGTTTCGCTGTGGTAGCGAt                                                                     <  4:86007/147‑1 (MQ=255)
                                                                                        tGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCtgt                                                                              >  6:424263/1‑129 (MQ=255)
                                                                                        tGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCtgt                                                                              <  5:424263/129‑1 (MQ=255)
                                                                                              ggATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGt                                                     >  7:282684/1‑148 (MQ=255)
                                                                                               gATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGt                                                     >  5:377833/1‑147 (MQ=255)
                                                                                               gATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTg                                                    <  7:352312/148‑1 (MQ=255)
                                                                                               gATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTg                                                    <  7:209339/148‑1 (MQ=255)
                                                                                                   atcTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGa                                                >  3:297710/1‑148 (MQ=255)
                                                                                                     ctcGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCAGCATTTCGGTTTCGCTGTGGTAGAGATTAAACACCTGATGGGTGAGGAtt                                              >  6:29188/1‑148 (MQ=255)
                                                                                                        ggcggcggcGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAg                                                                         <  8:61066/118‑1 (MQ=255)
                                                                                                        ggcggcggcGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAg                                                                         >  7:61066/1‑118 (MQ=255)
                                                                                                        ggcggcggcGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGAttt                                             <  4:297710/146‑1 (MQ=255)
                                                                                                                cgTGCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTcgcg                                   <  8:196250/148‑1 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAgcgc                                                                                                 >  1:6731/1‑83 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAgcgc                                                                                                 <  2:6731/83‑1 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCg                                                                                  <  1:323306/98‑1 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCg                                                                                  >  2:323306/1‑98 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGAttt                                             >  4:160295/1‑135 (MQ=255)
                                                                                                                   tCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGAttt                                             <  3:160295/135‑1 (MQ=255)
                                                                                                                    cAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCgcag                                >  6:420510/1‑147 (MQ=255)
                                                                                                                       tttCATGGTGCAGGAACCCAGCTGGATGATCGCCTGATTAAGCGCCAGATCTTTACGCTCAAGCGAGTGCATATAGAGCATCATTTCCGTTTCGCTGTGGTAGCGATTAAACACAGTATTGGTGAGGCTTTCGTCGTCGCGAAGAATc                            >  6:164425/1‑148 (MQ=255)
                                                                                                                       tttCATGGTGCAGGAACCCAGCGGGATGATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATc                            >  6:164421/1‑148 (MQ=255)
                                                                                                                        ttCATGGTGCAGGAACCCAGCGGGGTCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCg                           >  3:381587/1‑148 (MQ=255)
                                                                                                                           aTGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAg                        >  4:95593/1‑148 (MQ=255)
                                                                                                                             ggTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGc                      >  7:377735/1‑148 (MQ=255)
                                                                                                                             ggTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGc                      >  3:208531/1‑148 (MQ=255)
                                                                                                                                gCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAgg                                                 <  6:62815/118‑1 (MQ=255)
                                                                                                                                gCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAgg                                                 >  5:62815/1‑118 (MQ=255)
                                                                                                                                 cAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTgg                   <  5:420510/147‑1 (MQ=255)
                                                                                                                                   ggAACCCAGCGGGGTCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATa                <  4:381587/148‑1 (MQ=255)
                                                                                                                                    gAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGa                  >  5:66210/1‑145 (MQ=255)
                                                                                                                                       cccAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGc            >  7:388849/1‑148 (MQ=255)
                                                                                                                                            cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTg       >  2:152005/1‑148 (MQ=255)
                                                                                                                                            cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTg       >  6:92420/1‑148 (MQ=255)
                                                                                                                                              ggATGATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTc     <  5:164421/148‑1 (MQ=255)
                                                                                                                                              ggATGATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCAGCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCAGCGCAGGCTGGATAGAGCGGCTggc     <  5:164425/148‑3 (MQ=255)
                                                                                                                                               gATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTATCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCg    <  4:54453/148‑1 (MQ=255)
                                                                                                                                                 tcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTg  <  8:220302/148‑1 (MQ=255)
                                                                                                                                                 tcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGGCGTg  <  3:95593/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTGCAGTTCGGCAAATTCCGGCCAGGTGATTGGGATCATCTCGGCGGCGGCGTTCAGTTTCATGGTGCAGGAACCCAGCGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTG  >  NC_000913/3047378‑3047670

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