Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 639,560 G→C 11.1% intergenic (+52/‑193) ahpC → / → ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913639,5600GC11.1% 1630.6 / 43.5 514intergenic (+52/‑193)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (144/304);  new base C (64/2);  total (208/306)
Fisher's exact test for biased strand distribution p-value = 8.26e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 528 total aligned reads displayed.

AGTACGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGC  >  NC_000913/639411‑639709
                                                                                                                                                     |                                                                                                                                                     
aGTACTTAGCTTCTCACCCAGGTGAAGTTTGCCCGGATAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCg                                                                                                                                                       <  2:2971665/150‑1 (MQ=255)
    cGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCa                                                                                                                                                   <  2:4726650/150‑1 (MQ=255)
     gTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCac                                                                                                                                                  >  2:250460/1‑150 (MQ=255)
        gCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCact                                                                                                                                                 >  1:1138306/1‑148 (MQ=255)
          ttCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttact                                                                                                                                             <  1:7782840/150‑1 (MQ=255)
          ttCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttact                                                                                                                                             <  1:727916/150‑1 (MQ=255)
          ttCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttact                                                                                                                                             >  2:6621118/1‑150 (MQ=255)
             tcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtgt                                                                                                                                          <  2:7853446/150‑1 (MQ=255)
             tcACCCAGGTGAAGTTGGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtgt                                                                                                                                          <  1:6986186/150‑1 (MQ=255)
               aCCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACATACTTgtgt                                                                                                                                        >  2:3164423/1‑150 (MQ=255)
                cccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCCGTCACATACTATTGTa                                                                                                                                       >  1:1743635/1‑150 (MQ=255)
                 ccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAATTCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTTGCGTGGGCCGCTCACCCCGGTCACTTAATTGTGTaa                                                                                                                                      >  1:2676525/1‑150 (MQ=255)
                 ccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTaa                                                                                                                                      >  1:4703597/1‑150 (MQ=255)
                   aGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGc                                                                                                                                    <  2:6295157/150‑1 (MQ=255)
                        aaGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCg                                                                                                                               <  2:3822793/150‑1 (MQ=255)
                        aaGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCg                                                                                                                               <  1:4747963/150‑1 (MQ=255)
                         aGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCgg                                                                                                                              <  2:8226290/150‑1 (MQ=255)
                         aGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCgg                                                                                                                              <  1:2801753/150‑1 (MQ=255)
                          gTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGGGTCGGCCGCTCACCCCGGTGACTTTCTTGTGTTATCTCCCggg                                                                                                                             >  2:3991408/1‑150 (MQ=255)
                          gTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGGTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCGGTCACTTACTTGTGTAAGCTCCCggg                                                                                                                             >  2:5335860/1‑150 (MQ=255)
                           tttGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCACCCCCGGCACATACGTTTGTAATCTCCCcggg                                                                                                                            >  1:43282/1‑150 (MQ=255)
                                ccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTca                                                                                                                       >  2:2780968/1‑150 (MQ=255)
                                ccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGGGCGGCGCCCCCGCACCCCCGTCCCTTACGTGTGTAAGCGCCCCGGGATTca                                                                                                                       >  2:7649278/1‑150 (MQ=255)
                                 cGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcac                                                                                                                      <  2:3487480/150‑1 (MQ=255)
                                  ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                                                                                                                     <  1:5670345/150‑1 (MQ=255)
                                  ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGGCACTTAAATGTGTAAGCTCCCGGGGATTaaca                                                                                                                     >  1:5910000/1‑150 (MQ=255)
                                   gCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAg                                                                                                                    <  1:6349876/150‑1 (MQ=255)
                                   gCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAg                                                                                                                    <  1:1258396/150‑1 (MQ=255)
                                   gCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAg                                                                                                                    <  1:2477089/150‑1 (MQ=255)
                                     tAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGc                                                                                                                  <  1:364550/150‑1 (MQ=255)
                                     tAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGc                                                                                                                  <  1:4554306/150‑1 (MQ=255)
                                          ggAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGGCACATAAATGTGTAAGATCCCCGGGATTCACAAGCTAgcc                                                                                                             >  1:6482966/1‑150 (MQ=255)
                                          ggAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCTCACCCCCGTCACATACTTTTGTTATCTCCCGGGGGTTGACAGGCTAgcc                                                                                                             >  1:7686475/1‑150 (MQ=255)
                                          ggAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGGGCGGCGCCCCCGCACCCCCGGCCCATACGTGTGTAAGCGCCCCGGGCGTCACCGGCGAgcc                                                                                                             >  1:4972188/1‑150 (MQ=255)
                                             aagaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGGCGCCCCCGCACCCCCGGCACGTACTTGTGTTAGCTCCCGGGGAGTCACAGGCTCgccgcc                                                                                                          >  1:86002/1‑150 (MQ=255)
                                              agaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCt                                                                                                         <  2:5046618/150‑1 (MQ=255)
                                              agaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACATGGGTAAGCTCCCGGGGAATCACAGGCTAGCCGCCt                                                                                                         >  2:5151972/1‑150 (MQ=255)
                                               gaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCtt                                                                                                        <  2:1948871/150‑1 (MQ=255)
                                                 agGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGc                                                                                                      <  1:3531996/150‑1 (MQ=255)
                                                  ggTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGct                                                                                                     <  1:3229773/150‑1 (MQ=255)
                                                    tGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctct                                                                                                   <  2:4529330/150‑1 (MQ=255)
                                                      aaGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGa                                                                                                 <  2:6375584/150‑1 (MQ=255)
                                                       aGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGAc                                                                                                <  2:1659547/150‑1 (MQ=255)
                                                       aGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGGGCGTCGCCCGCACCCCCCCGGCACGTACATTTGTAAGCACCCCGGGCGTGACAGGCTAGCCGCCGTGcgcgggc                                                                                                >  1:803222/1‑144 (MQ=255)
                                                        gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGGGTCGCCCCCACACCCCCGTCACGTAATTTTGTAAGCACCCCGGGGTTCACAGGCTAGCCGCCGTGGTATGaag                                                                                               >  1:1039927/1‑148 (MQ=255)
                                                        gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGAcg                                                                                               <  1:6699763/150‑1 (MQ=255)
                                                        gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCCGTCACATACTTGTGTTAGCGCCCCGGGATTCACAGGCTAGCCCCCTTTCTCTGaag                                                                                               >  2:2125676/1‑148 (MQ=255)
                                                        gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTGGGCGTCGCCCGCGCACCCCCGGCACGAACGTGTGTAAGCGCCCCGGGATTCACAAGCTAGCCGCCATGCGCAGAcc                                                                                               >  1:760537/1‑149 (MQ=255)
                                                            ctctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGaa                                                                                           <  2:2903264/150‑1 (MQ=255)
                                                                 gCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACt                                                                                      <  1:6015664/150‑1 (MQ=255)
                                                                  cTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACtt                                                                                     <  2:5683221/150‑1 (MQ=255)
                                                                  cTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACtt                                                                                     <  2:5518624/150‑1 (MQ=255)
                                                                      gTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCGCACCCCCCGCACATAATTGTGTAAACTCCCCGGGGTTCACAAGCTAGCCGCCGTGCTCAGACCCCAAATaaatagga                                                                                 >  1:5290356/1‑143 (MQ=255)
                                                                       tctcTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCCCCCGCGCACCCCCGTCACTTACATGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACATCGGaa                                                                                >  2:7758921/1‑150 (MQ=255)
                                                                           tGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACATACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACTTCGGAAATTc                                                                            >  2:584909/1‑150 (MQ=255)
                                                                             gACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCAc                                                                          <  1:5355770/150‑1 (MQ=255)
                                                                                 tggttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCCGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACGTCGGAAATTCACCTaa                                                                      >  2:431739/1‑150 (MQ=255)
                                                                                 tggttggtAAAATCTAAATTTCCTTCGTCTTTCACCCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGACACTTACTTCTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTaa                                                                      <  1:2322454/150‑1 (MQ=255)
                                                                                   gttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTAACCCCCGTCACTTACTTGTGTAAGCTCCCGGGGAGTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTGCGGAAATTCACCTAAtt                                                                    >  2:1160009/1‑150 (MQ=255)
                                                                                   gttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCCCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACTTCGGAAAATCACCGAAtt                                                                    >  2:487714/1‑150 (MQ=255)
                                                                                      ggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttctt                                                                 <  2:2461436/150‑1 (MQ=255)
                                                                                       gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                <  2:2657103/150‑1 (MQ=255)
                                                                                       gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                <  1:3462362/150‑1 (MQ=255)
                                                                                       gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                <  1:2037840/150‑1 (MQ=255)
                                                                                       gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                <  2:1404449/150‑1 (MQ=255)
                                                                                       gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAAttcttc                                                                <  1:3995630/150‑1 (MQ=255)
                                                                                        tAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCCCTTACTTGTGTAAGCTCCCGGGGCTTCACAGGCTAGCCGCCTGGCTCGGACGCGAAATACTTCGGAAATTCACCGAATTCTTCg                                                               >  2:3316998/1‑150 (MQ=255)
                                                                                        tAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCg                                                               <  1:1937252/150‑1 (MQ=255)
                                                                                         aaaaTCTAAATTACCTTCGTCTTTCACTCCATAGCGGCGTTGGCGTCGCCCGCACACCCCCGTCACTTACTAGTGTAAGCACCCGGGGATTCACACGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCgg                                                              <  2:5935045/150‑1 (MQ=255)
                                                                                               tAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGtgctg                                                        <  2:3024565/150‑1 (MQ=255)
                                                                                                aaaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGtgctgc                                                       <  2:1025341/150‑1 (MQ=255)
                                                                                                aaaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCGCCCCCCGGTCAAAAAAATGTGTAAGCACCCGGGGATTCACAGGCTAGCCCCCGTGCTCGGACGCGAAAAACTACGGAAAATAACCTAAATAGTCGGGtgctgc                                                       >  2:3713682/1‑150 (MQ=255)
                                                                                                 aaTTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCg                                                      <  1:6827804/150‑1 (MQ=255)
                                                                                                  aggtCCTTCGTCGTTCACGCCAGCGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTCGTGGAAGCTCCCGGGGAGTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCgg                                                     <  1:7401137/147‑1 (MQ=255)
                                                                                                   tttCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACGTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGAAGCGAAATACTTCGGAAATTAACCGAATTCTTCGGGTGCTGCGGc                                                    >  2:2242372/1‑150 (MQ=255)
                                                                                                     tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCAc                                                  <  2:3670156/150‑1 (MQ=255)
                                                                                                       cTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCAACCCCGTCACTTACTTGTGTAAGCTCCAGGGGATTCACATGCTAGCCCCCTTGCACAGAAGCGAAATACTTCGGAAAATCACCGAATTCGTCGGGGTCTGCGGCcccc                                                >  2:7144113/1‑150 (MQ=255)
                                                                                                        ttCGTCTTTCACGCCATAGCGGCGTTGGCGTAGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCCGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGACATTCACCAAATTCTTAGGGTGCTGCGGCACCCg                                               >  2:4403306/1‑150 (MQ=255)
                                                                                                         tCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGa                                              <  1:5120468/150‑1 (MQ=255)
                                                                                                              tttCACGCCATCGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCGCCGGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTGGGGTGCTGCGGCACCCGATTttc                                         <  1:1544102/150‑1 (MQ=255)
                                                                                                              tttCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTttc                                         <  2:4111676/150‑1 (MQ=255)
                                                                                                                  aCGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTcc                                     <  1:6174755/150‑1 (MQ=255)
                                                                                                                   cGCCATAGCGGCGTTGGCGTCGCCCCCACACCCCGGTCACTTACTTGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCGTTCGCAGACGCGAAAAACTTCGGAAATTCACCAAAATATTCGGGGGCGGCGGCCCCCGAGTTTATTccc                                    >  2:1694859/1‑150 (MQ=255)
                                                                                                                      cATAGCGGCGTTGGCGTCGCCCCCTCACCCCGGTCACTTAATTGTGTAAGCGCCCGGGGATTCACAGGCTAGCCGCCGTGCTCTGACGCGAAATACTTCGGAAATTAACCTAATTATTCGGGTGCGGCGGCACCCGATTTTATTACCCGc                                 >  2:85074/1‑150 (MQ=255)
                                                                                                                        tAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCAc                               <  1:3345534/150‑1 (MQ=255)
                                                                                                                        tAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCAc                               <  2:1483572/150‑1 (MQ=255)
                                                                                                                        tAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCAc                               <  2:1455508/150‑1 (MQ=255)
                                                                                                                           cggcgTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCat                            >  2:6821145/1‑150 (MQ=255)
                                                                                                                                ttGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGCTTCACAGGCTAGCCGCCTTGCTCTGACGCGAAAAACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGTTGc                       >  1:2993261/1‑150 (MQ=255)
                                                                                                                                 tGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCa                      <  1:3843542/150‑1 (MQ=255)
                                                                                                                                   gCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAg                    <  2:5779534/150‑1 (MQ=255)
                                                                                                                                    cGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGc                   <  1:1381346/150‑1 (MQ=255)
                                                                                                                                    cGTCGCCCGCTCACCCCGGTCACTTACTGGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCGTCTGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGc                   <  2:7342546/150‑1 (MQ=255)
                                                                                                                                     gTCGCCCGCTCACCCCGGTCACTTACTTGTGTCAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGAGTTTCGTCCCCGCACCATGATGCAAGCt                  >  2:7475744/1‑150 (MQ=255)
                                                                                                                                             cTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAg          <  2:1867045/150‑1 (MQ=255)
                                                                                                                                                  cccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGc     <  1:4044599/150‑1 (MQ=255)
                                                                                                                                                  cccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGc     >  2:4169632/1‑150 (MQ=255)
                                                                                                                                                  cccGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGc     >  2:8097121/1‑150 (MQ=255)
                                                                                                                                                    cGGTCACTTACTGGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTCCTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCg   <  2:7982598/150‑1 (MQ=255)
                                                                                                                                                     ggTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGc  <  2:6115593/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTACGTAGCTTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGC  >  NC_000913/639411‑639709

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