Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,059 G→C 5.2% Q85H (CAG→CAC pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,0590GC 5.2% 371.1 / 9.8 115Q85H (CAG→CACpabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base G (56/53);  new base C (2/4);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

Only 100 of 123 total aligned reads displayed.

TCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGC  >  NC_000913/1894910‑1895208
                                                                                                                                                     |                                                                                                                                                     
agcGGCTATACGGATGATCGGCATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCTGTAAAGAAAACGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCcagcag                                                                                                                                                       >  1:1029783/3‑150 (MQ=255)
 ccGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTTAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCggg                                                                                                                                                      >  2:366514/1‑150 (MQ=255)
 ccGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAgg                                                                                                                                                      <  2:1200668/150‑1 (MQ=255)
     cTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCt                                                                                                                                                  >  1:472836/1‑150 (MQ=255)
      tATGCCGATCATCCGTCTAGCCGCTTTTATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAATAAACCGTTGTTAGTGAAAGCGAAAATCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTg                                                                                                                                                 <  1:602936/150‑1 (MQ=255)
         gCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGag                                                                                                                                              <  1:1868229/150‑2 (MQ=255)
         gCCGATCATCCGTATAGCCGCTATGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTCTTAGTGAAAGCGAAAAACGCACAACGACCACTCATGACCCGCTACAGCTGCTCCAGCAGGTGCTGGaa                                                                                                                                              >  2:1654984/1‑149 (MQ=255)
          ccGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGGTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATc                                                                                                                                             >  2:1922544/1‑150 (MQ=255)
          ccGATCATCCGTATAGCCGCTTTGATACTGTGCGCGCCGAGCCGATTTGCACTTTAACCACTATCGGTACAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagc                                                                                                                                             <  1:81333/150‑3 (MQ=255)
           cGATCATCCGTATAGCCGCTTTGATAATCTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTATAGAAACCGTTGTTAGTGAAAGCGAAATACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagca                                                                                                                                            <  1:109082/150‑4 (MQ=255)
            gATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagcac                                                                                                                                           <  2:104114/150‑5 (MQ=255)
              tcatcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCgcgc                                                                                                                                         <  1:991490/150‑1 (MQ=255)
               catcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCa                                                                                                                                        <  2:1701963/150‑1 (MQ=255)
               catcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCa                                                                                                                                        <  2:1186998/150‑1 (MQ=255)
                   cGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCTGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACTACCACTGAAGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCTCGCAGACa                                                                                                                                    >  2:872387/1‑150 (MQ=255)
                   cGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACa                                                                                                                                    <  1:362516/150‑1 (MQ=255)
                          ccGCTTTGATATTGTGGTCGCCGAGCCTATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGccc                                                                                                                             <  1:147103/150‑1 (MQ=255)
                             cTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAAc                                                                                                                          <  1:829022/150‑1 (MQ=255)
                              tttGATATTGTGGTCGCCGAGCCTATTTGCACTTTAACCACTTTCGGTAAAGAAAACGTTGTTAGTGAAAGCGAAAAACTCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACg                                                                                                                         <  1:1296204/150‑1 (MQ=255)
                                tGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCa                                                                                                                       >  2:1807506/1‑150 (MQ=255)
                                tGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCa                                                                                                                       <  2:1023775/150‑1 (MQ=255)
                                 gATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCAt                                                                                                                      <  1:245911/150‑1 (MQ=255)
                                      tgtgGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAGCCGTTGTTAGTGAAAGCGAAAAACGCCCAACGACCACTGATGCCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCAAACGCTTAACgc                                                                                                                 >  2:169447/1‑149 (MQ=255)
                                         ggTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGATAAACGCTCAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCAAACGCATAACGAAgt                                                                                                              >  2:357272/1‑149 (MQ=255)
                                         ggTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGCCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGa                                                                                                              >  1:501037/1‑150 (MQ=255)
                                          gTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGAt                                                                                                             <  1:1802592/150‑1 (MQ=255)
                                            cGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGAttt                                                                                                           >  2:1912850/1‑150 (MQ=255)
                                              ccgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGc                                                                                                         <  1:1536654/150‑1 (MQ=255)
                                               cgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGcc                                                                                                        <  1:1362007/150‑1 (MQ=255)
                                                  gcctaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCAtt                                                                                                     <  1:878532/150‑1 (MQ=255)
                                                  gccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCAtt                                                                                                     <  1:870985/150‑1 (MQ=255)
                                                  gccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCAtt                                                                                                     >  1:817884/1‑150 (MQ=255)
                                                   ccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCAttt                                                                                                    >  2:1709796/1‑150 (MQ=255)
                                                    cgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCTAAACACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTc                                                                                                   <  1:1302042/150‑1 (MQ=255)
                                                     gaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCa                                                                                                  <  1:514837/150‑1 (MQ=255)
                                                        ttGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCTGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCAAAACGAAGATTTGCCATTTCAggg                                                                                               <  2:947110/150‑1 (MQ=255)
                                                        ttGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAggg                                                                                               >  1:472034/1‑150 (MQ=255)
                                                           cACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGCCATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGggcgg                                                                                            >  1:1036155/1‑150 (MQ=255)
                                                            aCTTTAACCACTATCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGggcggc                                                                                           <  2:93659/150‑1 (MQ=255)
                                                               ttAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGgcgcc                                                                                        >  1:103904/1‑149 (MQ=255)
                                                                 aaCCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACt                                                                                      >  1:924674/1‑150 (MQ=255)
                                                                 aaCCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACt                                                                                      >  1:1125943/1‑150 (MQ=255)
                                                                   ccACTTTGCGTAAAGAAACCGTTGTTAGTGAAAGCCAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGATCCAGCAGGGGCTGGAGCGCGCAGACATTCGCCGAACGCATAACGAATATTTGCGATTTCAGGGCGGCGCACTgg                                                                                    <  2:111419/150‑1 (MQ=255)
                                                                     aCTTTCGGTAATGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCTGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCTGGGCGGCGCACTgggg                                                                                  <  1:1736550/150‑1 (MQ=255)
                                                                      cTTTCGTTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCTGGGCGGCGCACTGGGgt                                                                                 <  2:1117390/150‑1 (MQ=255)
                                                                       tttCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGgtt                                                                                <  2:1616309/150‑1 (MQ=255)
                                                                        ttCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACTCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGgttg                                                                               <  2:1091793/150‑1 (MQ=255)
                                                                           ggTAAAGAAACCGTTGTTAGTGAAAGCCAAAACCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTCGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGAACTGGGGGTGttt                                                                            >  2:1271746/1‑150 (MQ=255)
                                                                            gTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTg                                                                           <  1:205071/150‑1 (MQ=255)
                                                                             tAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTgg                                                                          >  2:438351/1‑150 (MQ=255)
                                                                               aaGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                                                        >  2:402024/1‑150 (MQ=255)
                                                                               aaGAAAACGTTGTTAGTTAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGTATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                                                        <  2:1812102/150‑1 (MQ=255)
                                                                                 taaaCCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTAc                                                                      <  1:547200/149‑1 (MQ=255)
                                                                                      cGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCt                                                                 >  2:1143012/1‑150 (MQ=255)
                                                                                      cGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCt                                                                 >  2:718169/1‑150 (MQ=255)
                                                                                      cGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGAATAACGAAGATTTGCCATTTCAGGGAGGCTCAATGTGGTTGTTTGGCTACGATCt                                                                 >  2:1249175/1‑150 (MQ=255)
                                                                                         tgttAGTGAAAGCGAAAAACGCACAACGCCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCACCGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                              >  2:1085433/1‑150 (MQ=255)
                                                                                         tgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                              <  2:1310497/150‑1 (MQ=255)
                                                                                            tAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGGGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccg                                                           <  1:1063543/150‑1 (MQ=255)
                                                                                              gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGCCATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCAGTGGGGTTGTTTGGCTACGATCTGGgccgcc                                                         >  1:387646/1‑150 (MQ=255)
                                                                                              gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATATGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccgcc                                                         >  1:34077/1‑150 (MQ=255)
                                                                                               tGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGccgccg                                                        >  1:1612281/1‑150 (MQ=255)
                                                                                                 aaaGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtt                                                      >  1:1748718/1‑150 (MQ=255)
                                                                                                  aaGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGttt                                                     >  1:588235/1‑150 (MQ=255)
                                                                                                  aaGCGAAAAACGAACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCCttt                                                     >  2:441077/1‑150 (MQ=255)
                                                                                                   aGCGAAAAACGCACAACGCCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCACCGCGTAACGAAGATTTGCCATTTCAGGGCGGCGCCCGGGGGTTGTTTGGGTACGAGCGGGGCCGCGGtttt                                                    >  2:165358/1‑150 (MQ=255)
                                                                                                   aGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCGGGATCGCGCATACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                                    <  2:743059/150‑1 (MQ=255)
                                                                                                      gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 >  2:922639/1‑150 (MQ=255)
                                                                                                      gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  2:472034/150‑1 (MQ=255)
                                                                                                      gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  1:366514/150‑1 (MQ=255)
                                                                                                      gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  2:1463999/150‑1 (MQ=255)
                                                                                                      gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTACAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  2:1029783/150‑1 (MQ=255)
                                                                                                        aaaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTc                                               >  1:497441/1‑150 (MQ=255)
                                                                                                          aaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCAc                                             <  1:1446678/150‑1 (MQ=255)
                                                                                                            cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTg                                           >  2:1184372/1‑150 (MQ=255)
                                                                                                             gCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGc                                          <  2:472836/150‑1 (MQ=255)
                                                                                                                caACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATACCGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAg                                       <  1:880859/150‑1 (MQ=255)
                                                                                                                 aaCGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGGATAACGAAGATTTGCCATTTCAGGGCGGGGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTGACTGCCAGa                                      <  2:674525/150‑1 (MQ=255)
                                                                                                                      ccACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGCCATTCGCCCAACGCATAACGAAGATTTGCCATTTCCTGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTGTGAGTCACTGCCAGAACTTg                                 >  2:304800/1‑150 (MQ=255)
                                                                                                                      ccACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTCCCAGAAATTg                                 >  1:1125481/1‑150 (MQ=255)
                                                                                                                       cACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGCTTTGCCCTTTCAGGGCGGCGCACTGGGCTTGAATGTCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGc                                >  1:383075/1‑150 (MQ=255)
                                                                                                                        ccTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCg                               <  1:1628214/149‑1 (MQ=255)
                                                                                                                        aCTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCg                               <  1:902390/150‑1 (MQ=255)
                                                                                                                         atgatgaCCCGCTTCAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCTCCCAACTCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCgg                              <  2:1176949/149‑1 (MQ=255)
                                                                                                                                         ggTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCt              >  2:584449/1‑150 (MQ=255)
                                                                                                                                          gTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTg             <  1:1912850/150‑1 (MQ=255)
                                                                                                                                           cgatCCAGCACCTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGc            >  1:1419846/4‑150 (MQ=255)
                                                                                                                                           cgatCCAGCACCTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGCGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCAAGAAATTTCGGAACAAGATATCGTTCTGc            >  2:1268831/4‑150 (MQ=255)
                                                                                                                                            gCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGcc           <  1:1709796/150‑1 (MQ=255)
                                                                                                                                            gCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGcc           >  1:1135255/1‑150 (MQ=255)
                                                                                                                                            gCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGcc           >  1:1983662/1‑150 (MQ=255)
                                                                                                                                            gCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGcc           >  2:1326560/1‑150 (MQ=255)
                                                                                                                                             cTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCATCGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGGTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCg          <  1:1082446/150‑1 (MQ=255)
                                                                                                                                              cccAGCAGGTGCTGGATCGCGCAGACATTGGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCGAGAATTTGCGGAACAAGATATCGTTCTGCCgg         <  1:1654984/149‑1 (MQ=255)
                                                                                                                                               ccAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGa        >  1:559254/1‑150 (MQ=255)
                                                                                                                                                cagcagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATTTCGTTCTGCCGGat       <  2:501037/150‑1 (MQ=255)
                                                                                                                                                cagcagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACGGGGGTTGATTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGat       <  1:1292733/150‑1 (MQ=255)
                                                                                                                                                   gagGTGCGGGATCGCGCTGAGATGCGCCAAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCAGTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg    <  2:1287399/149‑1 (MQ=255)
                                                                                                                                                   cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATATCGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg    <  1:1271746/150‑1 (MQ=255)
                                                                                                                                                     ggTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGc  <  1:724096/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGC  >  NC_000913/1894910‑1895208

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