Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,814,461 C→T L91L (CTG→CTA)  luxS ← S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,4610CT98.0% 1526.8 / 10.3 443L91L (CTG→CTA) luxSS‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base T (202/230);  total (206/237)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 443 total aligned reads displayed.

GCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGT  >  NC_000913/2814316‑2814603
                                                                                                                                                 |                                                                                                                                              
gcgcTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTc                                                                                                                                          >  7:686669/1‑151 (MQ=255)
   gttCCTTCAACGAGTGCATCTTTTAAGTGCCACACTGGTAGACGTTCAGTTCCGGCATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCata                                                                                                                                       <  8:125718/150‑1 (MQ=255)
   cTTCATGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGAATCTGATCCTGCACTATCAGAACGTCTTACATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATAAGACTCat                                                                                                                                        >  6:291407/1‑150 (MQ=255)
      ccTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCAtata                                                                                                                                     <  8:1343420/150‑1 (MQ=255)
      ccTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATaa                                                                                                                                    <  7:742283/151‑1 (MQ=255)
       cTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATaa                                                                                                                                    >  4:1056602/1‑150 (MQ=255)
       cTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATaa                                                                                                                                    >  3:1412808/1‑150 (MQ=255)
        tGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATaaa                                                                                                                                   >  2:1201573/1‑150 (MQ=255)
         gCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATaaaaa                                                                                                                                 <  4:1672175/151‑1 (MQ=255)
           aaCGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAAcc                                                                                                                               <  5:807228/151‑1 (MQ=255)
           aaCGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAAcc                                                                                                                               >  8:1533537/1‑151 (MQ=255)
            aCGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCg                                                                                                                              >  1:1578241/1‑151 (MQ=255)
             cGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGCTCCTGCACTTTCAGCACGGCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAAcccg                                                                                                                             >  2:1520063/1‑149 (MQ=255)
             cGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCgg                                                                                                                             >  6:1116686/1‑151 (MQ=255)
              gAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATAAGACTCATATAAAAACCGGt                                                                                                                            <  2:1429541/151‑1 (MQ=255)
                gTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGc                                                                                                                          <  8:1034227/151‑1 (MQ=255)
                    aTCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGcc                                                                                                                      >  5:1351374/1‑150 (MQ=255)
                        ggTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGcc                                                                                                                   <  2:1617181/150‑1 (MQ=255)
                        ggTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGcc                                                                                                                   <  5:363005/150‑1 (MQ=255)
                           aaGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCAtt                                                                                                               >  4:754821/1‑151 (MQ=255)
                           aaGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCAtt                                                                                                               <  5:506363/151‑1 (MQ=255)
                           aaGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTCCTGCAGCCCaat                                                                                                               >  2:1146450/1‑149 (MQ=255)
                                    acTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGat                                                                                                       >  6:1495357/1‑150 (MQ=255)
                                     cTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGatat                                                                                                     >  6:1696617/1‑151 (MQ=255)
                                      tGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATc                                                                                                    >  6:1445102/1‑151 (MQ=255)
                                       ggTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGatat                                                                                                     <  5:967985/149‑1 (MQ=255)
                                         tAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGAt                                                                                                 >  7:516420/1‑151 (MQ=255)
                                           gACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTTCGGCAGCCCATTGGCGAGATAACGATa                                                                                                >  6:465737/1‑150 (MQ=255)
                                            aCGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAAt                                                                                              <  8:203044/151‑1 (MQ=255)
                                             cGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGagactc                                                                                             >  6:880322/1‑146 (MQ=255)
                                             agTTCAGTTCCGTTTTCTGATTCTGATCCTGGTCTTTCAGGACGTGTTCCATTGCCGCTTTCCAGGCATCAGCAAGACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCGCATTGGCGAGATATCGATAAt                                                                                              <  4:900318/149‑1 (MQ=255)
                                              gTTCAGTTCCGGGATCTGATTCTGATACTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCAGTGGCGAGATATCGATAAtct                                                                                            >  5:1493706/1‑151 (MQ=255)
                                               ttCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAAtct                                                                                            >  4:539927/1‑150 (MQ=255)
                                                 cAGTTCCGTTATGTGATTCTGATCCTTCACTTTCAGGACGTCTTCCATTGCCGATTTCCAGGCATCAGCAACACGATGCTCATCTGGCGTACC‑ATTATACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATctct                                                                                          <  1:138787/150‑1 (MQ=255)
                                                 cAGTTCCGTGATCTGATTCTGATGCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTa                                                                                         <  8:600326/151‑1 (MQ=255)
                                                     tCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTacac                                                                                      <  5:1463367/150‑1 (MQ=255)
                                                        gggATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCAtt                                                                                  >  1:490140/1‑151 (MQ=255)
                                                          gATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCAACTGGCGTACC‑ATTAGACTCATATAAAAACCGCTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCAATAc                                                                                >  7:1552559/1‑151 (MQ=255)
                                                            tctgattctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAAACGGTGCGGCAGCACATTGGCGAGATATCGATAATCTCTACACCATTACCg                                                                              >  1:609704/1‑151 (MQ=255)
                                                             ctgattctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGt                                                                             >  1:1496855/1‑151 (MQ=255)
                                                              tgattctgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGtt                                                                            >  3:1636048/1‑151 (MQ=255)
                                                                attctggtCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTaa                                                                          <  2:514812/151‑1 (MQ=255)
                                                                attctgatCCTGCACTTTCAGCACGGCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTaa                                                                          <  3:473546/151‑1 (MQ=255)
                                                                    tgatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATg                                                                      >  6:742673/1‑151 (MQ=255)
                                                                     gatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGAt                                                                       <  4:1449191/149‑1 (MQ=255)
                                                                           gCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACg                                                                <  3:771145/150‑1 (MQ=255)
                                                                            cACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGc                                                               >  3:1083118/1‑150 (MQ=255)
                                                                            cACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGc                                                               >  3:1195949/1‑150 (MQ=255)
                                                                             aCTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACGCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCAt                                                             >  6:559930/1‑151 (MQ=255)
                                                                               tttCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATa                                                            >  7:1074238/1‑150 (MQ=255)
                                                                                ttCAGGACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaa                                                           <  7:35132/150‑1 (MQ=255)
                                                                                ttCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaa                                                           >  5:1025340/1‑150 (MQ=255)
                                                                                ttCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaa                                                           >  4:6112/1‑150 (MQ=255)
                                                                                ttCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaa                                                           >  1:1174101/1‑150 (MQ=255)
                                                                                 tCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaaa                                                          <  5:1138213/150‑1 (MQ=255)
                                                                                 tCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaaa                                                          <  3:1140864/150‑1 (MQ=255)
                                                                                  cAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATaaaa                                                         <  7:290307/150‑1 (MQ=255)
                                                                                   aGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAAcc                                                       <  4:1172093/151‑1 (MQ=255)
                                                                                   aGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAAcc                                                       >  5:171850/1‑151 (MQ=255)
                                                                                    gCACGTCTTCGATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCa                                                      <  2:786590/151‑1 (MQ=255)
                                                                                     cACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAg                                                     <  7:869326/151‑1 (MQ=255)
                                                                                      aCGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGc                                                    >  8:1457806/1‑151 (MQ=255)
                                                                                         tCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCaa                                                  >  7:610354/1‑150 (MQ=255)
                                                                                           ttCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCaaaat                                               >  1:1006467/1‑149 (MQ=255)
                                                                                            tCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACa                                               <  8:678303/150‑1 (MQ=255)
                                                                                            tCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACa                                               <  7:1471394/150‑1 (MQ=255)
                                                                                             cgaTTGCGGCTTTGCAGGCATCAGGAAGACGGTGCTCATCTGGCGTA‑CAATTAGACTCATATAAAAAGCGGTGCGGCAGCCCATTGGCGAGATATCGATAATGTCTACACCATTACGGTTAAGATGTTTACGCATAATACAAGCTAACAgg                                             <  6:390121/149‑1 (MQ=255)
                                                                                             ccATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAg                                              <  8:898940/150‑1 (MQ=255)
                                                                                                 tcccGCTTTCCAGGCATCAGCAACACGCTGCTCATATGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCTCATTGGCGAGATATCGATAATCTCTACGCAATTACCGATAAGATCGTTACGCATAAAACCAGCAAACAGCTGcg                                         >  3:245123/3‑150 (MQ=255)
                                                                                                 tGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGc                                          <  1:325518/150‑1 (MQ=255)
                                                                                                  gCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCt                                         <  4:803963/150‑1 (MQ=255)
                                                                                                    cGCTTTCCAGGCATCAGCAACACGCTTCTCATCTGGCGTACC‑ATTACACACAAATAAAAACCCGGGCGGCACCCCATTGGCGAGAAATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCa                                      <  2:636924/151‑1 (MQ=255)
                                                                                                     gCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGGGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAg                                     <  6:1407914/151‑1 (MQ=255)
                                                                                                        ttCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGt                                  >  5:1213402/1‑151 (MQ=255)
                                                                                                        ttCCAGGCAACAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAggg                                   <  3:1571972/150‑1 (MQ=255)
                                                                                                          ccAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTa                                 <  5:1398012/150‑1 (MQ=255)
                                                                                                          ccAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTa                                 <  3:900276/150‑1 (MQ=255)
                                                                                                          ccAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTAt                                <  8:444602/151‑1 (MQ=255)
                                                                                                          ccAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCACATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTAt                                <  1:233627/151‑1 (MQ=255)
                                                                                                           cAGGCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATg                               >  7:1225745/1‑151 (MQ=255)
                                                                                                              gCATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGAt                            >  3:1335887/1‑151 (MQ=255)
                                                                                                               cATCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATc                           >  7:1379954/1‑151 (MQ=255)
                                                                                                                 tCAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATccc                         >  2:92267/1‑151 (MQ=255)
                                                                                                                  cAGCAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCct                        >  4:411508/1‑151 (MQ=255)
                                                                                                                     cAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCtt                     >  3:378331/1‑151 (MQ=255)
                                                                                                                     cAACACGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCtt                     <  2:136545/151‑1 (MQ=255)
                                                                                                                       acacGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTc                   <  6:681278/151‑1 (MQ=255)
                                                                                                                         acGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTg                 <  8:1313745/151‑1 (MQ=255)
                                                                                                                          cGCTGCTCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTgg                <  1:983742/151‑1 (MQ=255)
                                                                                                                              gctCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCAt             <  3:658435/150‑1 (MQ=255)
                                                                                                                               ctCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATc            >  7:174026/1‑150 (MQ=255)
                                                                                                                                tCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAc          >  2:327566/1‑151 (MQ=255)
                                                                                                                                tCATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAc          >  4:951908/1‑151 (MQ=255)
                                                                                                                                 cATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTAAACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAct         >  6:1270593/1‑151 (MQ=255)
                                                                                                                                 cATCTGGCGTACC‑ATTAGACTCATATAAAAACCGGGGCGGCAGCCCATTCGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATCGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAct         >  2:259755/1‑151 (MQ=255)
                                                                                                                                  aTCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCActt        >  8:812173/1‑151 (MQ=255)
                                                                                                                                  aTCTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCActt        >  5:738973/1‑151 (MQ=255)
                                                                                                                                    cTGGCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCACTCTTTCTGGCATCActtct      >  4:276235/1‑151 (MQ=255)
                                                                                                                                      ggCGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCttt    >  1:1692937/1‑151 (MQ=255)
                                                                                                                                        cGTACC‑ATTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGt  <  8:1693391/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGT  >  NC_000913/2814316‑2814603

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