Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 349,615 T→C intergenic (+43/‑397) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA349,6150TC99.1% 322.0 / ‑5.9 107intergenic (+43/‑397)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (55/50);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAC  >  NC_000913_3_pae_tpiA/349468‑349720
                                                                                                                                                   |                                                                                                         
gTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCg                                                                                                         <  1:154115/150‑1 (MQ=255)
 tCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATACACAAGGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGcgccgc                                                                                                        >  3:329765/1‑150 (MQ=255)
 tCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGcgccgc                                                                                                        <  1:938254/150‑1 (MQ=255)
 tCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGcgccgc                                                                                                        >  2:551857/1‑150 (MQ=255)
 tCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGcgccgc                                                                                                        >  4:32216/1‑150 (MQ=255)
    tCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGCTGCCTGATGCGCCGCTGa                                                                                                     >  2:1309316/1‑150 (MQ=255)
      gACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcg                                                                                                   <  3:4503/150‑1 (MQ=255)
      gACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcg                                                                                                   >  3:660911/1‑150 (MQ=255)
      gACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcg                                                                                                   <  1:973432/150‑1 (MQ=255)
         acCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGGTGACGTGa                                                                                                >  3:154573/1‑150 (MQ=255)
         acCATGCAGACCCGCAACGAGCTGTACGAAAACATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcgtga                                                                                                >  4:1236117/1‑147 (MQ=255)
           cATGCAGACCCGCAACGAGCTGTACGAGAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACt                                                                                              >  1:1485404/1‑150 (MQ=255)
           cATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACt                                                                                              >  4:789052/1‑150 (MQ=255)
               cAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATc                                                                                          <  2:1477382/150‑1 (MQ=255)
               cAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATc                                                                                          <  1:1628510/150‑1 (MQ=255)
                aGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                         <  3:51409/150‑1 (MQ=255)
                aGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                         <  2:215890/150‑1 (MQ=255)
                aGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                         >  2:792858/1‑150 (MQ=255)
                aGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                         >  1:1828687/1‑150 (MQ=255)
                 gACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAg                                                                                        >  3:976757/1‑150 (MQ=255)
                   cccGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGc                                                                                      <  3:1138848/150‑1 (MQ=255)
                    ccGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGcc                                                                                     >  2:1389624/1‑150 (MQ=255)
                     cGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                           <  2:89613/143‑1 (MQ=255)
                     cGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                           >  1:89613/1‑143 (MQ=255)
                        aaCGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                 >  2:1876528/1‑150 (MQ=255)
                        aaCGAGCTGTACGAAAACATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                 <  3:1236117/150‑1 (MQ=255)
                         aCGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACgg                                                                                <  2:1828687/150‑2 (MQ=255)
                          cGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                               >  2:1334146/1‑148 (MQ=255)
                           gAGCTGTACCAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACgggg                                                                              >  1:1408970/1‑147 (MQ=255)
                               tGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                          >  3:408866/1‑150 (MQ=255)
                                gTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAc                                                                         <  1:952751/150‑1 (MQ=255)
                                 tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                        >  1:680336/1‑150 (MQ=255)
                                 tACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                        >  4:65551/1‑150 (MQ=255)
                                    gAAAGCATCAACTACTACCAGTACGGAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                     >  3:1244184/1‑150 (MQ=255)
                                    gAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                     <  2:1474117/150‑1 (MQ=255)
                                     aaaGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGGACTGaaa                                                                    >  1:1902967/1‑149 (MQ=38)
                                     aaaGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                    <  3:523907/150‑1 (MQ=255)
                                        gCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 <  2:125181/150‑1 (MQ=255)
                                        gCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGACGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 >  1:513325/1‑150 (MQ=255)
                                          aTCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                 <  1:355408/132‑1 (MQ=255)
                                          aTCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                 >  2:355408/1‑132 (MQ=255)
                                              actactACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGTACTGAACTGTAGGTCg                                                           >  2:1450508/1‑150 (MQ=37)
                                              actactACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 <  2:88874/144‑1 (MQ=255)
                                              actactACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 >  1:88874/1‑144 (MQ=255)
                                              actactACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                           <  1:1389624/150‑1 (MQ=255)
                                                tactacCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                         <  4:159135/150‑1 (MQ=255)
                                                 actacCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTTCGTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                        <  1:475541/150‑1 (MQ=37)
                                                 actacCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACtt                                                                                             <  4:167531/113‑1 (MQ=255)
                                                   tacCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAtgc                                                      >  2:218172/1‑148 (MQ=255)
                                                     ccAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                    <  4:660911/150‑1 (MQ=255)
                                                      cAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                   >  1:724719/1‑150 (MQ=255)
                                                       aGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGGCGGATAAGAca                                                  >  3:466590/1‑149 (MQ=25)
                                                       aGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                  <  2:569611/150‑1 (MQ=255)
                                                        gTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGACGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACgc                                                 <  1:1904084/150‑2 (MQ=255)
                                                           cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAta                                              <  1:1211892/150‑2 (MQ=255)
                                                            gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                <  2:1888214/131‑1 (MQ=255)
                                                            gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                >  1:1888214/1‑131 (MQ=255)
                                                            gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGACCTACGGGGTGCACTGAACTGTAGGTc                                                            >  3:1177657/1‑135 (MQ=25)
                                                            gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGACCTACGGGGTGCACTGAACTGTAGGTc                                                            <  4:1177657/135‑1 (MQ=25)
                                                             aaGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACgcatagc                                            <  1:1846234/150‑7 (MQ=21)
                                                               gagaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                          <  3:230143/150‑1 (MQ=21)
                                                               gagaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                          <  2:724719/150‑1 (MQ=21)
                                                               gagaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                          <  2:1485404/150‑1 (MQ=21)
                                                                agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                         <  1:865961/150‑1 (MQ=21)
                                                                   aGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                      <  1:1526287/150‑1 (MQ=18)
                                                                   aGCTCGACAACCTGTTTGCCCGTAGCCAGGCGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                      <  1:1334146/150‑1 (MQ=18)
                                                                    gCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                     >  3:796805/1‑150 (MQ=18)
                                                                      tCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                   >  4:537711/1‑150 (MQ=18)
                                                                       cGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                         <  1:1472837/95‑1 (MQ=255)
                                                                        gACAACCTTTTTGCCCGTAGCCAGGTTAAATAAAATACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACTGGGTGCACTGAACTGTAGGTCGTATAAGACGCATAGCGTCGCATCCGa                                 >  2:729492/1‑150 (MQ=17)
                                                                        gACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                 >  2:937680/1‑150 (MQ=18)
                                                                        gACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                 >  4:698233/1‑150 (MQ=18)
                                                                             ccTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                            >  1:376244/1‑150 (MQ=17)
                                                                                gTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                         <  4:121010/150‑1 (MQ=17)
                                                                                gTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                         >  3:998341/1‑150 (MQ=17)
                                                                                gTTTGCCCGTAGCCAGGTGAAACAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                         >  4:753163/1‑150 (MQ=17)
                                                                                 tttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                        >  2:514263/1‑150 (MQ=17)
                                                                                 tttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                        <  3:698233/150‑1 (MQ=17)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                  >  4:1104731/1‑121 (MQ=255)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                  <  3:1104731/121‑1 (MQ=255)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                     <  3:65551/150‑1 (MQ=17)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                     >  4:1218426/1‑150 (MQ=17)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                     <  4:56164/150‑1 (MQ=17)
                                                                                    gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                     <  1:885067/150‑1 (MQ=17)
                                                                                     cccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                    >  4:608649/1‑150 (MQ=17)
                                                                                     cccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                    <  3:608649/150‑1 (MQ=17)
                                                                                     cccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                    >  2:266202/1‑150 (MQ=17)
                                                                                      ccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                   <  1:858031/150‑1 (MQ=17)
                                                                                      ccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                   <  4:154573/150‑1 (MQ=17)
                                                                                       cGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                  <  1:1959390/150‑1 (MQ=17)
                                                                                           gCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                             >  3:35670/1‑135 (MQ=14)
                                                                                           gCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc              <  1:980622/150‑1 (MQ=14)
                                                                                             cAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGTACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg            <  1:1450508/150‑1 (MQ=14)
                                                                                                   aaaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcgtg                                                                                                 <  3:640259/59‑3 (MQ=25)
                                                                                                   aaaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGAcgtg                                                                                                 >  4:640259/1‑57 (MQ=25)
                                                                                                     aTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGACCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgcg    >  1:504539/1‑148 (MQ=11)
                                                                                                       aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                        >  3:1086854/1‑80 (MQ=17)
                                                                                                       aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 <  4:1051260/87‑1 (MQ=17)
                                                                                                       aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 >  3:1051260/1‑87 (MQ=17)
                                                                                                       aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc  >  3:502079/1‑150 (MQ=11)
                                                                                                                                                   |                                                                                                         
GTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAC  >  NC_000913_3_pae_tpiA/349468‑349720

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