Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,091 G→C 5.2% 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,0910GC 5.2% 542.6 / 17.1 156intergenic (+58/+170)glxK/allEglycerate kinase II/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (74/73);  new base C (5/3);  total (80/76)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-02

Only 100 of 157 total aligned reads displayed.

AGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAA  >  NC_000913/542946‑543238
                                                                                                                                                 |                                                                                                                                                   
  tGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCg                                                                                                                                                >  3:80348/1‑149 (MQ=255)
      cGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCAt                                                                                                                                            <  1:600999/149‑1 (MQ=255)
        ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa                                                                                                                                          <  4:1129722/149‑1 (MQ=255)
             ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                                                     <  4:851703/149‑1 (MQ=255)
                  aaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTTCTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                                                <  1:121329/149‑1 (MQ=255)
                   aCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcc                                                                                                                               >  3:515880/1‑149 (MQ=255)
                    ccAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcct                                                                                                                              >  4:423407/1‑149 (MQ=255)
                    ccAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcct                                                                                                                              >  3:756115/1‑149 (MQ=255)
                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                                                     >  2:184314/1‑141 (MQ=255)
                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                                                     <  1:184314/141‑1 (MQ=255)
                     cAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt                                                                                                                             >  2:61684/1‑149 (MQ=255)
                      aaTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                            >  4:618961/1‑149 (MQ=255)
                      aaTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                            >  1:196933/1‑149 (MQ=255)
                      aaTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                            <  3:352230/149‑1 (MQ=255)
                       aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcc                                                                                                                           >  2:359728/1‑149 (MQ=255)
                          tctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGa                                                                                                                        <  3:755439/149‑1 (MQ=255)
                            ttCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaa                                                                                                                      >  3:422054/1‑149 (MQ=255)
                                aCAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcg                                                                                                                  <  2:1280924/149‑1 (MQ=255)
                                aCAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcg                                                                                                                  <  1:1029068/149‑1 (MQ=255)
                                 cAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgc                                                                                                                 <  2:1120542/149‑1 (MQ=255)
                                       gcgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAAt                                                                                                           >  1:730021/1‑149 (MQ=255)
                                           aaTATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGAGGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCta                                                                                                       <  3:956117/149‑1 (MQ=255)
                                             tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                                                     <  1:133884/149‑1 (MQ=255)
                                             tatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                                                     >  2:1104100/1‑149 (MQ=255)
                                              atTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtatat                                                                                                    >  4:511776/1‑149 (MQ=255)
                                               ttGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtatata                                                                                                   >  3:272221/1‑149 (MQ=255)
                                                   cTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc                                                                                               >  1:396997/1‑149 (MQ=255)
                                                   cTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc                                                                                               <  4:729621/149‑1 (MQ=255)
                                                    tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATAttcag                                                                                              >  3:1036854/1‑146 (MQ=255)
                                                    tGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                              <  4:175687/149‑1 (MQ=255)
                                                          aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGTTGCGATTTATATCGCATAAGAGTGCAGTCCTCATGCCGTATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                        >  3:917975/1‑149 (MQ=255)
                                                          aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                        >  3:631981/1‑149 (MQ=255)
                                                          aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                                                        <  3:618961/149‑1 (MQ=255)
                                                            tAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                      <  4:416285/149‑1 (MQ=255)
                                                            tAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                                                                      <  1:244839/149‑1 (MQ=255)
                                                              aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                                                    <  2:730021/149‑1 (MQ=255)
                                                              aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                                                    >  1:967839/1‑149 (MQ=255)
                                                                      cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa                                                                            <  2:1127035/149‑1 (MQ=255)
                                                                      cAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa                                                                            <  3:976871/149‑1 (MQ=255)
                                                                       aGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                           >  1:167220/1‑149 (MQ=255)
                                                                         ggAATTGAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                                                         <  4:1380987/149‑1 (MQ=255)
                                                                         ggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATATAGGCCTGATaa                                                                         <  3:560781/149‑1 (MQ=255)
                                                                          gAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                                                         >  4:62590/1‑132 (MQ=255)
                                                                               aaaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaag                                                                   >  2:297215/1‑149 (MQ=255)
                                                                                 aaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                           >  1:396928/1‑139 (MQ=255)
                                                                                 aaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGAt                                                                           <  2:396928/139‑1 (MQ=255)
                                                                                  aaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAgcgc                                                                >  1:1246222/1‑149 (MQ=255)
                                                                                        cccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGc                                                                                <  4:83778/127‑1 (MQ=255)
                                                                                         ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                                                     <  2:1452599/105‑1 (MQ=255)
                                                                                         ccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGc                                                         >  4:1383791/1‑149 (MQ=255)
                                                                                           ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGc                                                                                <  2:65/124‑1 (MQ=255)
                                                                                             aCCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATg                                                     >  2:991587/1‑149 (MQ=255)
                                                                                                tttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc                                                                                                                            <  2:217230/75‑1 (MQ=255)
                                                                                                 ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTAc                                                 <  2:1014992/149‑1 (MQ=255)
                                                                                                 ttAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTAc                                                 <  4:190487/149‑1 (MQ=255)
                                                                                                      tcttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaa                                            >  4:254118/5‑149 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  4:917975/149‑1 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  3:193325/149‑1 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  4:403487/149‑1 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  4:272221/149‑1 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  1:1454848/149‑1 (MQ=255)
                                                                                                        ggATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  3:1107379/149‑1 (MQ=255)
                                                                                                           tGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGcc                                       >  3:173435/1‑149 (MQ=255)
                                                                                                               aTTTATATCGCATAGGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACgg                                   >  3:111549/1‑149 (MQ=255)
                                                                                                                   atatCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGtata                               >  2:1221062/1‑149 (MQ=255)
                                                                                                                     atCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATaaa                             >  3:457473/1‑149 (MQ=255)
                                                                                                                      tCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAAc                            >  2:1027213/1‑149 (MQ=255)
                                                                                                                        gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCa                                                                                              >  3:113506/1‑81 (MQ=255)
                                                                                                                        gCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTCCa                                                                                              <  4:113506/81‑1 (MQ=255)
                                                                                                                         cATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGt                         >  2:309628/1‑149 (MQ=255)
                                                                                                                          aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                        <  1:1027213/149‑1 (MQ=255)
                                                                                                                          aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                        <  3:254118/149‑1 (MQ=255)
                                                                                                                          aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                        <  1:61684/149‑1 (MQ=255)
                                                                                                                          aTAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTg                        <  2:230138/149‑1 (MQ=255)
                                                                                                                           tAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGcc                                       <  1:243424/133‑1 (MQ=255)
                                                                                                                           tAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTTTATACCGTgg                       <  1:1360181/149‑1 (MQ=255)
                                                                                                                            aaGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGc                      <  4:457473/149‑1 (MQ=255)
                                                                                                                            aaGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGATACAAAAAAAGCCACGGTATAAACCGTGGc                      >  3:1313885/1‑149 (MQ=255)
                                                                                                                              gagTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaagc                                                                  <  1:486657/103‑1 (MQ=255)
                                                                                                                                  gCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATc                <  2:1246222/149‑1 (MQ=255)
                                                                                                                                   cAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATcc               <  1:1441933/149‑1 (MQ=255)
                                                                                                                                   cAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATcc               <  2:528749/149‑1 (MQ=255)
                                                                                                                                     gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa             <  3:1061949/149‑1 (MQ=255)
                                                                                                                                     gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa             <  1:1167306/149‑1 (MQ=255)
                                                                                                                                     gTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa             <  1:359728/149‑1 (MQ=255)
                                                                                                                                     gTACTCATGCCGGAGGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCaa             >  1:747295/1‑149 (MQ=255)
                                                                                                                                      tACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                                                         >  1:484704/1‑72 (MQ=255)
                                                                                                                                      tACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                                                         <  2:484704/72‑1 (MQ=255)
                                                                                                                                         tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         >  3:1389409/1‑149 (MQ=255)
                                                                                                                                         tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         >  4:150269/1‑149 (MQ=255)
                                                                                                                                         tCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         <  3:18342/149‑1 (MQ=255)
                                                                                                                                         tCATGCCGCATCCGGCATGAGTACTATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         >  3:823460/1‑149 (MQ=255)
                                                                                                                                         tCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa                                          <  2:74230/116‑1 (MQ=255)
                                                                                                                                         tCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         >  4:638583/1‑149 (MQ=255)
                                                                                                                                           aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         <  2:1228120/147‑1 (MQ=255)
                                                                                                                                           aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATa         >  1:1228120/1‑147 (MQ=255)
                                                                                                                                           aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGc       >  3:518091/1‑149 (MQ=255)
                                                                                                                                           aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGc       >  2:549829/1‑149 (MQ=255)
                                                                                                                                            tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCt      >  2:1376152/1‑149 (MQ=255)
                                                                                                                                             gCCGGATTCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTa     >  2:795797/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAA  >  NC_000913/542946‑543238

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