Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,062 T→A 7.9% intergenic (+29/+199) glxK → / ← allE glycerate kinase II/S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,0620TA 7.9% 329.0 / 13.8 122intergenic (+29/+199)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (57/55);  new base A (7/3);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 123 total aligned reads displayed.

AGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTA  >  NC_000913/542914‑543207
                                                                                                                                                    |                                                                                                                                                 
aGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACATTTAAAGCGGATGCGATTtata                                                                                                                                                  >  2:688440/1‑150 (MQ=255)
  cATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATc                                                                                                                                                <  2:1410955/150‑1 (MQ=255)
  cATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATc                                                                                                                                                >  1:1438068/1‑150 (MQ=255)
      ttGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCAt                                                                                                                                            <  1:398632/150‑1 (MQ=255)
      ttGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAt                                                                                                                                            <  2:393702/150‑1 (MQ=255)
       tGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATa                                                                                                                                           <  1:1172765/150‑1 (MQ=255)
        gCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATaa                                                                                                                                          >  1:792338/1‑150 (MQ=255)
        gCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAtcc                                                                                                                                          <  1:1730742/150‑3 (MQ=255)
        gCCTCGTCTGGCACCTTTAGCCGAAGGGCGCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAtcc                                                                                                                                          <  2:200514/150‑3 (MQ=255)
         ccTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGATATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAag                                                                                                                                         >  1:1052267/1‑150 (MQ=255)
           tCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAagag                                                                                                                                       <  2:773211/150‑1 (MQ=255)
             gTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTg                                                                                                                                     >  1:812112/1‑150 (MQ=255)
             gTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGAGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTg                                                                                                                                     >  2:1704809/1‑150 (MQ=255)
               cTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCa                                                                                                                                   >  1:1512688/1‑150 (MQ=255)
                 ggCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTACCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCCGt                                                                                                                                 >  1:1158532/1‑150 (MQ=255)
                 ggCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACACTTACCTTTAAAGCGGATGCGAATTATATCTCATAAGAGTGCAGt                                                                                                                                 >  1:1164479/1‑150 (MQ=255)
                   cACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTACCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTAc                                                                                                                               >  2:129854/1‑150 (MQ=255)
                    aCCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACt                                                                                                                              <  2:1443243/150‑1 (MQ=255)
                     ccTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTc                                                                                                                             >  1:163346/1‑150 (MQ=255)
                      cTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCa                                                                                                                            <  2:512206/150‑1 (MQ=255)
                      cTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCa                                                                                                                            >  1:89464/1‑150 (MQ=255)
                          aGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGcc                                                                                                                        <  2:539041/150‑1 (MQ=255)
                            ccGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCgg                                                                                                                      >  1:1712089/1‑150 (MQ=255)
                              gAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGAt                                                                                                                    >  1:782121/1‑150 (MQ=255)
                                ccaatacGCCAGCGTTGAAACCAATCTCATCAACAGCGAGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                  <  1:1483791/144‑1 (MQ=255)
                                 gTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATATAACGCATAACAGTGCAGTACTCATTCCGGATGCg                                                                                                                 >  2:342866/1‑150 (MQ=255)
                                   gCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc                                                                                                               >  1:417447/1‑150 (MQ=255)
                                      cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATg                                                                                                            <  2:302698/150‑1 (MQ=255)
                                      cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATg                                                                                                            <  2:163346/150‑1 (MQ=255)
                                      cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGCAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATTCCGGATGCGGCAtc                                                                                                            >  2:654618/1‑149 (MQ=255)
                                         cAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGt                                                                                                         <  1:1350799/150‑1 (MQ=255)
                                           gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAc                                                                                                       <  2:1259978/150‑1 (MQ=255)
                                           gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAc                                                                                                       <  2:1438068/150‑1 (MQ=255)
                                             ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCa                                                                                                     <  1:984309/150‑1 (MQ=255)
                                              gTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCat                                                                                                    >  2:1079094/1‑150 (MQ=255)
                                              gTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCat                                                                                                    >  1:286349/1‑150 (MQ=255)
                                                gAAACCAATCTCATCAACAGCGCGCGAAATATTGACTGCGCCATTAAAATAGGGCAGGGAATTAAACACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCatat                                                                                                  <  1:49123/150‑1 (MQ=255)
                                                 aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                 <  1:128546/150‑1 (MQ=255)
                                                  aaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcc                                                                                                >  2:865369/1‑150 (MQ=255)
                                                  aaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcc                                                                                                >  2:187758/1‑150 (MQ=255)
                                                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                             >  2:1120054/1‑150 (MQ=255)
                                                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                             >  1:766344/1‑150 (MQ=255)
                                                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCAGTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                             <  1:100541/150‑1 (MQ=255)
                                                      aaTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                            <  2:1706714/150‑1 (MQ=255)
                                                       aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcct                                                                                           <  2:1493282/150‑1 (MQ=255)
                                                       aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcct                                                                                           <  2:792338/150‑1 (MQ=255)
                                                       aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTATACAGCGGATGCGATTGATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGAACCATAtccttcct                                                                                           >  1:240805/1‑150 (MQ=255)
                                                               aaCAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcg                                                                                   <  2:611484/150‑1 (MQ=255)
                                                                  aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCgcgc                                                                                <  1:564513/150‑1 (MQ=255)
                                                                     gcgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCaaa                                                                             >  1:1206136/1‑150 (MQ=255)
                                                                      cgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCTCGCAAAt                                                                            >  2:1697389/1‑150 (MQ=255)
                                                                      cgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAAt                                                                            >  1:43993/1‑150 (MQ=255)
                                                                          aaaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                        >  1:320966/1‑150 (MQ=255)
                                                                           aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGTATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtat                                                                       <  2:276100/150‑1 (MQ=255)
                                                                            atatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                      <  1:795753/150‑1 (MQ=255)
                                                                                  ccTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCTTATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc                                                                <  2:1401558/150‑1 (MQ=255)
                                                                                    tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCag                                                              <  2:1052267/150‑1 (MQ=255)
                                                                                    tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCag                                                              <  2:1634166/150‑1 (MQ=255)
                                                                                    tGCGACTTTAAAATAGTTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCag                                                              <  1:23790/150‑1 (MQ=255)
                                                                                     gcgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga                                                             >  2:238985/1‑150 (MQ=255)
                                                                                       gcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAgaga                                                           <  2:1005723/150‑1 (MQ=255)
                                                                                        ccATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCAGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                          >  2:457286/1‑150 (MQ=255)
                                                                                         cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                         >  1:1102752/1‑150 (MQ=255)
                                                                                         cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATACTTACTGAATATCGCGCAAATTCTATTTTTTGCagagtta                                                         >  2:741155/1‑147 (MQ=255)
                                                                                          aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                        >  2:622395/1‑150 (MQ=255)
                                                                                          aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                        <  2:701642/150‑1 (MQ=255)
                                                                                           ttAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                       <  1:1370403/150‑1 (MQ=255)
                                                                                            tAAAATAGGTCATTGAATTAAAAACAAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                      <  1:688440/150‑1 (MQ=255)
                                                                                              aaaTAGGTCAGGGAATTAAAAACTACCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                                    >  2:258234/1‑150 (MQ=255)
                                                                                              aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                                    <  2:89464/150‑1 (MQ=255)
                                                                                                   ggTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                               >  2:365867/1‑150 (MQ=255)
                                                                                                   ggTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATAGCAGAGATCATGTACGCCt                                               >  2:1657782/1‑150 (MQ=255)
                                                                                                      cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                            >  1:1115190/1‑150 (MQ=255)
                                                                                                      cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                            >  2:1064957/1‑150 (MQ=255)
                                                                                                      cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCAATATATTGCAGAGATCATGTAGGCCTGAt                                            >  2:1509126/1‑150 (MQ=255)
                                                                                                        gggAATTTAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATACATCCTGAAATTCGCGCAAAATCTATATATTTCAGAGATCATGTAGGCCTGATaa                                          >  1:1051768/1‑150 (MQ=255)
                                                                                                        gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                          >  2:1707612/1‑150 (MQ=255)
                                                                                                        gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                          <  1:1657851/150‑1 (MQ=255)
                                                                                                              tAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaag                                    <  2:1115190/150‑1 (MQ=255)
                                                                                                                aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGAGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCATAGATCATGTAGGCCTGATaagcgaagcg                                  >  2:1355264/1‑150 (MQ=255)
                                                                                                                    cTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATc                              <  1:1704809/150‑1 (MQ=255)
                                                                                                                         ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCTAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCa                         >  2:601280/1‑150 (MQ=255)
                                                                                                                         ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCa                         <  2:1411947/150‑1 (MQ=255)
                                                                                                                            tACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATg                      <  2:844801/150‑1 (MQ=255)
                                                                                                                               cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                   >  2:1185999/1‑150 (MQ=255)
                                                                                                                               cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                   >  2:1179201/1‑150 (MQ=255)
                                                                                                                                   aaaGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTTAAAATCGCGCAAATTCTATATATTGAAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaa               >  2:677756/1‑150 (MQ=255)
                                                                                                                                   aaaGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaa               >  2:1244009/1‑150 (MQ=255)
                                                                                                                                    actcttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaa              >  1:998895/7‑150 (MQ=255)
                                                                                                                                      gCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaa            >  1:407346/1‑150 (MQ=255)
                                                                                                                                       cttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa           >  1:1489034/4‑150 (MQ=255)
                                                                                                                                       cttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa           >  2:1566170/4‑150 (MQ=255)
                                                                                                                                       cttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATAATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa           >  1:1603464/4‑150 (MQ=255)
                                                                                                                                        ttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg          >  1:595136/3‑150 (MQ=255)
                                                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCTTGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg          <  1:840114/150‑1 (MQ=255)
                                                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAg          <  2:1508713/150‑1 (MQ=255)
                                                                                                                                          aTGCGATATAAATCGCATAAGATTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAACCTTCCTGAAAATCGCGCAAATTCTATATATTGCACAGATCATGTAGGCATGATAAGCGAAACGCATCACGCAATGTTACAAAAAAAGcc        >  2:835238/1‑150 (MQ=255)
                                                                                                                                                tttATATCGCATAAGAGTGCAGTACTCATGCCGTATGATGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta  <  1:1697389/150‑1 (MQ=255)
                                                                                                                                                tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAGCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta  <  2:1164479/150‑1 (MQ=255)
                                                                                                                                                tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta  <  2:1158532/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGCATTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTA  >  NC_000913/542914‑543207

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: