Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,353 A→T 6.7% L174Q (CTG→CAG)  gatD ← galactitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,3530AT 6.7% 459.4 / 16.5 150L174Q (CTG→CAG) gatDgalactitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (58/82);  new base T (6/4);  total (64/86)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 156 total aligned reads displayed.

AACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTC  >  NC_000913/2172205‑2172500
                                                                                                                                                    |                                                                                                                                                   
      ctGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAAt                                                                                                                                              <  1:458225/150‑1 (MQ=255)
        gCATTTGCGGCGCGTTCATTTCTCTACTGTTATATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGATCTAATGTCGATCGGCGTCACACTCTTTGCTCCCAGGGCTACAGCGGACTGAATGGCCAGCAGGcctgtgg                                                                                                                                            <  2:1395398/150‑6 (MQ=255)
        gCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattca                                                                                                                                            <  1:1891599/150‑5 (MQ=255)
        gCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                                                                             <  2:1143000/149‑4 (MQ=255)
           tttGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTc                                                                                                                                         <  1:1140406/150‑1 (MQ=255)
              gcggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCgg                                                                                                                                      <  2:873216/150‑1 (MQ=255)
               cggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGc                                                                                                                                     <  1:2375205/150‑1 (MQ=255)
               cggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGc                                                                                                                                     <  2:42540/150‑1 (MQ=255)
                ggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCa                                                                                                                                    >  1:2377934/1‑150 (MQ=255)
                 gcgcgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAc                                                                                                                                   <  2:1704956/150‑1 (MQ=255)
                  cgcgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCATACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCAcc                                                                                                                                  <  1:1364389/150‑1 (MQ=255)
                    cgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                                                                                <  2:485235/150‑1 (MQ=255)
                    cgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaa                                                                                                                                >  2:360424/1‑150 (MQ=255)
                     gcTCATTTCGCTACTGTTAAATGTTTGCATCGCAGCGAAAGATTTTTCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCTAGAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaat                                                                                                                               <  2:1155456/150‑1 (MQ=255)
                      cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                                                                              >  1:801222/1‑150 (MQ=255)
                      cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaata                                                                                                                              >  2:1749768/1‑150 (MQ=255)
                        cATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCaataat                                                                                                                            <  1:1481898/150‑1 (MQ=255)
                          tttCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACAGAATGGCCAGCAGGCCAATGCTTCCGGCACCaataataa                                                                                                                          >  2:2077360/1‑150 (MQ=255)
                            tCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAataataata                                                                                                                        <  1:55732/150‑1 (MQ=255)
                            tCGCTACTATTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAataataata                                                                                                                        <  2:31920/150‑1 (MQ=255)
                             cGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAAtaataataa                                                                                                                       <  2:1191513/150‑1 (MQ=255)
                              gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                                      >  1:542320/1‑150 (MQ=255)
                               cTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACg                                                                                                                     <  1:1749768/150‑1 (MQ=255)
                                 aCTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtt                                                                                                                   >  1:2120538/1‑150 (MQ=255)
                                  cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  <  1:1292198/150‑1 (MQ=255)
                                  cTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGttt                                                                                                                  <  2:2190375/150‑1 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                                                                                 <  1:717849/150‑1 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                                                                                 <  2:1472929/150‑1 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                                                                                 >  2:2090838/1‑150 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                                                                                 <  2:639550/150‑1 (MQ=255)
                                   tGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGtttt                                                                                                                 >  2:2022503/1‑150 (MQ=255)
                                     ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                               <  1:1433994/150‑1 (MQ=255)
                                     ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                               <  1:1479561/150‑1 (MQ=255)
                                     ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                               <  1:2057250/150‑1 (MQ=255)
                                       aaaTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTAtt                                                                                                             >  1:1815327/1‑150 (MQ=255)
                                        aaTGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtc                                                                                                            >  1:1218943/1‑150 (MQ=255)
                                          tGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtctc                                                                                                          >  2:1337158/1‑150 (MQ=255)
                                           gTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                         >  1:1763316/1‑150 (MQ=255)
                                           gTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                         <  2:1911115/150‑1 (MQ=255)
                                              tGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                                      >  1:673395/1‑150 (MQ=255)
                                               gCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAAc                                                                                                     >  1:387469/1‑150 (MQ=255)
                                                cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAAcc                                                                                                    <  2:377367/150‑1 (MQ=255)
                                                 aTCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCt                                                                                                   >  1:2443931/1‑150 (MQ=255)
                                                 aTCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCt                                                                                                   >  1:2435610/1‑150 (MQ=255)
                                                     cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcg                                                                                               <  1:247211/150‑1 (MQ=255)
                                                       ccGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCt                                                                                             >  1:152062/1‑150 (MQ=255)
                                                                 ttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCCCTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACATTGCGCTAAATGAAAAg                                                                                   >  2:2420563/1‑150 (MQ=255)
                                                                  tgccagtgccagTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGc                                                                                  >  2:1936396/1‑150 (MQ=255)
                                                                     cagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                               >  2:2248377/1‑150 (MQ=255)
                                                                     cagtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                               <  2:51949/150‑1 (MQ=255)
                                                                       gtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCa                                                                             >  1:1442947/1‑150 (MQ=255)
                                                                              tttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                      <  1:1808210/150‑1 (MQ=255)
                                                                               ttttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACg                                                                     <  2:1576216/150‑1 (MQ=255)
                                                                                tttCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACgg                                                                    <  2:2163820/150‑1 (MQ=255)
                                                                                  tCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAGTGAAAAGCATGCAGACCAACGGTa                                                                  >  2:796119/1‑150 (MQ=255)
                                                                                     taaCTAATGTGTATCGCCCTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCACGCCAATGGTTCCGGCACGAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                               <  2:1644594/149‑1 (MQ=255)
                                                                                     gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATc                                                               >  1:886956/1‑150 (MQ=255)
                                                                                          aaTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTc                                                          <  2:673395/150‑1 (MQ=255)
                                                                                           aTGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACAAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCa                                                         <  1:2181209/150‑1 (MQ=255)
                                                                                            tGTCGATCGCCGTCACACTCTTTGCTCCCATCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCaa                                                        <  1:1120093/150‑1 (MQ=255)
                                                                                             gTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAAt                                                       <  1:135032/150‑1 (MQ=255)
                                                                                              tCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATa                                                      <  2:532939/150‑1 (MQ=255)
                                                                                               cGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGACAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaa                                                     <  1:1009307/150‑1 (MQ=255)
                                                                                                gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATaaa                                                    >  1:2135473/1‑150 (MQ=255)
                                                                                                                tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                                    >  2:287384/1‑150 (MQ=255)
                                                                                                                tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                                    <  2:217234/150‑1 (MQ=255)
                                                                                                                tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                                    <  2:1763316/150‑1 (MQ=255)
                                                                                                                 ttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAg                                   <  2:2468733/150‑1 (MQ=255)
                                                                                                                  tGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAgg                                  <  2:152062/150‑1 (MQ=255)
                                                                                                                   gCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGc                                 <  2:726152/150‑1 (MQ=255)
                                                                                                                    cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCa                                >  2:1142460/1‑150 (MQ=255)
                                                                                                                    cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCa                                >  1:1607218/1‑150 (MQ=255)
                                                                                                                    cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCa                                >  2:1512410/1‑150 (MQ=255)
                                                                                                                    cTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCa                                >  1:840346/1‑150 (MQ=255)
                                                                                                                        cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                      >  1:1070243/1‑108 (MQ=255)
                                                                                                                        cAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAc                                                                      <  2:1070243/108‑1 (MQ=255)
                                                                                                                         aGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATcc                           >  2:1334385/1‑150 (MQ=255)
                                                                                                                          gcgcGACCGCCCACTGCATGGCCAGCAGGCCAATGGTTCCTGCAGCATTAATAATAAGGTTTTTATTCTGACGACGTTGCGCTAAATGAAAAGCATGCAGAGCGACGGGATTCGGCTGAATAAACGCCGCATCCTCAATAGGCATATCCg                          >  2:909946/1‑150 (MQ=255)
                                                                                                                            gcgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                        <  1:360424/150‑1 (MQ=255)
                                                                                                                            gcgACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTg                        <  1:327646/150‑1 (MQ=255)
                                                                                                                              gACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTggg                      <  2:801222/150‑1 (MQ=255)
                                                                                                                                cAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTa                    <  2:1218943/150‑1 (MQ=255)
                                                                                                                                  gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTaga                  >  1:117239/1‑150 (MQ=255)
                                                                                                                                   cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag                 <  2:486745/150‑1 (MQ=255)
                                                                                                                                   cgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCATCAATAGGCATATCCGTGGGTagag                 >  1:1655762/1‑150 (MQ=255)
                                                                                                                                   cgcACTGAATGGCAAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTCCGCTAAATGAAAAGCATGCAGACAAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag                 >  2:265671/1‑150 (MQ=255)
                                                                                                                                    gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGc                >  1:710051/1‑150 (MQ=255)
                                                                                                                                      aaccattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaa              >  1:1325332/7‑150 (MQ=255)
                                                                                                                                       cTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             >  1:267920/1‑150 (MQ=255)
                                                                                                                                       accattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCaaa             >  2:2102256/6‑150 (MQ=255)
                                                                                                                                        ccattGGCCTGCTGGCCTATGGTTCCGGCACCAATAATAATAACGTTTTTAGTGTCACAAGCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCTTATCCGTCGGTAGAGCAAAg            >  2:1441332/5‑150 (MQ=255)
                                                                                                                                         taaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa           <  2:2135473/149‑1 (MQ=255)
                                                                                                                                         gAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa           >  2:403303/1‑150 (MQ=255)
                                                                                                                                         cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa           >  2:1091254/4‑150 (MQ=255)
                                                                                                                                          aaTGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAc          <  1:2090838/150‑1 (MQ=255)
                                                                                                                                              gCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAttt      >  1:1250809/1‑150 (MQ=255)
                                                                                                                                              gCCAGCAGGCCAAAGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAttt      <  2:1210153/150‑1 (MQ=255)
                                                                                                                                               ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACAtttt     <  1:2248377/150‑1 (MQ=255)
                                                                                                                                                ctgcagGCAAATGGTTCCGCCACCAATAATAATAACGTTTTTATTCTCAGAACCTTGCGCTAAATGAAAAGCTTGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGGAAAGACAttttt    <  1:2436116/148‑1 (MQ=255)
                                                                                                                                                  ccagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAAGGGTAATCTGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCTTGGGTAGAGCAAAGACATTTTTTc  <  2:972999/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGCTCTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTC  >  NC_000913/2172205‑2172500

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