Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,025,579 T→C 8.7% S23S (AGT→AGC) ‡ yedP → putative mannosyl‑3‑phosphoglycerate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,025,5790TC 8.7% 213.8 / 10.1 69S23S (AGT→AGC) ‡yedPputative mannosyl‑3‑phosphoglycerate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base C (3/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GGATCGCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAA  >  NC_000913/2025432‑2025724
                                                                                                                                                   |                                                                                                                                                 
ggATCGCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGtt                                                                                                                                                  <  4:318642/149‑1 (MQ=255)
ggATCGCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGtt                                                                                                                                                  <  4:574309/149‑1 (MQ=255)
  aTCGCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTAt                                                                                                                                                <  1:407863/149‑1 (MQ=255)
     gCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGAc                                                                                                                                             <  3:588017/149‑1 (MQ=255)
       tCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTg                                                                                                                                           >  3:159842/1‑149 (MQ=255)
         ccGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAACTATGACTgtc                                                                                                                                         <  3:686983/149‑3 (MQ=255)
          cGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGCCAGTCATAACTATGACTgtcc                                                                                                                                        <  2:189600/149‑4 (MQ=255)
                 atcatcTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAACTATGACTgtcca                                                                                                                                       >  3:211958/1‑139 (MQ=255)
                 atcatcTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAACTATGACTgtcca                                                                                                                                       <  4:211958/143‑5 (MQ=255)
                  tcatcTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAg                                                                                                                                >  4:764020/1‑149 (MQ=255)
                       tttGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTAt                                                                                                                                                <  3:576097/128‑1 (MQ=255)
                       tgtGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACGACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTAt                                                                                                                                                >  4:576097/3‑128 (MQ=255)
                            tAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCt                                                                                                                      >  2:704880/1‑149 (MQ=255)
                                   tgCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGt                                                                                                               <  4:749696/149‑1 (MQ=255)
                                   tgCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGt                                                                                                               <  1:26963/149‑1 (MQ=255)
                                    gCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGtt                                                                                                              <  1:704880/149‑1 (MQ=255)
                                         tttGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTAcgc                                                                                                         <  3:704281/149‑1 (MQ=255)
                                            gagaATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGaa                                                                                                      >  4:270068/1‑149 (MQ=255)
                                               aaTCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCa                                                                                                   >  4:190489/1‑149 (MQ=255)
                                                atcatcCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCaa                                                                                                  >  3:287437/1‑149 (MQ=255)
                                                 tcatcCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGAc                                                                                                                                             <  1:82877/105‑1 (MQ=255)
                                                 tcatcCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGAc                                                                                                                                             >  2:82877/1‑105 (MQ=255)
                                                  catcCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAAt                                                                                                <  3:389634/149‑1 (MQ=255)
                                                      cTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTc                                                                                            <  3:764020/149‑1 (MQ=255)
                                                        ggTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTccc                                                                                          >  2:203262/1‑149 (MQ=255)
                                                         gTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCg                                                                                         >  1:270819/1‑149 (MQ=255)
                                                           aaGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTc                                                                                       <  4:159842/149‑1 (MQ=255)
                                                            agagGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCa                                                                                      >  4:732664/1‑149 (MQ=255)
                                                             gaggagTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGcccc                                                                                                                            >  3:541924/1‑110 (MQ=255)
                                                             gaggagTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCAt                                                                                     >  3:71112/1‑149 (MQ=255)
                                                             gaggagTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCAt                                                                                     >  3:571035/1‑149 (MQ=255)
                                                             gaggagTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCAGAGTTATGACTGGCAACCGGCAGcccc                                                                                                                            <  4:541924/110‑1 (MQ=255)
                                                                  gTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTctct                                                                                <  4:200581/149‑1 (MQ=255)
                                                                     cTGCGCGACTATGTTTTCAATTCAACAACCCCTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTa                                                                             <  3:111169/149‑1 (MQ=255)
                                                                        cgcgACTATGTTTTCAACTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGAGAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTGTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCa                                                                          <  4:501218/149‑1 (MQ=255)
                                                                             cTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTaaa                                                                     >  4:232629/1‑149 (MQ=255)
                                                                             cTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTaaa                                                                     >  2:519750/1‑149 (MQ=255)
                                                                                  ttttCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATc                                                                >  3:131242/1‑149 (MQ=255)
                                                                                     tcaattcaaCAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGc                                                             <  1:228166/149‑1 (MQ=255)
                                                                                     tcaattcaaCAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGc                                                             <  2:270819/149‑1 (MQ=255)
                                                                                        attcaaCAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGa                                                          >  3:410116/1‑149 (MQ=255)
                                                                                         ttcaaCAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGaa                                                         >  4:576119/1‑149 (MQ=255)
                                                                                             acaacCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGc                                                     <  2:334466/149‑1 (MQ=255)
                                                                                                  cACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTAc                                                >  3:579471/1‑149 (MQ=255)
                                                                                                  cACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTAc                                                >  4:162228/1‑149 (MQ=255)
                                                                                                   actactGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACt                                               >  1:412709/1‑149 (MQ=255)
                                                                                                    ctactGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACtt                                              >  1:163322/1‑149 (MQ=255)
                                                                                                    ctactGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACtt                                              <  3:732664/149‑1 (MQ=255)
                                                                                                    ctactGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACtt                                              <  4:131242/149‑1 (MQ=255)
                                                                                                        tGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGa                                                                                                                                              >  1:161454/1‑49 (MQ=255)
                                                                                                        tGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGa                                                                                                                                              <  2:161454/49‑1 (MQ=255)
                                                                                                          gTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCaaaa                                        >  4:696957/1‑149 (MQ=255)
                                                                                                           tttttAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGc                                                     <  2:460923/135‑1 (MQ=255)
                                                                                                           tttttAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGc                                                     >  1:460923/1‑135 (MQ=255)
                                                                                                            ttttAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCaaaaaa                                      <  4:287437/149‑1 (MQ=255)
                                                                                                                aGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGtt                                  <  3:232629/149‑1 (MQ=255)
                                                                                                                aGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGtt                                  <  1:519750/149‑1 (MQ=255)
                                                                                                                    aTCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTgggg                              >  4:623109/1‑149 (MQ=255)
                                                                                                                     tCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGc                             >  2:301063/1‑149 (MQ=255)
                                                                                                                              cACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGtt                    <  4:410116/149‑1 (MQ=255)
                                                                                                                                    gctgGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTa                                                 <  2:456210/114‑1 (MQ=255)
                                                                                                                                    gctgGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTa                                                 >  1:456210/1‑114 (MQ=255)
                                                                                                                                        gccAGTCATAACTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGAt          >  1:707780/3‑149 (MQ=255)
                                                                                                                                         ccAGTCATAACTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGAtt         >  1:791989/2‑149 (MQ=255)
                                                                                                                                         aCAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGAtt         <  1:203262/149‑1 (MQ=255)
                                                                                                                                             tCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAACCATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCg     >  2:582629/1‑149 (MQ=255)
                                                                                                                                             tCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCg     >  4:507149/1‑149 (MQ=255)
                                                                                                                                               aTAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAgg                      <  1:802949/130‑1 (MQ=255)
                                                                                                                                               aTAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAgg                      >  2:802949/1‑130 (MQ=255)
                                                                                                                                               aTAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGaa   >  3:60688/1‑149 (MQ=255)
                                                                                                                                                tAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGaaa  >  3:516248/1‑149 (MQ=255)
                                                                                                                                                   ttATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACaa                        <  4:238642/124‑1 (MQ=255)
                                                                                                                                                   ttATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACaa                        <  4:411014/124‑1 (MQ=255)
                                                                                                                                                   ttATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACaa                        >  3:238642/1‑124 (MQ=255)
                                                                                                                                                   ttATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACaa                        >  3:411014/1‑124 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGATCGCTCCCGAAAAAATCATCTTTGATAATTTGTGCGGTTTTGAGAATCATCCTGGTAAGAGGAGTACTGCGCGACTATGTTTTCAATTCAACAACCACTACTGGTTTTTAGCGATCTTGATGGCACCCTGCTGGACAGTCATAGTTATGACTGGCAACCGGCAGCCCCCTGGCTCACCCGTTTACGCGAAGCAAATGTTCCCGTCATTCTCTGTAGCAGTAAAACATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAA  >  NC_000913/2025432‑2025724

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