Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,901 C→T 5.7% intergenic (‑98/+27) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9010CT 5.7% 367.3 / 6.7 122intergenic (‑98/+27)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (70/44);  new base T (4/3);  total (75/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 127 total aligned reads displayed.

CGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGG  >  NC_000913/1260752‑1261050
                                                                                                                                                      |                                                                                                                                                     
 gTCGATCAGAGCGCGTAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACAAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgccg                                                                                                                                                      <  1:168875/150‑1 (MQ=255)
   cGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCa                                                                                                                                                    >  2:746704/1‑150 (MQ=255)
    gATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAg                                                                                                                                                   >  1:1047736/1‑150 (MQ=255)
       cAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCgg                                                                                                                                                >  1:107489/1‑150 (MQ=255)
        agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                                                                               >  1:393901/1‑150 (MQ=255)
        agagCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                                                                               <  2:1887330/150‑1 (MQ=255)
            cgcgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtt                                                                                                                                           >  1:1059320/1‑150 (MQ=255)
            cgcgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCtt                                                                                                                                           <  1:1857085/150‑1 (MQ=255)
              cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttt                                                                                                                                         <  1:2017486/150‑2 (MQ=255)
              cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                                                                         >  2:1247986/1‑150 (MQ=255)
              cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGGGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                                                                         <  1:1400221/150‑1 (MQ=255)
              cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAAAAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCAAACAAATGGATAATACAGACAAAAAAACCCGCAGCAGCGGCTCTTTg                                                                                                                                         >  2:201326/1‑150 (MQ=255)
              cagAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                                                                         <  1:848068/148‑1 (MQ=255)
                gAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                                                                       <  1:1522391/150‑1 (MQ=255)
                 aaaGGCATCAGATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAATTGGCAGGTGAAGGAGGCATAGGTCATACTAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttttgt                                                                                                                                      <  1:1604795/150‑5 (MQ=255)
                 aaaGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGc                                                                                                                                      <  1:1313324/150‑1 (MQ=255)
                    ggCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCgg                                                                                                                                   >  2:229459/1‑150 (MQ=255)
                     gCATCACATTTGAGGAAAATATTTTGCTCACAAAAATGTTTCATCCGTGAGCGCATCATCTGTCAATTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCggg                                                                                                                                  <  2:372103/150‑1 (MQ=255)
                      cATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                                                                                 >  2:945413/1‑150 (MQ=255)
                      cATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGt                                                                                                                                 >  2:899299/1‑150 (MQ=255)
                       aTCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                                                                                <  2:1723527/150‑1 (MQ=255)
                          acaTGTTAGGAAAATATTTTGATCACAATAATCTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGa                                                                                                                             <  1:335514/150‑1 (MQ=255)
                            aTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                                                                           <  2:1866141/150‑1 (MQ=255)
                             tgtgAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCCACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTAAAGACAAAAAAACCCGCCGCAGCGGGTCTTTTCTCCGGGTGCGTTTa                                                                                                                          >  2:1941555/1‑150 (MQ=255)
                              gtgAGGAAAATATTTTGTTCACAATAATTTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                                                                         <  2:242666/150‑1 (MQ=255)
                               tgAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                                                                        >  2:1082019/1‑150 (MQ=255)
                                  ggAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                                                                     >  2:738996/1‑150 (MQ=255)
                                    aaaaTATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCg                                                                                                                   >  1:608248/1‑150 (MQ=255)
                                    aaaaTATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGACGCAGCGGGTCTTTCAGCAGGGTTCGAATAGTGTTCg                                                                                                                   >  2:495517/1‑150 (MQ=255)
                                       atatTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAAc                                                                                                                >  2:1484377/1‑150 (MQ=255)
                                           tttGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATgg                                                                                                            <  2:1106656/150‑1 (MQ=255)
                                           tttGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATgg                                                                                                            <  2:1894154/150‑1 (MQ=255)
                                            ttGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGc                                                                                                           <  1:158395/150‑1 (MQ=255)
                                            ttGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGc                                                                                                           >  2:544582/1‑150 (MQ=255)
                                                tCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCagag                                                                                                       >  2:1784439/1‑150 (MQ=255)
                                                tCACAATAATGTTTCATCCGTGAGCGCATCATATGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAAAGGATATTACACAAAAAAAAACACGACGCAGCTGGTCTTTGATCCCGGTTCGATTATTGTTCGATCATGCCagag                                                                                                       >  1:202117/1‑150 (MQ=255)
                                                 cacaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAgaga                                                                                                      >  2:964580/1‑150 (MQ=255)
                                                   caATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATc                                                                                                    >  2:1590416/1‑150 (MQ=255)
                                                      taatGTTTCATCCGTGAGCGCATCATCTGTCCACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAACGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAt                                                                                                 <  2:1742233/150‑1 (MQ=255)
                                                       aatGTTTGAGCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                                <  1:1106307/150‑1 (MQ=255)
                                                       aatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                                <  2:1677551/150‑1 (MQ=255)
                                                       aatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCCGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAtt                                                                                                >  2:614775/1‑150 (MQ=255)
                                                             tCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGt                                                                                          <  2:78935/150‑1 (MQ=255)
                                                              cATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                                                                         >  2:1494074/1‑150 (MQ=255)
                                                                tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                                                                       <  2:1968361/150‑1 (MQ=255)
                                                                  cGTGATCGCATCATCTGTCAACTAAACGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTg                                                                                     <  2:619660/150‑1 (MQ=255)
                                                                  cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCCCCCCAGCGGGGCAATGAGCCGGGTTCGATTAGTTTTCCAAAATGGCAGAGAACGATTCATCGTTGCTg                                                                                     >  2:607872/1‑150 (MQ=255)
                                                                     gAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATa                                                                                  >  2:214236/1‑150 (MQ=255)
                                                                      aGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATac                                                                                 >  2:17885/1‑150 (MQ=255)
                                                                         gcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATAATACAGACAAAAAAACCCGCCGAAGCGGGTCTTTGAGCCGGGTTCGATAATTGTTCGATCATGCCAGAGATCGAATCTTCGTTGCTGATacgac                                                                              >  1:882100/1‑150 (MQ=255)
                                                                          catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgacg                                                                             <  2:395836/150‑1 (MQ=255)
                                                                              atcTGTCAACTAAATGGCCGGCGAATGAGGCATAGGTCCTACAAATGGATATTACAGCCAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATc                                                                         <  2:380638/150‑1 (MQ=255)
                                                                                   tCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                                                                    <  1:1487508/150‑1 (MQ=255)
                                                                                    cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGAACGATTCTACGTTGCTGATACGACGAATCGCTACg                                                                   >  2:294937/1‑150 (MQ=255)
                                                                                         aaaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcaac                                                              >  1:1500141/1‑147 (MQ=255)
                                                                                         aaaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg                                                              >  1:1894794/1‑150 (MQ=255)
                                                                                          aaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc                                                             <  1:1484377/150‑1 (MQ=255)
                                                                                             ggCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGACGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATa                                                          >  2:1972791/1‑150 (MQ=255)
                                                                                                aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt                                                       >  1:1853077/1‑150 (MQ=255)
                                                                                                aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt                                                       >  1:625573/1‑150 (MQ=255)
                                                                                                  gTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGa                                                     >  1:1412181/1‑150 (MQ=255)
                                                                                                   tGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAc                                                    >  1:1720769/1‑150 (MQ=255)
                                                                                                     aaGGAGGCATAGGTCATACAAATGGATATTACAGGCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGATCAAGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAg                                                  >  1:868614/1‑150 (MQ=255)
                                                                                                     aaGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAg                                                  <  1:371976/150‑1 (MQ=255)
                                                                                                      aggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAgg                                                 >  2:452988/1‑150 (MQ=255)
                                                                                                       ggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAggg                                                <  1:1784439/150‑1 (MQ=255)
                                                                                                         aggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTc                                              <  2:407270/150‑1 (MQ=255)
                                                                                                         aggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGATTCGGCCAGCATACCTGACAGGGTc                                              >  1:1224834/1‑150 (MQ=255)
                                                                                                          ggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCa                                             >  2:50782/1‑150 (MQ=255)
                                                                                                              tAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGt                                         <  2:1047736/150‑1 (MQ=255)
                                                                                                                  tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGc                                     >  2:175392/1‑150 (MQ=255)
                                                                                                                   cATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGcc                                    >  1:1675093/1‑149 (MQ=255)
                                                                                                                    aTACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCAc                                   <  1:945413/150‑1 (MQ=255)
                                                                                                                     tACATATGGATATTACTGACAAAAAAACCCGCCGCCGCGGGTCTTTGAGCCGGGCTCGATTCGTGTGCGAACATGGCAGAGATCGATTCTGCGATGCTGATACGACGACTCGCTTCGGCCAGCATACCTGACAGCGTCAGAGTACGCAcc                                  <  1:286921/150‑2 (MQ=255)
                                                                                                                     tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGCGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACg                                  >  1:1480017/1‑150 (MQ=255)
                                                                                                                     tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACg                                  <  2:558141/150‑1 (MQ=255)
                                                                                                                      aCAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGt                                 >  1:705150/1‑150 (MQ=255)
                                                                                                                       cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGtt                                >  2:1401147/1‑150 (MQ=255)
                                                                                                                         aaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCg                              >  1:1929681/1‑150 (MQ=255)
                                                                                                                            ggATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCa                           <  2:494726/150‑1 (MQ=255)
                                                                                                                               ttttACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGCGCCGGGTGCGAGTAGTGTTCGAACATGGCAGCGATCGAGTGTTCGTTGCTGATACGCCGATTAGCTACGGCCAGCATACCTGACAGGGTCAGACTACGCACGTAAGGCAggg                        >  2:320343/3‑148 (MQ=255)
                                                                                                                               tatTACAGACAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGt                         <  1:1674872/150‑1 (MQ=255)
                                                                                                                               tatTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtg                        >  1:2009731/1‑150 (MQ=255)
                                                                                                                                   aCAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgattt                    >  2:266721/1‑150 (MQ=255)
                                                                                                                                   aCAGACAAAAAAACCCGCCGCAGCGGGGATTTGAGCCGGGTTCGATTAGTGTTCGCACATGGCAGAGATCGATTCTTCGTTGCTGAGACGACGAATCGCATCGGCCATCATACATGACAGGGTAAGAGTACGCAGGTTCGGCCGtgattt                    >  2:1493345/1‑150 (MQ=255)
                                                                                                                                      gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                 >  1:585567/1‑150 (MQ=255)
                                                                                                                                      gACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgat                 >  2:637351/1‑150 (MQ=255)
                                                                                                                                        gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGtgatttgatt                >  1:372170/5‑150 (MQ=255)
                                                                                                                                        cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt               >  1:1401265/1‑150 (MQ=255)
                                                                                                                                        cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt               >  2:1988271/1‑150 (MQ=255)
                                                                                                                                        cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt               >  2:1490337/1‑150 (MQ=255)
                                                                                                                                        cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt               >  1:1114918/1‑150 (MQ=255)
                                                                                                                                         ctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTgatttgattt               >  1:1729331/4‑150 (MQ=255)
                                                                                                                                          aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATAGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGATCGGCAGTGATTTGATTTc              >  2:1378662/1‑150 (MQ=255)
                                                                                                                                           caaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCa             >  2:1938987/2‑150 (MQ=255)
                                                                                                                                             aaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACAGGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATc           <  1:1655200/150‑1 (MQ=255)
                                                                                                                                                 cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGAAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTc       >  2:83067/1‑150 (MQ=255)
                                                                                                                                                  ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCa      >  2:240868/1‑150 (MQ=255)
                                                                                                                                                   cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAg     >  1:1151248/1‑150 (MQ=255)
                                                                                                                                                      cgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCgg  <  1:1082019/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGACAGGGTCAGAGTACGCACGTTCGGCAGTGATTTGATTTCATCGCTCAGCGG  >  NC_000913/1260752‑1261050

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: