Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 339,756 T→C 8.2% intergenic (+13/‑409) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,7560TC 8.2% 314.7 / 13.7 97intergenic (+13/‑409)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base C (4/4);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 101 total aligned reads displayed.

ATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATC  >  NC_000913/339609‑339899
                                                                                                                                                   |                                                                                                                                               
aTTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGTATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATg                                                                                                                                                <  4:835548/149‑1 (MQ=255)
 ttCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGc                                                                                                                                               >  2:1404016/1‑149 (MQ=255)
          tGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcgg                                                                                                                                      >  4:1263640/1‑149 (MQ=255)
           gATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggc                                                                                                                                     >  3:910795/1‑149 (MQ=255)
              cAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATccggcatg                                                                                                                                  >  1:19425/1‑146 (MQ=255)
               aaCACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATAAGGCGttc                                                                                                                                 >  4:475202/1‑147 (MQ=255)
                acacCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                                                                                                >  3:916092/1‑149 (MQ=255)
                acacCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGaa                                                                                                                                <  2:19425/149‑1 (MQ=255)
                  acCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGcc                                                                                                                                              <  3:1036210/133‑1 (MQ=255)
                  acCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGcc                                                                                                                                              >  4:1036210/1‑133 (MQ=255)
                  acCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                                                                                                >  2:1273461/1‑147 (MQ=255)
                  acCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                                                                                                <  1:1273461/147‑1 (MQ=255)
                   cccccATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGAAcat                                                                                                                             <  3:969460/146‑3 (MQ=255)
                    cGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCg                                                                                                                                             >  2:1299424/1‑132 (MQ=255)
                    cGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCg                                                                                                                                             <  1:1299424/132‑1 (MQ=255)
                    cGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                            >  2:23125/1‑149 (MQ=255)
                     gCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                           <  2:880115/149‑1 (MQ=255)
                        aTCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                        <  1:23125/149‑1 (MQ=255)
                         tCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                       <  1:1404016/149‑1 (MQ=255)
                            cccATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                    >  1:1181646/1‑149 (MQ=255)
                             ccATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                   <  3:475202/149‑1 (MQ=39)
                               aTAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                       >  2:419556/1‑143 (MQ=255)
                               aTAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                       <  1:419556/143‑1 (MQ=255)
                                 aaaGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcg                                                                                                                                       <  2:1414383/125‑1 (MQ=255)
                                 aaaGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcg                                                                                                                                       >  1:1414383/1‑125 (MQ=255)
                                          ggCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcgc                                                                                                      <  2:1181646/149‑1 (MQ=255)
                                           gCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGt                                                                                                                                   >  3:1443933/1‑119 (MQ=255)
                                           gCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGt                                                                                                                                   <  4:1443933/119‑1 (MQ=255)
                                                 gCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGCTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATc                                                                                               >  1:774168/1‑149 (MQ=255)
                                                    tGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGt                                                                                            <  3:1263640/149‑1 (MQ=255)
                                                      aTTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATccggcat                                                                                                                                   <  2:1235554/108‑3 (MQ=255)
                                                      aTTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATccggcat                                                                                                                                   >  1:1235554/1‑106 (MQ=255)
                                                          gcgcGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                      <  4:916092/149‑1 (MQ=255)
                                                               gtcATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTCCGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAATGGCGCAAGAATCTGTAGGCCTGATaa                                                                                 >  3:1146827/3‑149 (MQ=255)
                                                               ggCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATaa                                                                                 <  4:910795/149‑1 (MQ=255)
                                                                   ttGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                      <  1:575485/140‑1 (MQ=255)
                                                                   ttGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                      >  2:575485/1‑140 (MQ=255)
                                                                     gtgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAg                                                                           <  4:1146827/149‑1 (MQ=255)
                                                                         aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                            <  3:761241/96‑1 (MQ=255)
                                                                         aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                            >  4:761241/1‑96 (MQ=255)
                                                                         aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                                                          <  2:1214841/146‑1 (MQ=255)
                                                                         aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                                                          >  1:1214841/1‑146 (MQ=255)
                                                                            cACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                        <  1:1095074/129‑1 (MQ=255)
                                                                            cACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                        >  2:1095074/1‑129 (MQ=255)
                                                                             aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                            >  4:612890/1‑92 (MQ=255)
                                                                             aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                            <  3:612890/92‑1 (MQ=255)
                                                                             aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg                                                                                          <  4:506662/126‑1 (MQ=255)
                                                                             aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg                                                                                          >  3:506662/1‑126 (MQ=255)
                                                                              ccTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAg                                                                           >  4:1431654/1‑140 (MQ=255)
                                                                              ccTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAg                                                                           <  3:1431654/140‑1 (MQ=255)
                                                                               cTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                               >  2:785034/1‑103 (MQ=255)
                                                                               cTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                               <  1:785034/103‑1 (MQ=255)
                                                                                ttttGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggc                                                                                                                                     >  3:826077/1‑80 (MQ=255)
                                                                                ttttGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggc                                                                                                                                     <  4:826077/80‑1 (MQ=255)
                                                                                   tGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                             <  2:774168/149‑1 (MQ=255)
                                                                                       tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAAc                                                                                                                               <  1:1004700/79‑1 (MQ=255)
                                                                                       tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAAc                                                                                                                               >  2:1004700/1‑79 (MQ=255)
                                                                                       tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATg                                                                                                                                  >  1:344423/1‑76 (MQ=38)
                                                                                       tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATg                                                                                                                                  <  2:344423/76‑1 (MQ=38)
                                                                                            cGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAtt                                                    <  4:1214444/149‑1 (MQ=255)
                                                                                            cGAGAAAGCCATTCTTGGCCGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAtt                                                    >  3:1214444/1‑149 (MQ=255)
                                                                                              agaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                     >  1:766890/1‑82 (MQ=255)
                                                                                              agaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                     <  2:766890/82‑1 (MQ=255)
                                                                                              agaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                                                           >  4:259077/1‑92 (MQ=255)
                                                                                              agaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                                                           <  3:259077/92‑1 (MQ=255)
                                                                                                  aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTg                                                                                     <  3:638600/110‑1 (MQ=255)
                                                                                                  aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTg                                                                                     >  4:638600/1‑110 (MQ=255)
                                                                                                   gCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    <  1:717439/110‑1 (MQ=255)
                                                                                                   gCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    <  1:1039172/110‑1 (MQ=255)
                                                                                                   gCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    >  2:717439/1‑110 (MQ=255)
                                                                                                   gCCATTATTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    >  2:1039172/1‑110 (MQ=255)
                                                                                                    ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg                                                                                          <  4:534554/103‑1 (MQ=255)
                                                                                                    ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg                                                                                          >  3:534554/1‑103 (MQ=255)
                                                                                                       ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                                                             >  3:552602/1‑97 (MQ=255)
                                                                                                       ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                                                             <  4:552602/97‑1 (MQ=255)
                                                                                                        tCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCAt                                                                                                                                                 >  4:186474/1‑44 (MQ=255)
                                                                                                        tCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCAt                                                                                                                                                 <  3:186474/44‑1 (MQ=255)
                                                                                                               tGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                           <  2:1087855/91‑1 (MQ=255)
                                                                                                               tGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                           >  1:1087855/1‑91 (MQ=255)
                                                                                                                   gCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGc                                <  4:1122786/146‑1 (MQ=255)
                                                                                                                   gCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGc                                >  3:1122786/1‑146 (MQ=255)
                                                                                                                         gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                             >  4:68354/1‑47 (MQ=255)
                                                                                                                         gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                             <  3:68354/47‑1 (MQ=255)
                                                                                                                                 gTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGaa               <  3:608578/149‑1 (MQ=255)
                                                                                                                                 gTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGaa               >  4:608578/1‑149 (MQ=255)
                                                                                                                                        cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                               <  1:165120/94‑1 (MQ=255)
                                                                                                                                        cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                               >  2:165120/1‑94 (MQ=255)
                                                                                                                                        cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGc                                                  >  3:156697/1‑107 (MQ=255)
                                                                                                                                        cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGc                                                  <  4:156697/107‑1 (MQ=255)
                                                                                                                                            aTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATaa                                                                                 <  3:574614/72‑1 (MQ=40)
                                                                                                                                            aTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCGGTAGGCCTGATaa                                                                                 >  4:574614/1‑72 (MQ=255)
                                                                                                                                              gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATc  >  2:1347928/1‑149 (MQ=255)
                                                                                                                                               ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    <  1:190158/66‑1 (MQ=37)
                                                                                                                                               ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    >  2:190158/1‑66 (MQ=37)
                                                                                                                                                tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                   >  2:306720/1‑82 (MQ=33)
                                                                                                                                                tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                   <  1:306720/82‑1 (MQ=33)
                                                                                                                                                  aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    <  4:540243/63‑1 (MQ=31)
                                                                                                                                                  aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                    >  3:540243/1‑63 (MQ=31)
                                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAt                                                     >  2:290019/1‑93 (MQ=34)
                                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAt                                                     <  1:290019/93‑1 (MQ=34)
                                                                                                                                                   |                                                                                                                                               
ATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATC  >  NC_000913/339609‑339899

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