Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,658 C→T 5.9% intergenic (‑560/+38) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,6580CT 5.9% 271.7 / 6.1 86intergenic (‑560/+38)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (55/25);  new base T (2/3);  total (57/29)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAG  >  NC_000913/4153512‑4153805
                                                                                                                                                  |                                                                                                                                                   
gCGGCACTGGAGTTTGGCAACAGTGCCGTATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGcc                                                                                                                                                  <  1:358103/150‑1 (MQ=255)
gCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGcc                                                                                                                                                  >  2:1330518/1‑150 (MQ=39)
        ggAGTTTGGCAACCGTGCCGGATGCGGCGCGAGCGCCTGAGCCGGCCTACGGTGGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATgcgg                                                                                                                                          <  2:71075/150‑1 (MQ=255)
          aGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCACTGTAGGTCGGAGAAGATGCGCACGCTGCGCATCCGACATTATTGGGGCACGGGAGTTTGGCAACAGCGCCGGATgcggcg                                                                                                                                        >  2:1377082/1‑150 (MQ=14)
             ttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCTGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTGTGGCAACAGCGCCGGATGCGGCGCGa                                                                                                                                     >  2:328062/1‑150 (MQ=255)
                gCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAgcg                                                                                                                                  >  1:1168253/1‑150 (MQ=39)
                      gagactgATGCGGCGCGCGCCCCTCATCCGGCCTACGGGTGGGCATCGTTCGATACTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                            <  2:351404/144‑1 (MQ=14)
                       tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                           >  2:696242/1‑150 (MQ=39)
                        gCCGGATGCGGCGCGAGCGCCATATACGGCCAACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAACCAACGCATCCGACATAAAAGCGGCACAGGAGTTTGGCAAAAGCGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                          >  1:860881/1‑150 (MQ=14)
                          cgcttGCGGCGCGAGCGCGAAATCCGGCCCACGGTGGGGCAGCGACCGCGCCCGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                        <  1:509444/146‑1 (MQ=11)
                            gATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGAATCGTTTGAGTCAGTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                      <  1:281328/150‑1 (MQ=255)
                            gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                      >  2:459326/1‑150 (MQ=39)
                            gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                      <  1:281299/150‑1 (MQ=39)
                            gATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGGTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                      >  2:429482/1‑150 (MQ=255)
                                 ggcgCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGt                                                                                                                 >  1:1385906/1‑150 (MQ=25)
                                  gcgcGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtt                                                                                                                >  2:1284825/1‑150 (MQ=38)
                                  gcgcGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtt                                                                                                                >  1:730358/1‑150 (MQ=40)
                                   cgcgAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGttt                                                                                                               >  2:36057/1‑150 (MQ=255)
                                    gcgAGCGCCTTATCCGGCCTACTGTTGGGCATCGTTTGAGTCCCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGtttt                                                                                                              <  2:1359867/150‑1 (MQ=255)
                                      gAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTaa                                                                                                            >  2:1148010/1‑150 (MQ=255)
                                       cgcgcCTTATCCGTCCTAATCGTGGGCAGATTTAGATTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAAt                                                                                                           <  2:1335225/149‑1 (MQ=14)
                                         cgcCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGc                                                                                                         >  1:1031826/1‑150 (MQ=255)
                                          gcCTTATCCGGCCTACGGTTGGACATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGcc                                                                                                        >  1:853340/1‑150 (MQ=255)
                                           ccTTATCCGGCCTACGGGTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCa                                                                                                       <  1:716314/150‑1 (MQ=255)
                                            cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAg                                                                                                      >  2:35299/1‑150 (MQ=255)
                                            cTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTAATCCGGCCTACGTTTTAATGCCAg                                                                                                      >  2:1294507/1‑150 (MQ=255)
                                             ttATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGc                                                                                                     >  1:871631/1‑150 (MQ=255)
                                                     ccatcGGGTGGGAATAGTGTGAGTCTCTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCTCTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATgg                                                                                             <  1:819379/146‑1 (MQ=17)
                                                      cTACTGTTGGGCATCGTTTGCGTCGGTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGt                                                                                            <  2:990261/150‑1 (MQ=17)
                                                           gTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAtt                                                                                       >  2:765229/1‑150 (MQ=255)
                                                              gggCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTAcc                                                                                    >  1:775234/1‑150 (MQ=255)
                                                                gCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAAATGTGAATTACCTg                                                                                  >  1:531088/1‑150 (MQ=255)
                                                                 cATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTgg                                                                                 >  1:781099/1‑150 (MQ=255)
                                                                      tttGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTAt                                                                            <  1:1168102/150‑1 (MQ=255)
                                                                          aGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGt                                                                        >  2:1149464/1‑150 (MQ=255)
                                                                                tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                                                                  >  1:1195820/1‑150 (MQ=255)
                                                                                tGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCgg                                                                  >  1:1210745/1‑150 (MQ=255)
                                                                                    ggTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTgg                                                              <  1:1279969/150‑1 (MQ=255)
                                                                                     gTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCGAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                                                             >  1:518194/1‑150 (MQ=255)
                                                                                        ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCtt                                                          <  2:682952/150‑1 (MQ=255)
                                                                                              gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                                                             >  2:606477/1‑141 (MQ=255)
                                                                                              gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTggg                                                             <  1:606477/141‑1 (MQ=255)
                                                                                               aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                   >  1:1092042/1‑150 (MQ=255)
                                                                                               aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                   >  1:1092044/1‑150 (MQ=255)
                                                                                               aTGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAAc                                                   >  1:1092045/1‑150 (MQ=255)
                                                                                                 gcgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCg                                                 >  1:549659/1‑150 (MQ=255)
                                                                                                  cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                >  1:386995/1‑150 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                >  2:501556/1‑150 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                >  2:501553/1‑150 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                <  1:849219/150‑1 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                <  1:288541/150‑1 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                >  1:195508/1‑150 (MQ=255)
                                                                                                  cgcAAGCATCGCATCCGACAGTCTTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgt                                                <  1:537456/150‑1 (MQ=255)
                                                                                                   gcAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCgtg                                               >  1:1167887/1‑150 (MQ=255)
                                                                                                       gcatcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGGTTc                                           >  1:1067020/1‑150 (MQ=255)
                                                                                                         atcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCa                                         <  2:386995/150‑1 (MQ=255)
                                                                                                          tcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAg                                        <  1:485964/150‑1 (MQ=255)
                                                                                                           ngcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGa                                       >  1:501892/2‑150 (MQ=255)
                                                                                                              atcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGatat                                    >  2:977449/1‑150 (MQ=255)
                                                                                                                ccGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtc                                  <  1:1284825/150‑1 (MQ=255)
                                                                                                                 cGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtca                                 >  1:797179/1‑150 (MQ=255)
                                                                                                                  gACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcat                                >  2:216420/1‑150 (MQ=255)
                                                                                                                   aCATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATtcatc                               >  1:1039295/1‑150 (MQ=255)
                                                                                                                       tattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                           >  1:1100817/1‑150 (MQ=255)
                                                                                                                       tattGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGt                           >  1:397191/1‑150 (MQ=255)
                                                                                                                         ttGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtg                         >  2:1251360/1‑150 (MQ=255)
                                                                                                                            cGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtgatg                      >  2:983445/1‑150 (MQ=255)
                                                                                                                             ggCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTtgatga                     >  2:76308/1‑150 (MQ=255)
                                                                                                                                    gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAt              <  1:1200178/150‑1 (MQ=255)
                                                                                                                                    gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAt              >  1:76365/1‑150 (MQ=255)
                                                                                                                                    gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAt              >  1:826322/1‑150 (MQ=255)
                                                                                                                                    gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAt              >  1:183174/1‑150 (MQ=255)
                                                                                                                                    gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGAt              <  1:1200166/150‑1 (MQ=255)
                                                                                                                                       tttGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTaa           <  1:696242/150‑1 (MQ=255)
                                                                                                                                           gCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGa       <  1:420940/150‑1 (MQ=255)
                                                                                                                                           gCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGa       <  2:797179/150‑1 (MQ=255)
                                                                                                                                            cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAg      <  1:844517/150‑1 (MQ=255)
                                                                                                                                            cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAg      >  2:1174516/1‑150 (MQ=255)
                                                                                                                                            cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAg      >  1:1093697/1‑150 (MQ=255)
                                                                                                                                            cAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCAGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCAACAGGTTGATGAGCCTGATTAATTTAg      >  1:1084476/1‑150 (MQ=255)
                                                                                                                                               tcgcgcCGCATCCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCa   >  1:1210803/3‑150 (MQ=255)
                                                                                                                                               cAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCa   <  2:421846/150‑1 (MQ=255)
                                                                                                                                               cAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCa   <  2:421655/150‑1 (MQ=255)
                                                                                                                                                aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATTGTGAATTAACAGGGTTATCAGTTAGCGGGTGGGCTTGATAAACCGGGTTTCCAGATATTCATAAGGTTGATGAGCCTGATTAATTTAGCCAg  >  1:533249/1‑150 (MQ=39)
                                                                                                                                                aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCCTCAGGTTGTTGAGCCTGATTAATTGAGCCAg  >  1:533227/1‑150 (MQ=255)
                                                                                                                                                aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAg  >  2:875693/1‑150 (MQ=255)
                                                                                                                                                aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAAGTACCTCGGTTATCAGTTCGCCGGTGGGCTTGATAAACCGTGTTTCCAGATAGTCATCAGGTTGATGAGCCTGATGAATGGAGCCAg  >  2:1079526/1‑150 (MQ=39)
                                                                                                                                                  |                                                                                                                                                   
GCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAG  >  NC_000913/4153512‑4153805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: