Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,973,862 T→G 5.6% G85G (GGT→GGG rffH → glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,973,8620TG 5.6% 227.6 / 2.3 92G85G (GGT→GGGrffHglucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (31/55);  new base G (2/4);  total (33/59)
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.86e-01

TATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTAT  >  NC_000913/3973713‑3974012
                                                                                                                                                     |                                                                                                                                                      
tATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTACAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGt                                                                                                                                                        >  1:1088220/1‑150 (MQ=255)
   ccGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTg                                                                                                                                                     <  2:415258/150‑1 (MQ=255)
   ccGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTg                                                                                                                                                     >  1:27449/1‑150 (MQ=255)
    cGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTgg                                                                                                                                                    <  2:495217/150‑1 (MQ=255)
      cTGGCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAACGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGACCGTCCGGGCTgggc                                                                                                                                                  <  2:1120180/150‑3 (MQ=255)
        gTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCa                                                                                                                                                >  1:861965/1‑150 (MQ=255)
          cGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAgg                                                                                                                                              >  2:323978/1‑150 (MQ=255)
          cGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAgg                                                                                                                                              <  2:771169/150‑1 (MQ=255)
           ggTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGc                                                                                                                                             >  1:939704/1‑150 (MQ=255)
            gTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGcc                                                                                                                                            <  1:1207465/150‑1 (MQ=255)
              gCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGACCGGACGGTCTGGCGCAGGACtt                                                                                                                                          >  1:15440/1‑150 (MQ=255)
               cTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCttt                                                                                                                                         >  2:298905/1‑150 (MQ=255)
                tGATGCTGGACGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGcggtta                                                                                                                                        >  1:656933/1‑145 (MQ=255)
                   ggCTGGCCGGTAGCCGCGAAATTCTCATCCTCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGCGTGCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTatca                                                                                                                                     <  1:94247/149‑1 (MQ=255)
                     taggCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAACGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGGCTGGCGCAGGCCTTTatcatc                                                                                                                                   <  2:807169/148‑1 (MQ=255)
                       ggCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGTTAACGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGTGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCgg                                                                                                                                 <  1:11929/150‑1 (MQ=255)
                       ggCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGTCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCgg                                                                                                                                 <  2:763902/150‑1 (MQ=255)
                       ggCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCgg                                                                                                                                 <  1:1251991/150‑1 (MQ=255)
                       ggCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCGTTATCATCgg                                                                                                                                 >  2:792534/1‑150 (MQ=255)
                        gCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGt                                                                                                                                <  1:1305164/150‑1 (MQ=255)
                         ccGGTATCCACGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTg                                                                                                                               <  2:744228/150‑1 (MQ=255)
                           ggTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGGTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACAGTCTGGCGCAGGCCTTTATCATCGGTGaa                                                                                                                             <  2:885119/150‑1 (MQ=255)
                                  gcgAAATTCTCATCATCACTACGCCGGAAGATTAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAAccttcc                                                                                                                      <  2:342077/150‑1 (MQ=255)
                                  gcgAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATTGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAAccttcc                                                                                                                      <  2:711700/150‑1 (MQ=255)
                                    gAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAAccttcctt                                                                                                                    <  1:834289/150‑1 (MQ=255)
                                    gAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAAccttcctt                                                                                                                    <  1:984906/150‑1 (MQ=255)
                                      aaTTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTaa                                                                                                                  <  1:525215/150‑1 (MQ=255)
                                       aTTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAAt                                                                                                                 >  2:1316280/1‑150 (MQ=255)
                                           tcatcatcaCTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTg                                                                                                             <  1:869872/150‑1 (MQ=255)
                                            catcatcaCTACGCCGGAAGATAAAGGATATTTCCAGCGCCTGCTTCGCGATGGTAGTGAGTTCTGTATCCAGCAGCAATATGCCGAACAGCCCAGCCCGGACAGTCTGCCGCAGGCCATTATCAACTGTGAAACCTTCCTTAATGGTGa                                                                                                            >  2:77886/1‑150 (MQ=255)
                                               catcaCTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAAcc                                                                                                         <  2:80388/150‑1 (MQ=255)
                                               aatcaCTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAAcc                                                                                                         <  2:80386/149‑1 (MQ=255)
                                                  cccTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTAAAACCTTCCTTAATGGTGAACcttc                                                                                                      <  2:329644/148‑1 (MQ=255)
                                                    ctccGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGAGGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACcttctt                                                                                                    <  2:984631/147‑1 (MQ=255)
                                                     tACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTg                                                                                                   <  1:247498/150‑1 (MQ=255)
                                                         ccGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTg                                                                                               >  1:1305710/1‑150 (MQ=255)
                                                         ccGGAAGATAAAGGTTATTTCCAGCGCCGGCTGGGCGATGGTAGTGAGGTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTg                                                                                               <  2:42932/150‑1 (MQ=255)
                                                          cGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTGGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTgg                                                                                              >  2:730875/1‑150 (MQ=255)
                                                                aTAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTgg                                                                                        >  1:944250/1‑150 (MQ=255)
                                                                 tAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAACCCTTCCTTAATGGTGATCCTTCTTGTCTGGTGCTggg                                                                                       >  2:1225037/1‑150 (MQ=255)
                                                                 aaaaGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTggg                                                                                       <  2:778143/149‑1 (MQ=255)
                                                                  aaaGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGc                                                                                      <  2:629966/150‑1 (MQ=255)
                                                                  aaaGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGc                                                                                      <  2:902428/150‑1 (MQ=255)
                                                                  aaaGGTTATTTCCAGCGCCTGCTGGGCGAGGGGAGGGAGTTCGGTATCCAGCTGGTATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGc                                                                                      <  2:1303151/150‑1 (MQ=255)
                                                                  aaaGGTTATATCCAGCGCCTGCTGGGCGATGGGAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGc                                                                                      <  2:80008/150‑1 (MQ=255)
                                                                     ggTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGAt                                                                                   >  1:967662/1‑150 (MQ=255)
                                                                      gTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATa                                                                                  <  1:162994/150‑1 (MQ=255)
                                                                          tttCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGCCGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAt                                                                              >  2:497408/1‑150 (MQ=255)
                                                                           ttCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtc                                                                             <  1:274721/150‑1 (MQ=255)
                                                                           ttCCAGCGCCGGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCGGGAATATGCCGAACAGCCCAGCCCGGACGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtc                                                                             <  1:247675/150‑1 (MQ=255)
                                                                            tCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGTGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtct                                                                            <  1:349564/150‑1 (MQ=255)
                                                                            tCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtct                                                                            >  2:297998/1‑150 (MQ=255)
                                                                            tCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtct                                                                            >  1:1247282/1‑150 (MQ=255)
                                                                             ccAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtctt                                                                           >  2:147332/1‑150 (MQ=255)
                                                                               aGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACAtcttct                                                                         <  1:1309062/150‑1 (MQ=255)
                                                                                gcgcCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATcttctt                                                                        <  1:496910/150‑1 (MQ=255)
                                                                                 ggcCTGCTGGGCGTTGGTAGTGATTTCGGTATCCATCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTAGTGCTGGGCGATAACATCttcttc                                                                       <  2:656933/149‑1 (MQ=255)
                                                                                  gcCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCg                                                                      <  2:1225621/150‑1 (MQ=255)
                                                                                  gcCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCg                                                                      >  1:385206/1‑150 (MQ=255)
                                                                                  gcCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCg                                                                      <  2:1151647/150‑1 (MQ=255)
                                                                                         gggCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAggg                                                               >  2:357520/1‑150 (MQ=255)
                                                                                         gggCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCGTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAggg                                                               <  1:970395/150‑1 (MQ=255)
                                                                                          ggCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGc                                                              <  1:520303/150‑1 (MQ=255)
                                                                                           gCGATGGTAGTGAGTTCGGTATCCAGCGGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCt                                                             <  2:257624/150‑1 (MQ=255)
                                                                                               tGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAg                                                         >  2:649992/1‑150 (MQ=255)
                                                                                                     tGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCATAATGGTGAACCATCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGaa                                                   >  2:328853/1‑150 (MQ=255)
                                                                                                       tgTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGc                                                 <  2:969851/149‑1 (MQ=255)
                                                                                                        gTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCt                                                >  1:683580/1‑150 (MQ=255)
                                                                                                         ttCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTg                                               >  1:1141041/1‑150 (MQ=255)
                                                                                                            gggaTCCTGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGt                                            <  1:822225/147‑1 (MQ=255)
                                                                                                            ggTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGt                                            >  2:159854/1‑150 (MQ=255)
                                                                                                                      tGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTgcgg                                  <  2:119640/150‑1 (MQ=255)
                                                                                                                      tGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTgcgg                                  <  2:126702/150‑1 (MQ=255)
                                                                                                                      tGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTgcgg                                  <  2:1271229/150‑1 (MQ=255)
                                                                                                                       ggAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATCGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTgcggc                                 >  1:232426/1‑150 (MQ=255)
                                                                                                                        gAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTAGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTgcggcg                                <  2:735971/150‑1 (MQ=255)
                                                                                                                            atGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCa                            <  1:459505/150‑1 (MQ=255)
                                                                                                                               cccaaCAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGGGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCTTCATGTTTCGGCGCGCACCg                         <  2:402295/147‑1 (MQ=255)
                                                                                                                                  cctgcgccagaccGTCCGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGCAGCTGCGTCATGTTGCGGCGCGCACCGAAg                      >  2:374681/7‑150 (MQ=255)
                                                                                                                                  aaCAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAg                      <  2:445911/150‑1 (MQ=255)
                                                                                                                                   ctgcgccagaccGTCCGGGCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAgg                     >  1:963716/6‑150 (MQ=255)
                                                                                                                                   aCAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGGCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAgg                     <  1:694068/150‑1 (MQ=255)
                                                                                                                                     agcccagcccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAgggg                   <  2:936171/150‑1 (MQ=255)
                                                                                                                                      gccccgcccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGGCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGc                  <  1:485997/150‑1 (MQ=255)
                                                                                                                                      gcccagcccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGc                  <  2:23471/150‑1 (MQ=255)
                                                                                                                                         cagcccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTATCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGAc               >  1:24230/1‑150 (MQ=255)
                                                                                                                                          agcccTGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCTCACCGAAGGGGCGACg              <  2:1088220/150‑1 (MQ=255)
                                                                                                                                           gcccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACgg             >  2:1072808/1‑150 (MQ=255)
                                                                                                                                             ccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGtt           >  2:877503/1‑150 (MQ=255)
                                                                                                                                             ccGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGGCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGtt           <  1:605231/150‑1 (MQ=255)
                                                                                                                                                    gTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTACTGGTGAACCTTCTTGGCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCGGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCt    <  1:1258304/150‑1 (MQ=255)
                                                                                                                                                    gTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCt    <  2:27449/150‑1 (MQ=255)
                                                                                                                                                    gTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCt    >  1:666417/1‑150 (MQ=255)
                                                                                                                                                     tCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTa   <  1:649992/150‑1 (MQ=255)
                                                                                                                                                     tCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATACATCTTCTTCGGTCAGGGCTTCAGGCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTAt  >  2:1245064/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TATCCGCTGTCGGTGCTGATGCTGGCCGGTATCCGCGAAATTCTCATCATCACTACGCCGGAAGATAAAGGTTATTTCCAGCGCCTGCTGGGCGATGGTAGTGAGTTCGGTATCCAGCTGGAATATGCCGAACAGCCCAGCCCGGACGGTCTGGCGCAGGCCTTTATCATCGGTGAAACCTTCCTTAATGGTGAACCTTCTTGTCTGGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTAT  >  NC_000913/3973713‑3974012

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