Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,895,060 G→C 6.2% V86L (GTG→CTG)  pabB → aminodeoxychorismate synthase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,0600GC 6.2% 275.1 / 6.6 81V86L (GTG→CTG) pabBaminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (33/43);  new base C (2/3);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGC  >  NC_000913/1894913‑1895208
                                                                                                                                                   |                                                                                                                                                    
ggCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTg                                                                                                                                                    <  2:972604/150‑1 (MQ=255)
   tATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTg                                                                                                                                                 <  2:240966/150‑1 (MQ=255)
    aTGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTANAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGCCCCGCTACAGGTGCTCCAGCAGGTGCTgg                                                                                                                                                <  2:538901/150‑1 (MQ=255)
     tGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGa                                                                                                                                               >  1:837055/1‑150 (MQ=255)
       ccGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATc                                                                                                                                             >  1:472992/1‑150 (MQ=255)
       ccGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATc                                                                                                                                             <  1:458196/150‑1 (MQ=255)
         gATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATcgc                                                                                                                                           >  1:526359/1‑150 (MQ=255)
         gATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagcac                                                                                                                                           <  1:397446/150‑5 (MQ=255)
          atcatcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagcacc                                                                                                                                          <  2:524809/150‑6 (MQ=255)
           tcatcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCACCTGCTGGagcacct                                                                                                                                         <  2:487466/150‑7 (MQ=255)
             atcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCTCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAg                                                                                                                                       <  2:803331/150‑1 (MQ=255)
             atcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAg                                                                                                                                       <  1:361462/150‑1 (MQ=255)
             atcCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAg                                                                                                                                       <  2:1335530/150‑1 (MQ=255)
               ccGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACTACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGAc                                                                                                                                     <  2:553518/150‑1 (MQ=255)
               ccGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCAGTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGAc                                                                                                                                     <  2:558108/150‑1 (MQ=255)
                cGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACa                                                                                                                                    >  2:891446/1‑150 (MQ=255)
                             tGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCa                                                                                                                       >  2:781515/1‑150 (MQ=255)
                             tGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCa                                                                                                                       <  1:781515/150‑1 (MQ=255)
                               atatTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATa                                                                                                                     >  2:531544/1‑150 (MQ=255)
                               atatTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATa                                                                                                                     <  1:1253541/150‑1 (MQ=255)
                                      ggTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGa                                                                                                              <  1:6904/150‑1 (MQ=255)
                                       gTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGAt                                                                                                             <  2:526359/150‑1 (MQ=255)
                                        tCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGAtt                                                                                                            >  1:756672/1‑150 (MQ=255)
                                         cGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGAttt                                                                                                           >  1:549278/1‑150 (MQ=255)
                                          gccgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTg                                                                                                          >  1:869628/1‑150 (MQ=255)
                                           ccgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGc                                                                                                         <  1:49554/150‑1 (MQ=255)
                                            cgagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGcc                                                                                                        <  2:1300974/150‑1 (MQ=255)
                                             gagccgaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCa                                                                                                       >  1:363651/1‑150 (MQ=255)
                                                  gaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCa                                                                                                  >  2:595438/1‑150 (MQ=255)
                                                  gaTTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCa                                                                                                  >  2:42769/1‑150 (MQ=255)
                                                    tttGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGCCATTCGCCCAACGCATAACGAAGATTTGCCATTTCAgg                                                                                                >  2:1091581/1‑150 (MQ=255)
                                                      tGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGggc                                                                                              <  2:822506/150‑1 (MQ=255)
                                                               aCCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTg                                                                                     >  1:289812/1‑150 (MQ=255)
                                                                ccACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgg                                                                                    >  1:346252/1‑150 (MQ=255)
                                                                        ggTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGttt                                                                            >  1:172803/1‑150 (MQ=255)
                                                                         gTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTg                                                                                     <  1:826841/140‑1 (MQ=255)
                                                                         gTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTg                                                                                     >  2:826841/1‑140 (MQ=255)
                                                                         gTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTg                                                                           <  1:380636/150‑1 (MQ=255)
                                                                           aaaGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGCCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGCTCGCGCAGACATTCGCCCACCGCATAACGAAGATTTGCCCTTTCAGGGCGGCGCACTGGGGTTGTTTGGc                                                                         >  2:470761/1‑150 (MQ=255)
                                                                            aaGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGGTCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                                                        <  1:884458/150‑1 (MQ=255)
                                                                            aaGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                                                        <  2:346252/150‑1 (MQ=255)
                                                                            aaGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCt                                                                        <  1:807873/150‑1 (MQ=255)
                                                                             aGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTa                                                                       >  2:802125/1‑150 (MQ=255)
                                                                               aatccGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACg                                                                     <  2:149757/147‑1 (MQ=255)
                                                                               aaaCCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACg                                                                     <  2:154902/150‑1 (MQ=255)
                                                                                aaCCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGa                                                                    >  1:329653/1‑150 (MQ=255)
                                                                                      tgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                              <  1:508563/150‑1 (MQ=255)
                                                                                      tgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                              <  1:483777/150‑1 (MQ=255)
                                                                                           gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccgcc                                                         >  2:832725/1‑150 (MQ=255)
                                                                                           gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccgcc                                                         >  1:409642/1‑150 (MQ=255)
                                                                                             gAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGt                                                       <  1:1045186/150‑1 (MQ=255)
                                                                                                 gCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATCACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   <  2:661557/150‑1 (MQ=255)
                                                                                                 gCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   <  2:680390/150‑1 (MQ=255)
                                                                                                 gCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   >  2:497514/1‑150 (MQ=255)
                                                                                                 gCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   >  1:582353/1‑150 (MQ=255)
                                                                                                   gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  2:582353/150‑1 (MQ=255)
                                                                                                   gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  2:869628/150‑1 (MQ=255)
                                                                                                   gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  1:595438/150‑1 (MQ=255)
                                                                                                   gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                 <  1:891446/150‑1 (MQ=255)
                                                                                                        aCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                            >  1:683293/1‑150 (MQ=255)
                                                                                                         cGCACAACTACCAATGATGACCAGCTACAGGTGCTACAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTg                                           >  2:380183/1‑150 (MQ=255)
                                                                                                         cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCa                                              <  1:659073/147‑1 (MQ=255)
                                                                                                         cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCa                                              >  2:659073/1‑147 (MQ=255)
                                                                                                         cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTg                                           <  2:756672/150‑1 (MQ=255)
                                                                                                                  aCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAtt                                  <  2:472992/150‑1 (MQ=255)
                                                                                                                    cACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGc                                <  1:789202/150‑1 (MQ=255)
                                                                                                                      cTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGTCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCgg                              <  1:787009/150‑1 (MQ=255)
                                                                                                                          tgaCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACa                          <  2:837055/150‑1 (MQ=255)
                                                                                                                           gaCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACaa                         >  1:655464/1‑150 (MQ=255)
                                                                                                                            aCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   >  1:844899/1‑123 (MQ=255)
                                                                                                                            aCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                   <  2:844899/123‑1 (MQ=255)
                                                                                                                            aCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAg                        >  1:637166/1‑150 (MQ=255)
                                                                                                                              ccGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGat                      <  1:204070/150‑1 (MQ=255)
                                                                                                                               cGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                     >  2:1145480/1‑150 (MQ=255)
                                                                                                                                 cTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATc                   <  2:329653/150‑1 (MQ=255)
                                                                                                                                        cgatCCAGCACCTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGc            >  1:871152/4‑150 (MQ=255)
                                                                                                                                        cgatCCAGCACCTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGc            >  2:582618/4‑150 (MQ=255)
                                                                                                                                              agcagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGata      >  1:431417/1‑150 (MQ=255)
                                                                                                                                                cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg    <  1:1307627/150‑1 (MQ=255)
                                                                                                                                                cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg    <  1:34427/150‑1 (MQ=255)
                                                                                                                                                cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg    >  2:1120169/1‑150 (MQ=255)
                                                                                                                                                  ggTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGc  <  2:1395944/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGC  >  NC_000913/1894913‑1895208

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: