Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,217,943 T→C 7.4% Y279H (TAT→CAT)  aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,217,9430TC 7.4% 240.3 / 7.7 81Y279H (TAT→CAT) aceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (38/37);  new base C (2/4);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAA  >  NC_000913/4217798‑4218086
                                                                                                                                                 |                                                                                                                                               
tGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATgc                                                                                                                                             >  1:439045/1‑150 (MQ=255)
   ccGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGcc                                                                                                                                          <  2:1166413/150‑1 (MQ=255)
      tACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATACGCCagg                                                                                                                                       <  2:920868/150‑3 (MQ=255)
       aCCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATACGCCaggc                                                                                                                                      <  2:838927/150‑4 (MQ=255)
       aCCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATACGCCaggc                                                                                                                                      <  1:44897/150‑4 (MQ=255)
         cGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGc                                                                                                                                    <  2:1487363/150‑1 (MQ=255)
         cGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGc                                                                                                                                    <  2:87991/150‑1 (MQ=255)
          gATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCt                                                                                                                                   <  1:237453/150‑1 (MQ=255)
            tGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGa                                                                                                                                 >  1:641118/1‑150 (MQ=255)
                 aTGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctgg                                                                                                                            >  2:55292/1‑150 (MQ=255)
                        ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtg                                                                                                                     <  1:1252698/150‑1 (MQ=255)
                        ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTGCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGTCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtg                                                                                                                     <  2:1518306/150‑1 (MQ=255)
                         gATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtgt                                                                                                                    <  1:106506/150‑1 (MQ=255)
                            cTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGaa                                                                                                                 >  1:63670/1‑150 (MQ=255)
                            cTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGaa                                                                                                                 >  1:906612/1‑150 (MQ=255)
                                 cACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTc                                                                                                            <  1:360805/150‑1 (MQ=255)
                                  aCCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTcc                                                                                                           <  1:390076/150‑1 (MQ=255)
                                         aTTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCgg                                                                                                    <  2:456526/150‑1 (MQ=255)
                                           tGCGACCCGTATGACAGCGAGTTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGAt                                                                                                  >  1:384566/1‑150 (MQ=255)
                                           tGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGAt                                                                                                  >  1:122880/1‑150 (MQ=255)
                                            gCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATc                                                                                                 <  1:1218512/150‑1 (MQ=255)
                                              gACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTg                                                                                               >  2:1468432/1‑150 (MQ=255)
                                              gACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTg                                                                                               >  2:1473803/1‑150 (MQ=255)
                                                cccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGa                                                                                             >  2:551310/1‑150 (MQ=255)
                                                cccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGa                                                                                             <  1:424156/150‑1 (MQ=255)
                                                  cGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAAc                                                                                           <  1:1281735/150‑1 (MQ=255)
                                                      tGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgc                                                                                       >  2:1252487/1‑150 (MQ=255)
                                                      tGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgc                                                                                       >  1:1289974/1‑150 (MQ=255)
                                                       gACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcg                                                                                      <  1:162362/150‑1 (MQ=255)
                                                        aCAGCGAATTTATTACCGGCGAGCGTACCAGGGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                     <  2:140847/150‑1 (MQ=255)
                                                         cAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTACCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGcgcg                                                                                    <  1:138608/150‑1 (MQ=255)
                                                            cGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCg                                                                                 >  2:118983/1‑150 (MQ=255)
                                                               aTTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCtt                                                                              <  2:818676/150‑1 (MQ=255)
                                                                  tattaCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGc                                                                           <  2:606448/150‑1 (MQ=255)
                                                                     taCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGcaca                                                                        >  2:875306/1‑150 (MQ=255)
                                                                      aCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACaa                                                                       >  1:1318481/1‑150 (MQ=255)
                                                                       ccGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAg                                                                      >  1:147105/1‑150 (MQ=255)
                                                                        cGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGc                                                                     <  1:55292/150‑1 (MQ=255)
                                                                        cGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGc                                                                     <  1:1309482/150‑1 (MQ=255)
                                                                          gCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTa                                                                   <  2:854801/150‑1 (MQ=255)
                                                                             aGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATcc                                                                >  2:678772/1‑150 (MQ=255)
                                                                                gTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcg                                                             >  2:208868/1‑150 (MQ=255)
                                                                                 tACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcgc                                                            >  2:1200298/1‑150 (MQ=255)
                                                                                  aCCAGTGAAGGCTTCTTCCGTACTAATTCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATATGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcgcg                                                           >  2:1248846/1‑150 (MQ=255)
                                                                                     aGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGaaa                                                        <  1:1252487/150‑1 (MQ=255)
                                                                                        gAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAatat                                                     >  1:1041876/1‑150 (MQ=255)
                                                                                        gAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAatat                                                     >  1:534088/1‑150 (MQ=255)
                                                                                          aGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATcc                                                   >  2:112871/1‑150 (MQ=255)
                                                                                           ggCTTCTTCCGTACTCATGCTGGCATTGAGCAAGCGATCAGCCGTTGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCg                                                  <  2:1165800/150‑1 (MQ=255)
                                                                                             cttcttCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCggg                                                >  2:455938/1‑150 (MQ=255)
                                                                                               tcttcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCa                                              <  1:1443354/150‑1 (MQ=255)
                                                                                                  tcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAAc                                           <  1:1468432/150‑1 (MQ=255)
                                                                                                  tcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAAc                                           <  1:1473803/150‑1 (MQ=255)
                                                                                                       aCTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAActgctg                                      <  2:534088/150‑1 (MQ=255)
                                                                                                          cATGCGGGCATTGAGCAAGCGATCAGTCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                   >  2:496225/1‑150 (MQ=255)
                                                                                                          cATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                   >  2:950863/1‑150 (MQ=255)
                                                                                                          cATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                   >  2:1295726/1‑150 (MQ=255)
                                                                                                          cATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCt                                   >  1:881541/1‑150 (MQ=255)
                                                                                                           aTGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCtt                                  <  2:137750/150‑1 (MQ=255)
                                                                                                            tGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTta                                 <  1:1056717/150‑1 (MQ=255)
                                                                                                             gCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTtat                                <  2:1279712/150‑1 (MQ=255)
                                                                                                                ggCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAAc                             <  1:187369/150‑1 (MQ=255)
                                                                                                                 gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                            >  1:613953/1‑150 (MQ=255)
                                                                                                                 gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                            >  1:838520/1‑150 (MQ=255)
                                                                                                                 gCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                            >  1:1204309/1‑150 (MQ=255)
                                                                                                                 atcctGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACt                            <  1:548343/146‑1 (MQ=255)
                                                                                                                  cATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTg                           >  2:933323/1‑150 (MQ=255)
                                                                                                                   aTTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGCCCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGc                          >  2:932441/1‑150 (MQ=255)
                                                                                                                     tGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTc                        >  2:401130/1‑150 (MQ=255)
                                                                                                                       tgCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGc                      <  2:1299280/149‑1 (MQ=255)
                                                                                                                       aGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGc                      >  1:1515413/1‑150 (MQ=255)
                                                                                                                       aGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGc                      >  1:1510275/1‑150 (MQ=255)
                                                                                                                              gATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCTGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGtt               <  2:1204309/150‑1 (MQ=255)
                                                                                                                               aTCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTc              <  1:416472/150‑1 (MQ=255)
                                                                                                                                     cGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTgg        >  2:162990/1‑150 (MQ=255)
                                                                                                                                      gTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGc       <  2:439045/150‑1 (MQ=255)
                                                                                                                                       gcatatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCa                                              <  2:129903/105‑1 (MQ=255)
                                                                                                                                       gcatatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCa                                              >  1:129903/6‑110 (MQ=255)
                                                                                                                                       gcatatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCa      >  2:217832/6‑150 (MQ=255)
                                                                                                                                        ggCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAg     <  1:208868/150‑1 (MQ=255)
                                                                                                                                           cTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTAGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGaaa  >  1:1001207/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TGCCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAA  >  NC_000913/4217798‑4218086

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