Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,207,038 G→A 5.3% P192P (CCG→CCA ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,207,0380GA 5.3% 302.9 / 3.1 93P192P (CCG→CCAttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (40/48);  new base A (2/3);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

TGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACTCAACCGTCTG  >  NC_000913/3206894‑3207185
                                                                                                                                                |                                                                                                                                                   
tGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACt                                                                                                                                                <  2:408527/150‑1 (MQ=255)
tGCGGAAATCGAAGTTTACATGGCAGACGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCAGTAcc                                                                                                                                                <  2:1600515/150‑2 (MQ=255)
     aaaTCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACtggtgg                                                                                                                                           <  1:753027/150‑1 (MQ=255)
     aaaTCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTAACTGGCCGCTCGAAAGTGTTAATGACGTCAGAAGGCTACGAAGGCGTTGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACtggtgg                                                                                                                                           <  2:630843/150‑1 (MQ=255)
      aaTCGAAGTTTACATGGCAGGTGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGgtggg                                                                                                                                          <  1:1351376/150‑1 (MQ=255)
          gAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCAGTACcggtggacagg                                                                                                                                      >  1:889950/1‑145 (MQ=255)
          gAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCAGTACcggtggacagg                                                                                                                                      <  2:889950/150‑6 (MQ=255)
          caaGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCACCACCCTCGCCGTAAACGCCTGCCCACCAGTACcggtggacagg                                                                                                                                      <  1:1584852/149‑6 (MQ=255)
           aaGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGgtgggcgtgg                                                                                                                                     <  1:181395/150‑1 (MQ=255)
           aaGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGgtgggcgtgg                                                                                                                                     <  2:1305168/150‑1 (MQ=255)
            aGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGgtgggcgtggg                                                                                                                                    <  1:1313967/150‑1 (MQ=255)
              tttACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCa                                                                                                                                  <  1:1491962/150‑1 (MQ=255)
                 aCATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCg                                                                                                                               <  2:765629/150‑1 (MQ=255)
                 aCATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCg                                                                                                                               <  2:1496559/150‑1 (MQ=255)
                    tGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGccc                                                                                                                            >  2:396404/1‑149 (MQ=255)
                       cAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCt                                                                                                                         >  2:624335/1‑150 (MQ=255)
                        aGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTc                                                                                                                        >  1:872116/1‑150 (MQ=255)
                          gcggcggcTGCACGCTACCTGGCCGCTCGAACGTGCTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCgg                                                                                                                      <  1:323175/150‑1 (MQ=255)
                          gcggcggcTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCgg                                                                                                                      >  2:636208/1‑150 (MQ=255)
                             gcggcTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTgg                                                                                                                   <  1:1378930/150‑1 (MQ=255)
                              cggcTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGa                                                                                                                  <  2:348882/150‑1 (MQ=255)
                              cggcTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGa                                                                                                                  >  2:1325698/1‑150 (MQ=255)
                                gcTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGaaa                                                                                                                <  2:452249/150‑1 (MQ=255)
                                 cTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAAc                                                                                                               <  1:1374240/150‑1 (MQ=255)
                                 cTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAAc                                                                                                               >  1:1449405/1‑150 (MQ=255)
                                    cACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAAccgc                                                                                                            <  1:727605/150‑1 (MQ=255)
                                     aCGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAATGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCCCCTCGGTGGCAAccccc                                                                                                           >  1:755027/1‑150 (MQ=255)
                                      cGCTACCTGGCCGCGCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCCCCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCATGGGCATCGCCACCTCGGTGGAAAccgccg                                                                                                          <  2:826729/150‑1 (MQ=255)
                                           ccTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTa                                                                                                     >  1:470929/1‑150 (MQ=255)
                                            cTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAACATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTAc                                                                                                    >  1:115458/1‑150 (MQ=255)
                                                ccGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACtctc                                                                                                >  2:160297/1‑150 (MQ=255)
                                                 cGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCCCCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTcg                                                                                               <  1:310923/150‑1 (MQ=255)
                                                  gCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTcgc                                                                                              <  1:738255/150‑1 (MQ=255)
                                                     cGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTa                                                                                           >  2:1348674/1‑150 (MQ=255)
                                                     cGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTa                                                                                           <  2:872116/150‑1 (MQ=255)
                                                     cGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTa                                                                                           >  1:407737/1‑150 (MQ=255)
                                                      gAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTaa                                                                                          >  2:807585/1‑150 (MQ=255)
                                                      gAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTaa                                                                                          <  1:1011481/150‑1 (MQ=255)
                                                         aGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACACTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCAGTACTCACGAGTAAAGc                                                                                       >  1:256487/1‑150 (MQ=255)
                                                            gtTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACGCTCGCGTAAAGCCAt                                                                                    >  2:399648/1‑150 (MQ=255)
                                                            gtTAATACCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCAt                                                                                    <  1:1392785/150‑1 (MQ=255)
                                                              tAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCAttt                                                                                  >  2:1373716/1‑150 (MQ=255)
                                                               aaTGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCAtttt                                                                                 >  2:1435906/1‑150 (MQ=255)
                                                                 tGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTgc                                                                               >  2:1372715/1‑150 (MQ=255)
                                                                  gCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTgcg                                                                              <  1:1336346/150‑1 (MQ=255)
                                                                     gTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGccc                                                                           <  2:1517716/150‑1 (MQ=255)
                                                                       cAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGa                                                                         >  2:109215/1‑150 (MQ=255)
                                                                        aGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGAt                                                                        <  2:1180532/150‑1 (MQ=255)
                                                                            ggCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGc                                                                    >  1:968642/1‑150 (MQ=255)
                                                                              cTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAg                                                                  <  1:775249/150‑1 (MQ=255)
                                                                              cTACGAAGGCGTGGAGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAg                                                                  <  1:297172/150‑1 (MQ=255)
                                                                                aCGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAgcc                                                                >  1:933889/1‑150 (MQ=255)
                                                                                aCGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGGGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAgcc                                                                <  1:1435906/150‑1 (MQ=255)
                                                                                 cGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAgccg                                                               >  1:546261/1‑150 (MQ=255)
                                                                                  gAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAgccgc                                                              >  2:429127/1‑150 (MQ=255)
                                                                                   aaGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAgccgcc                                                             >  1:1604662/1‑150 (MQ=255)
                                                                                      gCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATc                                                          <  1:1353770/150‑1 (MQ=255)
                                                                                       cGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATcc                                                         >  2:642862/1‑150 (MQ=255)
                                                                                        gtggtgAAATTCGTCTTCGAAAACATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATccc                                                        <  2:115458/150‑1 (MQ=255)
                                                                                        gtggggAAATTCGGCTTCGAAAATATCTCCACCCGCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATccc                                                        <  2:1498535/150‑1 (MQ=255)
                                                                                            tgATATTCGTCTTCGAAAATATCTCGACCCTCGCCGAAAACGCCTGTCGACCGGTACTGGTGGGCCTGGGCATCGCCACCTCTGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGGCCGATTGGCAGCCGCGCTCCCACTc                                                    >  2:759770/1‑150 (MQ=255)
                                                                                            tgAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAAtg                                                    >  2:1259834/1‑149 (MQ=255)
                                                                                            tgAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCAGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCCATc                                                    >  1:1253612/1‑150 (MQ=255)
                                                                                             gAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATcc                                                   <  1:1356020/150‑1 (MQ=255)
                                                                                             gAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATcc                                                   >  2:176788/1‑150 (MQ=255)
                                                                                                aTTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCaaa                                                <  1:1323230/150‑1 (MQ=255)
                                                                                                  tCGTCTTCGAAAATATCTCCACCCTCGCCGTACACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAg                                              <  2:715962/150‑1 (MQ=255)
                                                                                                  tCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAg                                              <  1:1664826/150‑1 (MQ=255)
                                                                                                      cTTCGAAAATATCTCCACCCTCGCCGTAAACGCATGTCCACAGGTACTGGTGGGCGTGCGCATCGCCACCTCGGTGGAAACCGCCGCCGTACACTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGACGCCATCCCAATCCAAAAGAGGc                                          >  2:1033317/1‑150 (MQ=255)
                                                                                                         cGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCACTCCAAAAGCGGCAGa                                       >  1:26628/1‑150 (MQ=255)
                                                                                                         cGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGa                                       <  2:1383958/150‑1 (MQ=255)
                                                                                                          gAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGaa                                      >  2:1242002/1‑150 (MQ=255)
                                                                                                          gAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGaa                                      >  2:140120/1‑150 (MQ=255)
                                                                                                          gAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGaa                                      >  1:1414284/1‑150 (MQ=255)
                                                                                                          gAAAATATCTCAACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGGGGAAAACGCCGCCGGACTCTCGCGTAACGCCATGTTGCGCCCGCTTGGCAGCCGCCATCCCAATCCAAAAGCGGCAgca                                      >  1:1197236/1‑148 (MQ=255)
                                                                                                                 tctcCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGc                               <  2:968642/150‑1 (MQ=255)
                                                                                                                 tctcCACCCTCGCCGTAAACGCCTGTCCACCGGGACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGc                               <  2:1449405/150‑1 (MQ=255)
                                                                                                                      aCCCTCGCCGTAAACGCCTGGCCAGCGGTACTGGTGGGCGTGGGCTTCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTATAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGc                          <  2:256487/150‑1 (MQ=255)
                                                                                                                       cccTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGcc                         <  2:699558/150‑1 (MQ=255)
                                                                                                                         cTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTg                       >  2:458556/1‑150 (MQ=255)
                                                                                                                         cTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTg                       <  1:1372715/150‑1 (MQ=255)
                                                                                                                          tCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTgg                      <  1:624335/150‑1 (MQ=255)
                                                                                                                             ccGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATAGCCACCTCGGTGGCAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGCGCTACGCCTGgaag                   >  1:1657613/1‑150 (MQ=255)
                                                                                                                              cGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGgaaga                  >  2:420795/1‑150 (MQ=255)
                                                                                                                                 aaaCGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAgg               <  1:160297/150‑1 (MQ=255)
                                                                                                                                  aaCGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGa              >  1:542561/1‑150 (MQ=255)
                                                                                                                                   aCGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAACGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAACGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGAc             >  1:726059/1‑150 (MQ=255)
                                                                                                                                   aCGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGAc             <  1:176788/150‑1 (MQ=255)
                                                                                                                                    cGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACt            <  2:997873/150‑1 (MQ=255)
                                                                                                                                      ccTGTCCACCGGTACTTGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACTCa          <  2:768979/150‑1 (MQ=255)
                                                                                                                                          gcccccGGTACTGGTGGGCGTGGGCAGCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACTCAACCg      <  2:1072716/146‑1 (MQ=255)
                                                                                                                                            cACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAAATGGAGCTACGCCTGGAAGAAGGACTACACCGTc    >  1:1112086/1‑150 (MQ=255)
                                                                                                                                              ccGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACTCAACCGTCTg  <  2:755027/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGGCCGCTCGAAAGTGTTAATGCCGTCAGAAGGCTACGAAGGCGTGGTGAAATTCGTCTTCGAAAATATCTCCACCCTCGCCGTAAACGCCTGTCCACCGGTACTGGTGGGCGTGGGCATCGCCACCTCGGTGGAAACCGCCGCCGTACTCTCGCGTAAAGCCATTTTGCGCCCGATTGGCAGCCGCCATCCCAATCCAAAAGCGGCAGAACTGGAGCTACGCCTGGAAGAAGGACTCAACCGTCTG  >  NC_000913/3206894‑3207185

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