Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,916,621 T→C 5.1% F122L (TTT→CTT)  yebS → inner membrane PqiA domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,916,6210TC 5.1% 241.0 / 2.7 78F122L (TTT→CTT) yebSinner membrane PqiA domain protein
Reads supporting (aligned to +/- strand):  ref base T (31/43);  new base C (2/2);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

TTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGATATCTACCTGGTCGG  >  NC_000913/1916477‑1916769
                                                                                                                                                |                                                                                                                                                    
ttCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCttttt                                                                                                                                                  >  8:347380/1‑149 (MQ=255)
   aCCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCt                                                                                                                                                <  7:82203/148‑1 (MQ=255)
   aCCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTg                                                                                                                                               <  5:372948/149‑1 (MQ=255)
     cATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCg                                                                                                                                             >  4:183381/1‑149 (MQ=255)
         cTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTa                                                                                                                                          <  8:175505/148‑1 (MQ=255)
          tGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTAtt                                                                                                                                        >  1:347951/1‑149 (MQ=255)
          tGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTAtt                                                                                                                                        >  6:74203/1‑149 (MQ=255)
           gTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTg                                                                                                                                       >  4:193049/1‑149 (MQ=255)
            ttGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAACCAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTgg                                                                                                                                      >  4:299349/1‑149 (MQ=255)
              gCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTg                                                                                                                                       <  6:249681/146‑1 (MQ=255)
              gCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTg                                                                                                                                       >  5:249681/1‑146 (MQ=255)
               cTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTg                                                                                                                                    <  1:302528/148‑1 (MQ=255)
               cTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGc                                                                                                                                   <  8:254061/149‑1 (MQ=255)
                ttATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGCTGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGcc                                                                                                                                  >  4:133516/1‑149 (MQ=255)
                  aTGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGcc                                                                                                                                  <  1:45810/147‑1 (MQ=255)
                    gCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATAGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTACTGCGTTATTGGTGCCCCGc                                                                                                                              >  6:168211/1‑149 (MQ=255)
                     ccGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGCCGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGc                                                                                                                              <  5:357609/148‑1 (MQ=255)
                       gTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAACCAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCa                                                                                                                           >  6:214011/1‑149 (MQ=255)
                        tttGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCAt                                                                                                                          >  5:49337/1‑149 (MQ=255)
                            ccTGGGGTGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTg                                                                                                                      <  4:216399/149‑1 (MQ=255)
                            ccTGGGGCGAACCGCTGTTGCATATCTGGCTTTGAGGCATTCGTATTGACGCCAACGTAATGCAAGGCATCTGGCAAATGACCAAAGAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTg                                                                                                                      <  6:25463/149‑1 (MQ=255)
                            ccTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGGCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTg                                                                                                                      <  5:74203/149‑1 (MQ=255)
                                gggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGt                                                                                                                    <  3:193049/147‑1 (MQ=255)
                                gggCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGt                                                                                                                    <  7:336647/147‑1 (MQ=255)
                                        cGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCa                                                                                                            >  7:66168/1‑147 (MQ=255)
                                        cGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATa                                                                                                          <  5:152781/149‑1 (MQ=255)
                                          cTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGc                                                                                                        <  3:378599/149‑1 (MQ=255)
                                             ttGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCtta                                                                                                     <  2:467983/149‑1 (MQ=255)
                                             atGCATATCTGGCATTGATGCATTCGTAATGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATCACGGGGTCGATGGTCTATATCTTCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATatc                                                                                                        <  6:114293/145‑3 (MQ=255)
                                               gCATATCTGGCTGGTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAATCAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCttat                                                                                                    <  4:407743/148‑1 (MQ=255)
                                                  tatCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCttatttat                                                                                                <  8:85695/149‑1 (MQ=255)
                                                    tCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATg                                                                                               <  2:9647/148‑1 (MQ=255)
                                                        gCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGtt                                                                                            <  3:183381/147‑1 (MQ=255)
                                                        gCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTg                                                                                          <  2:347951/149‑1 (MQ=255)
                                                                 cATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGttt                                                                                           >  2:445402/1‑139 (MQ=255)
                                                                 cATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGttt                                                                                           <  1:445402/139‑1 (MQ=255)
                                                                   ttCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTAtt                                                                                                                                        <  8:406247/92‑1 (MQ=255)
                                                                       tATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaa                                                                           <  3:133516/149‑1 (MQ=255)
                                                                         gtGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaat                                                                          <  8:331773/147‑1 (MQ=255)
                                                                           gACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGGCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaatgaa                                                                       >  4:109373/1‑149 (MQ=255)
                                                                            ccGCCATCGTTTTGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGGCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTg                                                                                          <  6:154306/128‑1 (MQ=255)
                                                                            aCGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGCCCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaatgaat                                                                      >  7:208481/1‑149 (MQ=255)
                                                                            aCGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaatgaa                                                                       <  3:299349/148‑1 (MQ=255)
                                                                            aCGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTATTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGCCCTCCATAGCTTATTTATGGTTTg                                                                                          >  5:154306/1‑129 (MQ=255)
                                                                              gCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAAtt                                                                     >  3:133722/1‑148 (MQ=255)
                                                                                   cGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTc                                                               <  5:214011/149‑1 (MQ=255)
                                                                                      tATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTAc                                                                  <  1:412051/143‑1 (MQ=255)
                                                                                      tATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTAc                                                                  >  2:412051/1‑143 (MQ=255)
                                                                                           aaGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTActg                                                       >  1:424213/1‑149 (MQ=255)
                                                                                           aaGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTActg                                                       >  5:449447/1‑149 (MQ=255)
                                                                                           aaGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTActg                                                       >  4:155846/1‑149 (MQ=255)
                                                                                              gCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTActgctg                                                    <  5:168211/149‑1 (MQ=255)
                                                                                               cATCTGGCAAATGACCACACTCGGCGAATCGATTACTGGGGAGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGAGACCTCCATATCTTATTTATGGTTGGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGa                                                   <  8:177786/149‑1 (MQ=255)
                                                                                                aTCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGAt                                                  >  4:267041/1‑149 (MQ=255)
                                                                                                 tCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaa                                                                           >  6:12605/1‑123 (MQ=255)
                                                                                                 tCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGgaa                                                                           <  5:12605/123‑1 (MQ=255)
                                                                                                   tGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTgg                                                                              <  2:9637/118‑1 (MQ=255)
                                                                                                   tGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTgg                                                                              >  1:9637/1‑118 (MQ=255)
                                                                                                   tGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCCTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGc                                                >  4:11769/1‑148 (MQ=255)
                                                                                                     gCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTaa                                                                                      <  4:218456/108‑1 (MQ=255)
                                                                                                     gCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTaa                                                                                      >  3:218456/1‑108 (MQ=255)
                                                                                                     gCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCtt                                              <  8:208481/148‑1 (MQ=255)
                                                                                                     gCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCCTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGc                                                <  3:11769/146‑1 (MQ=255)
                                                                                                      cAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGa                                            <  7:347380/149‑1 (MQ=255)
                                                                                                       aaaTGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAg                                           >  2:314938/1‑149 (MQ=255)
                                                                                                       aaaTGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCCTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGa                                                   <  4:229805/141‑1 (MQ=255)
                                                                                                       aaaTGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCCTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGa                                                   >  3:229805/1‑141 (MQ=255)
                                                                                                              cAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACtt                                     >  7:177944/1‑148 (MQ=255)
                                                                                                              cAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACtt                                     >  8:184624/1‑148 (MQ=255)
                                                                                                               aaaCAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTa                                    >  8:97798/1‑148 (MQ=255)
                                                                                                                aaCAGGGCGATGCGATAACGGGGTCGATGGGCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATCGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTaaa                                  <  6:49337/149‑1 (MQ=255)
                                                                                                                     ggCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTg                             <  7:97798/149‑1 (MQ=255)
                                                                                                                     ggCGATGCGATAACGGGGTCGATGGTCTTTTTATGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTg                             <  3:109373/149‑1 (MQ=255)
                                                                                                                      gcgatgcgatAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTgg                            <  7:184624/149‑1 (MQ=255)
                                                                                                                      gcgatgcgatAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTgg                            <  1:314938/149‑1 (MQ=255)
                                                                                                                             gatAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCACCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGTTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGc                     >  3:373409/1‑149 (MQ=255)
                                                                                                                               tAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTc                   >  8:442122/1‑149 (MQ=255)
                                                                                                                                  cGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGa                 >  7:196067/1‑148 (MQ=255)
                                                                                                                                        cGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGATATCTAc          <  4:133722/149‑1 (MQ=255)
                                                                                                                                                tttttCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCAAATATATGGTTTGGTAACCGACTGGGAATGAATTTACGTCAGGTACTGCTGATGATTGAGCGACTTAAAGAGTCGGTAATGCTCGATATCTCCCTGGTCgg  >  6:293684/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTCTGCGTTATTGGTGCCCCGCTCATTCTGGTGACCTCCATAGCTTATTTATGGTTTGGTAACCGACTGGGAATGAATTTACGTCCGGTACTGCTGATGCTTGAGCGACTTAAAGAGTGGGTAATGCTCGATATCTACCTGGTCGG  >  NC_000913/1916477‑1916769

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