Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,749 T→A 6.4% T100S (ACA→TCA) ‡ allE ← (S)‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,7490TA 6.4% 324.8 / 8.1 102T100S (ACA→TCA) ‡allE(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (51/44);  new base A (3/4);  total (54/48)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGA  >  NC_000913/543600‑543896
                                                                                                                                                     |                                                                                                                                                   
gCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGt                                                                                                                                                     >  1:895355/1‑150 (MQ=255)
gCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGt                                                                                                                                                     >  1:275789/1‑150 (MQ=255)
gCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGt                                                                                                                                                     >  1:194427/1‑150 (MQ=255)
 ccAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGtt                                                                                                                                                    >  2:518737/1‑150 (MQ=255)
  cAGAAACCAGCCCCGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGCGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGgtt                                                                                                                                                   >  2:440988/1‑150 (MQ=255)
      aaCCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGcc                                                                                                                                               <  1:1534823/150‑1 (MQ=255)
       aCCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCt                                                                                                                                              >  1:992263/1‑150 (MQ=255)
        ccAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCtt                                                                                                                                             >  2:463115/1‑150 (MQ=255)
        ccAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCtt                                                                                                                                             >  1:612578/1‑150 (MQ=255)
        ccAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCtt                                                                                                                                             >  1:450956/1‑150 (MQ=255)
         gtggctcGGTGCATAGCCTTCTAGCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCGTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCtta                                                                                                                                            <  2:1317025/144‑2 (MQ=255)
         cAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTc                                                                                                                                            <  1:1064505/150‑1 (MQ=255)
           gCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTTGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATATAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCttaag                                                                                                                                          <  2:911987/150‑4 (MQ=255)
           gCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCgg                                                                                                                                          >  2:1350130/1‑150 (MQ=255)
            ccACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGc                                                                                                                                         <  2:174801/150‑1 (MQ=255)
            ccACGGTGCATAGCCTTCTACCCGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCttaagc                                                                                                                                         <  1:884389/150‑5 (MQ=255)
             cACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTTGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGACCATGGCTGGCAATAAAGATAGCCACCTTCGCTTAATCCAAATGTTTTTCCTTCGGCt                                                                                                                                        >  2:549305/1‑150 (MQ=255)
             cACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCttaagcg                                                                                                                                        <  1:850518/150‑6 (MQ=255)
               cGGTGCATAGCCTTCTACCGGAACATAGCGGCTCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGATCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCttt                                                                                                                                      <  1:549305/150‑1 (MQ=255)
                     aTAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGt                                                                                                                                <  1:53507/150‑1 (MQ=255)
                       aGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGa                                                                                                                              >  2:519122/1‑150 (MQ=255)
                              tACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTcc                                                                                                                       <  1:1088078/150‑1 (MQ=255)
                               aCCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCa                                                                                                                      >  2:1398622/1‑150 (MQ=255)
                               aCCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCa                                                                                                                      <  1:250012/150‑1 (MQ=255)
                                ccGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCag                                                                                                                     >  1:72793/1‑150 (MQ=255)
                                ccGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCag                                                                                                                     >  2:295161/1‑150 (MQ=255)
                                  ggAACATAGCGGCTCTTATATAAAAAGATTTCGCTGTCTTCGGCCAGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGATTTGCCTTCGGCTTTGGCAGTGATATTTCCagag                                                                                                                   >  2:874204/1‑150 (MQ=255)
                                  ggAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCagag                                                                                                                   >  2:1510436/1‑150 (MQ=255)
                                  ggAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCagag                                                                                                                   >  2:1505331/1‑150 (MQ=255)
                                     aCATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATc                                                                                                                <  1:300216/150‑1 (MQ=255)
                                      cATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATca                                                                                                               <  1:1398622/150‑1 (MQ=255)
                                        tAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATcaca                                                                                                             <  1:54204/150‑1 (MQ=255)
                                         aGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACAt                                                                                                            <  2:1270790/150‑1 (MQ=255)
                                             gcgcTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAg                                                                                                        <  2:1438583/150‑1 (MQ=255)
                                              cgcTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAgg                                                                                                       <  2:1357350/150‑1 (MQ=255)
                                                  tatataAAAAGATTTTGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACg                                                                                                   <  1:964392/150‑1 (MQ=255)
                                                  tatataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACg                                                                                                   >  2:596662/1‑150 (MQ=255)
                                                    tataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGtt                                                                                                 <  1:1526209/150‑1 (MQ=255)
                                                     ataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGTGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACATGAACGttt                                                                                                <  1:396843/150‑1 (MQ=255)
                                                      taAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTc                                                                                               >  1:1164335/1‑150 (MQ=255)
                                                      taAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTc                                                                                               >  2:1205085/1‑150 (MQ=255)
                                                      taAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGACTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACCTACAGGAACGTTTc                                                                                               >  1:1154772/1‑150 (MQ=255)
                                                       aaaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCa                                                                                              >  1:588264/1‑150 (MQ=255)
                                                        aaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCaa                                                                                             >  2:142965/1‑150 (MQ=255)
                                                        aaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCaa                                                                                             <  1:479117/150‑1 (MQ=255)
                                                        aaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCaa                                                                                             <  2:275789/150‑1 (MQ=255)
                                                                 ggCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGc                                                                                    >  2:822086/1‑150 (MQ=255)
                                                                   cTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGttt                                                                                                >  1:1170272/1‑136 (MQ=255)
                                                                   cTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGttt                                                                                                <  2:1170272/136‑1 (MQ=255)
                                                                   cTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCa                                                                                  >  1:1239327/1‑150 (MQ=255)
                                                                     gTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTGGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCAcc                                                                                <  2:1190252/150‑1 (MQ=255)
                                                                           ggCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGcc                                                                          >  1:588667/1‑150 (MQ=255)
                                                                           ggCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGcc                                                                          >  1:363460/1‑150 (MQ=255)
                                                                              cTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTg                                                                       >  1:223576/1‑150 (MQ=255)
                                                                                     ttATCAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCATATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTgttggttg                                                                <  1:973978/150‑1 (MQ=255)
                                                                                       aaCAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGcc                                                              >  2:1186484/1‑150 (MQ=255)
                                                                                       aaCAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGcc                                                              <  1:1350130/150‑1 (MQ=255)
                                                                                       aaCAAACGTCATTAAGGAGCCAGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTGCGGCTTTGGCAGTGATATTTACAGAGATCACATACAGGAACGTTTCAATCCCTTCGCACCCGACGCCCTgttggttgtt                                                              >  2:1173105/1‑148 (MQ=255)
                                                                                         cAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCAc                                                            >  2:671934/1‑150 (MQ=255)
                                                                                             cGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGtt                                                        >  2:32452/1‑150 (MQ=255)
                                                                                             cGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGtt                                                        <  2:1025614/150‑1 (MQ=255)
                                                                                             cGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGtt                                                        >  2:1441518/1‑150 (MQ=255)
                                                                                               tCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGtttt                                                      >  1:1092856/1‑150 (MQ=255)
                                                                                                cATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACCTACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTg                                                     <  2:1154772/150‑1 (MQ=255)
                                                                                                cATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTg                                                     >  2:812014/1‑150 (MQ=255)
                                                                                                cATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTg                                                     <  2:342242/150‑1 (MQ=255)
                                                                                                   tAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATg                                                  >  1:189338/1‑150 (MQ=255)
                                                                                                    aaGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGc                                                 <  2:612578/150‑1 (MQ=255)
                                                                                                      ggAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAg                                               >  2:182343/1‑150 (MQ=255)
                                                                                                       gAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAgt                                              >  1:208367/1‑150 (MQ=255)
                                                                                                        aGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAgtg                                             >  1:723749/1‑150 (MQ=255)
                                                                                                          ccTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGtgtg                                           <  2:624250/150‑1 (MQ=255)
                                                                                                           cTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCAACGTTTTGATGCAGtgtgc                                          >  2:1453947/1‑149 (MQ=255)
                                                                                                            tGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGAc                                         >  1:716696/1‑150 (MQ=255)
                                                                                                             ggCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGAcc                                        <  2:1092856/150‑1 (MQ=255)
                                                                                                              gCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCa                                       <  1:812014/150‑1 (MQ=255)
                                                                                                                 ggCAATAAAGATAGCCACCTTCTCTAAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGAt                                    <  2:803097/150‑1 (MQ=255)
                                                                                                                 ggCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGAt                                    <  2:1105978/150‑1 (MQ=255)
                                                                                                                 ggCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGAt                                    <  2:895355/150‑1 (MQ=255)
                                                                                                                 ggCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGAAGTGTGACCAGAt                                    >  1:673145/1‑150 (MQ=255)
                                                                                                                     aTAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATc                                >  2:663654/1‑150 (MQ=255)
                                                                                                                     aTAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATc                                >  1:1317106/1‑150 (MQ=255)
                                                                                                                       aaaGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATcaa                              >  1:1061505/1‑150 (MQ=255)
                                                                                                                             aGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAag                        <  1:32452/150‑1 (MQ=255)
                                                                                                                              gCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAaga                       <  2:378642/150‑1 (MQ=255)
                                                                                                                                cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGc                                                 <  2:379119/122‑1 (MQ=255)
                                                                                                                                cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGc                                                 >  1:338574/1‑122 (MQ=255)
                                                                                                                                cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGc                                                 <  2:338574/122‑1 (MQ=255)
                                                                                                                                cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAgg                     >  1:1009328/1‑150 (MQ=255)
                                                                                                                                       gCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAg              <  1:1205085/150‑1 (MQ=255)
                                                                                                                                        ccgaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGc             >  1:1267841/4‑150 (MQ=255)
                                                                                                                                         cgaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCt            >  1:23795/3‑150 (MQ=255)
                                                                                                                                          tAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCtt           >  2:6937/1‑150 (MQ=255)
                                                                                                                                            aGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTg         <  2:588264/150‑1 (MQ=255)
                                                                                                                                             ggCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTgg        >  2:409908/1‑150 (MQ=255)
                                                                                                                                             ggCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTgg        >  1:1437477/1‑150 (MQ=255)
                                                                                                                                             ggCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTgg        >  2:983623/1‑150 (MQ=255)
                                                                                                                                             ggCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGGTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTgg        <  1:519122/150‑1 (MQ=255)
                                                                                                                                               cAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCg      <  1:175452/150‑1 (MQ=255)
                                                                                                                                                   tGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGga  >  2:225175/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGA  >  NC_000913/543600‑543896

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