Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,644,540 A→C 5.4% intergenic (‑236/+111) ydfU ← / ← rem Qin prophage; uncharacterized protein/Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,644,5400AC 5.4% 313.4 / 4.0 105intergenic (‑236/+111)ydfU/remQin prophage; uncharacterized protein/Qin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (36/62);  new base C (4/3);  total (40/65)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CGAGACTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTC  >  NC_000913/1644398‑1644688
                                                                                                                                              |                                                                                                                                                    
cGAGACTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCg                                                                                                                                               >  1:497603/1‑150 (MQ=255)
 gagaCTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTTCTTGCCTGCACGCACGGGCgg                                                                                                                                              <  1:262829/150‑2 (MQ=255)
  agaCTGAGAGCTTTCTTAACCTCGACCCCGAGCCTGCGGTAATAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGaa                                                                                                                                             <  2:50628/150‑1 (MQ=255)
  agaCTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGATTGCGTGAGAATACCGCCCGTTCTTGCCTGCACGCACGGGCggt                                                                                                                                             <  1:361792/150‑3 (MQ=255)
   gaCTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAAt                                                                                                                                            <  1:223982/150‑1 (MQ=255)
     cTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTa                                                                                                                                          >  2:1030025/1‑150 (MQ=255)
       gagaGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATc                                                                                                                                        >  2:1126548/1‑150 (MQ=255)
           gCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGaa                                                                                                                                    >  2:438669/1‑150 (MQ=255)
             tttCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAAtt                                                                                                                                  <  2:635583/150‑1 (MQ=255)
              ttCTTACCCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGTGGGTGCTGGAACAATACAGACTTGAATGCGTGAGTATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTg                                                                                                                                 >  2:355866/1‑150 (MQ=255)
                   aaCCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAata                                                                                                                            <  1:1168269/150‑1 (MQ=255)
                      cTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATAttt                                                                                                                         >  1:435004/1‑150 (MQ=255)
                       tCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATAtttt                                                                                                                        >  1:860343/1‑150 (MQ=255)
                       tCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATAttat                                                                                                                        >  2:1643917/1‑148 (MQ=255)
                        cGAACCCGCGCCTTCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTttg                                                                                                                       <  1:1418210/150‑1 (MQ=255)
                        cGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTttg                                                                                                                       <  1:1814278/150‑1 (MQ=255)
                            cccGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTttgcgtt                                                                                                                   >  1:169874/1‑150 (MQ=255)
                              cgcgCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTttgcgttgc                                                                                                                 >  1:310790/1‑150 (MQ=255)
                               ccgcCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCGGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTttgcgttgcg                                                                                                                <  2:453533/149‑1 (MQ=255)
                                cgcCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTtgcgttgcgt                                                                                                               >  2:452909/1‑150 (MQ=255)
                                 gcCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTg                                                                                                              <  2:865366/150‑1 (MQ=255)
                                 gcCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTg                                                                                                              <  1:69309/150‑1 (MQ=255)
                                     gCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGAGAGAGAATACCGCCCGTGCGTGCAGGCAAGAACGTGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCAt                                                                                                          >  2:729520/1‑150 (MQ=255)
                                        gTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCgg                                                                                                       <  1:1252842/150‑1 (MQ=255)
                                         tAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGt                                                                                                      <  1:1620623/150‑1 (MQ=255)
                                           agtatgaatCAGCCCTTCGGCCTGGTCTGCAGTGCATGGAGGGTGCCGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGttt                                                                                                    <  2:191007/147‑1 (MQ=255)
                                             aatgaatCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTtc                                                                                                  >  2:1504347/1‑150 (MQ=255)
                                                gaatCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTg                                                                                                                                 >  1:304281/1‑116 (MQ=255)
                                                gaatCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTg                                                                                                                                 <  2:304281/116‑1 (MQ=255)
                                                gaatCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTcc                                                                                               <  2:716123/150‑1 (MQ=255)
                                                   tCAGCAATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCATGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGt                                                                                            <  2:758406/150‑1 (MQ=255)
                                                    cAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTg                                                                                           >  2:1083009/1‑150 (MQ=255)
                                                      gCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGc                                                                                         <  2:954938/150‑1 (MQ=255)
                                                          ttCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACgg                                                                                     <  1:228599/150‑1 (MQ=255)
                                                            cGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGgtg                                                                                   <  2:219990/150‑1 (MQ=255)
                                                            cGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGgtg                                                                                   <  2:1647497/150‑1 (MQ=255)
                                                              gCCTGTTCTGCAGTGCATGGAGGGTGCTGGATCCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCTGTTTTCTCCGGTGGCACGGTGtt                                                                                 <  1:263526/150‑1 (MQ=255)
                                                              gCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGtt                                                                                 <  2:1812857/150‑1 (MQ=255)
                                                              gCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGtt                                                                                 <  1:726415/150‑1 (MQ=255)
                                                                cgggtCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGGGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTAc                                                                               <  1:1601509/146‑1 (MQ=255)
                                                                cTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTAc                                                                               >  1:91582/1‑150 (MQ=255)
                                                                cTGTTCTGCAGGGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTAc                                                                               <  1:1085291/150‑1 (MQ=255)
                                                                     cTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGAGAGAATACCGCCCGTGCGTGCAGGAAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGc                                                                          >  1:973769/1‑150 (MQ=255)
                                                                     cTGCAGTGCATGGAGGGTGCGGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGc                                                                          <  1:212358/150‑1 (MQ=255)
                                                                      tGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCg                                                                         >  2:1296826/1‑150 (MQ=255)
                                                                       gCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGATTTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCgg                                                                        <  1:1191150/150‑1 (MQ=255)
                                                                       gCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCgg                                                                        >  1:933198/1‑150 (MQ=255)
                                                                        cAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCggg                                                                       >  1:1426175/1‑150 (MQ=255)
                                                                        cAGTGCATGGAGGGGGCTGGAACCATTCAGACTAGACTGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCggg                                                                       <  1:884955/150‑1 (MQ=255)
                                                                         aGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCgggt                                                                      >  1:361969/1‑149 (MQ=255)
                                                                         aGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGa                                                                      >  1:1100614/1‑150 (MQ=255)
                                                                          gTGCATGGTGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGGGCGGGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTGTTCTCCGGTGGCACGGTGTTACTCAGCggagt                                                                     <  1:1414968/150‑4 (MQ=255)
                                                                            gCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCCGCGGGAGtt                                                                   >  2:751964/1‑150 (MQ=255)
                                                                             cATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTc                                                                  <  1:456188/150‑1 (MQ=255)
                                                                             cATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTc                                                                  >  1:762095/1‑150 (MQ=255)
                                                                                ggAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTAAGc                                                               >  1:1147594/1‑150 (MQ=255)
                                                                                 gAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGcc                                                              <  2:1704227/150‑1 (MQ=255)
                                                                                  aGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGccc                                                             <  1:991405/150‑1 (MQ=255)
                                                                                      tGCTGGAACCATTCAGTCTTGAATGCGTGAGAAAACCGCCCGTGCGTGCAGGCAAGAACCGGCGAATTATCAGAATTTTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCcgcg                                                         <  2:188330/150‑1 (MQ=255)
                                                                                      tGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCcgcg                                                         <  2:1713914/150‑1 (MQ=255)
                                                                                      ggCTGGAACCATTCAGACTTGAATGCGTGAGACTACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTGTTCTCCGGTGTCACGGTGTTACTCAGCGGGAGTTCAGCCCcgcg                                                         <  2:1706864/149‑1 (MQ=255)
                                                                                       gCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCTCGGTGTTACTCAGCGGGAGTTCAGCCCCgcgc                                                        >  1:156964/1‑150 (MQ=255)
                                                                                       gCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCgcgc                                                        <  1:1126548/150‑1 (MQ=255)
                                                                                          ggAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTTTCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAg                                                     <  1:1023476/150‑1 (MQ=255)
                                                                                          ggAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGAGCAGGCAAGAACTGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAg                                                     >  2:387967/1‑150 (MQ=255)
                                                                                            aaCCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAACAt                                                   >  2:725863/1‑150 (MQ=255)
                                                                                             gtccggtaGACTTGACTTCGTGAGAATACCGCCCGTGCGGGCAGGCAAGAACGGGCTAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGAGCAAGAtt                                                  <  2:1636422/143‑1 (MQ=255)
                                                                                             aCCATTCAGACTGGAATGCGTGAGAATACCGCCCGGGCGGGCAGGCAAGAACGGGCGAATTATCAGGAGTGTAATATTTTGCGTAGCGTGCCATCGGTTTGCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGAtt                                                  <  1:290614/150‑1 (MQ=255)
                                                                                              ccATTCAGACTTGAATGCGTGAGAATTCCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATATTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTg                                                 <  1:1320304/150‑1 (MQ=255)
                                                                                               cATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGACTTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGt                                                <  1:1650273/150‑1 (MQ=255)
                                                                                                aTTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGCACGGGCGAATTATCAGAATTGGAATATTTTGCGTTGCGTGCCATCGGGGTTCTCCGGTGGCCCGGTGTTACTCCGCGGGAGTTCTGCCCCGCGCAAGATTGTa                                               >  2:1563339/1‑150 (MQ=255)
                                                                                                 ttCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAg                                              <  2:1355199/150‑1 (MQ=255)
                                                                                                  tCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTTTCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGGGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGa                                             <  1:989913/150‑1 (MQ=255)
                                                                                                  tCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTgaaag                                             <  2:974851/150‑5 (MQ=255)
                                                                                                   cAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCAGCGGTTTTCTCCGGTGGCACGGTGTTACACAGCGGGAGTTCTGCCCCGCGCAAGATTGTAGAt                                            >  2:1454388/1‑150 (MQ=255)
                                                                                                        ttGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGtt                                       <  2:497603/150‑1 (MQ=255)
                                                                                                               cGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATtctc                                >  1:898951/1‑150 (MQ=255)
                                                                                                               cGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATtctc                                >  1:1304161/1‑150 (MQ=255)
                                                                                                                gTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTcc                               <  2:1751793/150‑1 (MQ=255)
                                                                                                                gTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTcc                               <  2:860343/150‑1 (MQ=255)
                                                                                                                gTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTcc                               <  2:3739/150‑1 (MQ=255)
                                                                                                                  gagaATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTg                             >  1:1800269/1‑150 (MQ=255)
                                                                                                                  gagaATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTg                             <  1:438669/150‑1 (MQ=255)
                                                                                                                   gataattcGCCCGTTCTTGCCTGCACGCACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGa                            >  1:797532/8‑150 (MQ=255)
                                                                                                                      aattcGCCCGTTCTTGCCTGCACGCACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTACCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGaaaa                         >  2:1393728/5‑150 (MQ=255)
                                                                                                                      aTACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTAAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGaaaa                         >  2:618091/1‑150 (MQ=255)
                                                                                                                        aCCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAg                       <  2:1100614/150‑1 (MQ=255)
                                                                                                                          cGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCa                     <  1:529609/150‑1 (MQ=255)
                                                                                                                              cgtgcgtgCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGaaa                 >  2:1594925/1‑150 (MQ=255)
                                                                                                                               gtgcgtgcAGGCAAGCACGGGCGAAATTTAAGAATTGTAATATATTGCGTTGCGAGCCATCGTTTTTCTCTCGTGGCACGGGGTTACTCAGCGGGAGGTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCATGAAAAAGCAGaaaa                >  2:1829078/1‑150 (MQ=255)
                                                                                                                                  cgtgcAGGCAAGACCGGGCGAATTATCAGAATTGTTATATTTTGCGTTGCGTGCCATCGGTTTTCTGCGGTGGCACGGGGTTACTCAGCGGGAGTTCAGGCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGcc             <  1:242107/150‑1 (MQ=255)
                                                                                                                                  cgtgcAGGCAACCACGGGCGAATTATCAGAATTTTAATATTTTGATTTGCGGGACTTCTGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGcc             <  2:77717/150‑1 (MQ=255)
                                                                                                                                     gcAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGc          <  2:1538112/150‑1 (MQ=255)
                                                                                                                                     gcAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGc          <  1:452909/150‑1 (MQ=255)
                                                                                                                                      cAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCt         <  1:1504347/150‑1 (MQ=255)
                                                                                                                                        ggCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCttt       <  2:169874/150‑1 (MQ=255)
                                                                                                                                         gCAAGACCGGGCGAATTATCAGTATTTGAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCATCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCtttt      <  1:144768/150‑1 (MQ=255)
                                                                                                                                         gCAAGAACGGGCGAATTATCAGATTTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTCGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCtttt      <  1:139252/150‑1 (MQ=255)
                                                                                                                                            aGAACGGGCGAATTATCGGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTAtt   <  1:769948/150‑1 (MQ=255)
                                                                                                                                             gAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTc  <  2:18579/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CGAGACTGAGAGCTTTCTTAACCTCGAACCCGCGCCTGCGGTAAGAATGAATCAGCCATTCGGCCTGTTCTGCAGTGCATGGAGGGTGCTGGAACCATTCAGACTTGAATGCGTGAGAATACCGCCCGTGCGTGCAGGCAAGAACGGGCGAATTATCAGAATTGTAATATTTTGCGTTGCGTGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTC  >  NC_000913/1644398‑1644688

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: