Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,091 G→C 10.4% intergenic (+58/+170) 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,0910GC10.4% 325.5 / 22.9 96intergenic (+58/+170)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (45/41);  new base C (3/7);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAA  >  NC_000913/542944‑543239
                                                                                                                                                   |                                                                                                                                                    
gAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGa                                                                                                                                                     <  1:664300/149‑1 (MQ=255)
gAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGa                                                                                                                                                     >  2:68404/1‑149 (MQ=255)
      cTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGc                                                                                                                                               >  1:53652/1‑149 (MQ=255)
       tCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCa                                                                                                                                              >  3:832431/1‑149 (MQ=255)
       tCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATCGCCTGCGCCATTAAAATAGGTCAGGCAATTAAAAACTAACCCTTACCTTGAAAGCGGATGCGATTTATATCGCATAAGAGTGCTGTACTCATGCCGCATCCGTCa                                                                                                                                              <  4:574700/149‑1 (MQ=255)
        cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                             >  1:169456/1‑149 (MQ=255)
         gCCAGCGGTGAAACCAACCTCTTCAACAGCGCGCGAAATATTGCCCGCGCCATTAAAATAAGTCAGGTAATTAAAAACTAACCCATACCTTTAAAGCGGATGCGATTTATAGCGCATAAGAGTGCAGAACTCATGCCGCATCCGGCATg                                                                                                                                            <  2:263638/149‑1 (MQ=255)
              cGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAc                                                                                                                                       <  2:475613/149‑1 (MQ=255)
               ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAct                                                                                                                                      <  2:53652/149‑2 (MQ=255)
               ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTAct                                                                                                                                      <  4:884236/149‑2 (MQ=255)
                    aaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                                                 >  2:460231/1‑149 (MQ=255)
                      ccAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcct                                                                                                                               <  1:717439/149‑1 (MQ=255)
                       cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                           >  1:125602/1‑136 (MQ=255)
                       cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                           <  2:125602/136‑1 (MQ=255)
                          tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcct                                                                                                                           <  1:722755/149‑1 (MQ=255)
                                   cAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgc                                                                                                                  <  3:783693/149‑1 (MQ=255)
                                    aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                            <  4:295884/138‑1 (MQ=255)
                                    aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                            >  3:295884/1‑138 (MQ=255)
                                    aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                            >  1:516538/1‑138 (MQ=255)
                                    aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                            <  2:516538/138‑1 (MQ=255)
                                    aGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg                                                                                                                 >  3:159200/1‑149 (MQ=255)
                                     gcgcgcGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                          <  4:501449/123‑1 (MQ=255)
                                     gcgcgcGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAg                                                                                                                                          >  3:501449/1‑123 (MQ=255)
                                       gcgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCaa                                                                                                              >  2:195661/1‑149 (MQ=255)
                                        cgcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCaaa                                                                                                             <  1:68404/149‑1 (MQ=255)
                                            aaaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCt                                                                                                         <  2:92996/149‑1 (MQ=255)
                                             aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                        <  4:159200/149‑1 (MQ=255)
                                             aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                        <  2:169456/149‑1 (MQ=255)
                                              atatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtat                                                                                                       <  3:232711/149‑1 (MQ=255)
                                               tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                                                      >  1:764521/1‑149 (MQ=255)
                                                   gCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATAtt                                                                                                  <  1:844475/149‑1 (MQ=255)
                                                      tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                               >  3:160981/1‑149 (MQ=255)
                                                      tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                               <  2:71114/149‑1 (MQ=255)
                                                      tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                               <  1:748987/149‑1 (MQ=255)
                                                         gcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCagag                                                                                            >  1:539875/1‑149 (MQ=255)
                                                            aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                         >  1:674909/1‑149 (MQ=255)
                                                            aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                         >  1:457072/1‑149 (MQ=255)
                                                                 aaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                                                                    <  4:511539/149‑1 (MQ=255)
                                                                      gTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                                                               >  4:513216/1‑149 (MQ=255)
                                                                       tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                                                  >  2:894141/1‑81 (MQ=255)
                                                                       tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                                                  <  1:894141/81‑1 (MQ=255)
                                                                         aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                            >  4:575253/1‑149 (MQ=255)
                                                                          gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATa                                                                           >  2:778822/1‑149 (MQ=255)
                                                                          gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                           >  3:519135/1‑85 (MQ=255)
                                                                          gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa                                                                                                                                           <  4:519135/85‑1 (MQ=255)
                                                                           ggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                                                          >  4:681500/1‑149 (MQ=255)
                                                                           ggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                                                          <  3:575253/149‑1 (MQ=255)
                                                                                 aaaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATATGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTCAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaag                                                                    <  3:646357/149‑1 (MQ=255)
                                                                                  aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaagc                                                                   >  3:202506/1‑149 (MQ=255)
                                                                                       tAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATc                                                              <  4:160981/149‑1 (MQ=255)
                                                                                          cccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAgg                                                           <  1:778822/149‑1 (MQ=255)
                                                                                            cTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                         >  2:362393/1‑149 (MQ=255)
                                                                                             ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                        >  1:844819/1‑149 (MQ=255)
                                                                                             ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa                                                        <  3:490192/149‑1 (MQ=255)
                                                                                               aCCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATg                                                      >  1:875782/1‑149 (MQ=255)
                                                                                                 cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaag                                                                         <  1:932839/128‑1 (MQ=255)
                                                                                                 cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaag                                                                         >  2:932839/1‑128 (MQ=255)
                                                                                                 cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGtt                                                    >  4:89566/1‑149 (MQ=255)
                                                                                                  tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                                                     <  2:804243/115‑1 (MQ=255)
                                                                                                  tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                                                     >  1:804243/1‑115 (MQ=255)
                                                                                                  tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                                                   <  2:311100/149‑1 (MQ=255)
                                                                                                  tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                                                   >  2:410306/1‑149 (MQ=255)
                                                                                                  tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAGAAGCGAAGCGCATCAGGCAATGTTa                                                   >  3:604418/1‑149 (MQ=255)
                                                                                                   ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTGCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTAc                                                  <  4:604418/149‑1 (MQ=255)
                                                                                                      aaGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaa                                               >  4:347181/1‑149 (MQ=255)
                                                                                                      aaGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaa                                               >  3:480607/1‑149 (MQ=255)
                                                                                                          ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                           <  4:818339/149‑1 (MQ=255)
                                                                                                          ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                           <  1:460231/149‑1 (MQ=255)
                                                                                                            aTGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGc                                         <  1:195661/149‑1 (MQ=255)
                                                                                                                gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACg                                     >  2:385088/1‑149 (MQ=255)
                                                                                                                gATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACg                                     >  2:659013/1‑149 (MQ=255)
                                                                                                                 aTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACgg                                    >  3:480646/1‑149 (MQ=255)
                                                                                                                 aTTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACgg                                    >  2:272892/1‑149 (MQ=255)
                                                                                                                   ttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta                                  <  4:107126/149‑1 (MQ=255)
                                                                                                                    tataTCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtat                                 <  2:457072/149‑1 (MQ=255)
                                                                                                                      tatCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaa                               >  4:847562/1‑149 (MQ=255)
                                                                                                                        tCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAAc                             >  2:198448/1‑149 (MQ=255)
                                                                                                                        tCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAAc                             >  2:768194/1‑149 (MQ=255)
                                                                                                                             tAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTgg                        <  1:410306/149‑1 (MQ=255)
                                                                                                                             tAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTgg                        <  1:198448/149‑1 (MQ=255)
                                                                                                                                gagTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaa                     >  4:473299/1‑149 (MQ=255)
                                                                                                                                  gTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCaaaa                   <  4:480607/149‑1 (MQ=255)
                                                                                                                                   tGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt                  >  2:438595/1‑149 (MQ=255)
                                                                                                                                   tGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt                  <  1:272892/149‑1 (MQ=255)
                                                                                                                                   tGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAAt                  >  4:281586/1‑149 (MQ=255)
                                                                                                                                     cAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATcc                >  3:910880/1‑149 (MQ=255)
                                                                                                                                      aGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCa               >  1:449790/1‑149 (MQ=255)
                                                                                                                                       gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa              <  2:539875/149‑1 (MQ=255)
                                                                                                                                       gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa              <  3:347181/149‑1 (MQ=255)
                                                                                                                                       gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa              <  1:907037/149‑1 (MQ=255)
                                                                                                                                         aCTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                       <  1:803985/74‑1 (MQ=255)
                                                                                                                                         aCTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                       >  2:803985/1‑74 (MQ=255)
                                                                                                                                         aCTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACa            <  1:659013/149‑1 (MQ=255)
                                                                                                                                           tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa          >  4:508909/1‑149 (MQ=255)
                                                                                                                                           tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa          <  1:385088/149‑1 (MQ=255)
                                                                                                                                              tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCt       >  1:613756/1‑149 (MQ=255)
                                                                                                                                                   gATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTaaaaa  >  1:607672/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAA  >  NC_000913/542944‑543239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: