Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,174 T→A 11.0% intergenic (+20/+44) emrB → / ← luxS multidrug efflux system protein/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,1740TA11.0% 312.4 / 25.2 100intergenic (+20/+44)emrB/luxSmultidrug efflux system protein/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (35/54);  new base A (5/6);  total (40/60)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 105 total aligned reads displayed.

GGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTC  >  NC_000913/2814028‑2814322
                                                                                                                                                  |                                                                                                                                                    
ggggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                                    <  4:925647/149‑1 (MQ=255)
ggggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                                    <  3:770330/149‑1 (MQ=255)
 gggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTg                                                                                                                                                   <  4:676015/149‑1 (MQ=255)
   gCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGaa                                                                                                                                                 >  3:304174/1‑149 (MQ=255)
    cTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGaa                                                                                                                                                 >  2:58269/1‑148 (MQ=255)
       attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                             <  4:888311/149‑1 (MQ=255)
       attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAAGTGGCGGCGGAGGCGGTGCGCACCAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                             <  4:743259/149‑1 (MQ=255)
        ttattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGc                                                                                                                                            <  3:218949/149‑1 (MQ=255)
        ttattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGc                                                                                                                                            <  3:291342/149‑1 (MQ=255)
          attTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCt                                                                                                                                           <  1:836375/148‑1 (MQ=255)
           tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttt                                                                                                                                         >  3:47736/1‑149 (MQ=255)
           tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttt                                                                                                                                         >  4:524078/1‑149 (MQ=255)
           tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtta                                                                                                                                         <  4:359671/149‑2 (MQ=255)
            ttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTAGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttt                                                                                                                                        <  4:327893/149‑1 (MQ=255)
             tCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttttt                                                                                                                                       <  3:945832/149‑1 (MQ=255)
               cGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTc                                                                                                                                     <  3:877023/149‑1 (MQ=255)
                 ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaa                                                                                                                                   <  3:132468/149‑1 (MQ=255)
                 ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCa                                                                                                                                    <  1:58269/148‑1 (MQ=255)
                 ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCa                                                                                                                                    <  1:664543/148‑1 (MQ=255)
                  cAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaa                                                                                                                                   <  1:394068/148‑1 (MQ=255)
                   aaTGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaatt                                                                                                                                 >  4:238063/1‑149 (MQ=255)
                      gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                               >  2:628750/1‑148 (MQ=255)
                      gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                               >  2:901952/1‑148 (MQ=255)
                       agaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaat                                                                                                                              <  2:316569/148‑1 (MQ=255)
                        gaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtg                                                                                                                            >  4:384719/1‑149 (MQ=255)
                         aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                                                                           >  3:549569/1‑149 (MQ=255)
                           cttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                                                                         <  4:467568/149‑1 (MQ=255)
                                   tGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCATTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                                                                 <  3:887769/149‑1 (MQ=255)
                                   tGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                                                                 <  3:228648/149‑1 (MQ=255)
                                       aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTAc                                                                                                              <  1:382759/148‑1 (MQ=255)
                                       aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTAc                                                                                                              <  2:451646/148‑1 (MQ=255)
                                       aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                                                             <  3:230337/149‑1 (MQ=255)
                                        gCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                                                                         <  3:950400/136‑1 (MQ=255)
                                         ccGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAGCTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGa                                                                                                           <  4:700925/149‑1 (MQ=255)
                                         ccGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                                                                            <  2:323032/148‑1 (MQ=255)
                                          cGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGa                                                                                                           <  1:628750/148‑1 (MQ=255)
                                           gggATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                                                          >  1:594295/1‑148 (MQ=255)
                                             gATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAg                                                                                                       <  4:198834/149‑1 (MQ=255)
                                             gATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAg                                                                                                       >  4:874743/1‑149 (MQ=255)
                                              atatTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                                                                      <  4:787829/149‑1 (MQ=255)
                                               tatTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                                                                      <  2:466762/148‑1 (MQ=255)
                                                    cTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                                >  3:517835/1‑149 (MQ=255)
                                                     tCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTgg                                                                                                                                             >  2:438789/1‑103 (MQ=255)
                                                     tCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTgg                                                                                                                                             <  1:438789/103‑1 (MQ=255)
                                                         ccTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTccc                                                                                           <  3:574568/149‑1 (MQ=255)
                                                              tgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCaa                                                                                      <  4:13512/149‑1 (MQ=255)
                                                                 ggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACt                                                                                    >  2:534216/1‑148 (MQ=255)
                                                                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                                                                    >  1:513628/1‑130 (MQ=255)
                                                                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                                                                    <  2:513628/130‑1 (MQ=255)
                                                                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtc                                                                                  >  1:562801/1‑148 (MQ=255)
                                                                       tgtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                               >  2:916287/1‑99 (MQ=255)
                                                                       tgtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaa                                                                                                                               <  1:916287/99‑1 (MQ=255)
                                                                          ggTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCgg                                                                          <  3:384719/149‑1 (MQ=255)
                                                                            tttGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                         >  1:438361/1‑148 (MQ=255)
                                                                               gCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTTATTGTGAAGATAGtt                                                                                                                 <  3:437283/105‑1 (MQ=255)
                                                                               gCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                                                                 >  4:437283/1‑105 (MQ=255)
                                                                                 tAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                             <  2:668497/139‑1 (MQ=255)
                                                                                 tAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                             >  1:668497/1‑139 (MQ=255)
                                                                                  aaaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                                                                   >  2:502465/1‑148 (MQ=255)
                                                                                   aaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                                                                  >  1:868678/1‑148 (MQ=255)
                                                                                    aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACAGACTAGTTGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGt                                                                >  4:802981/1‑149 (MQ=255)
                                                                                      cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTGCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc                                                              >  3:138243/1‑149 (MQ=255)
                                                                                      cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGtt                                                               >  2:881953/1‑148 (MQ=255)
                                                                                       gCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttct                                                             <  3:861187/149‑1 (MQ=255)
                                                                                         cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTTAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttctt                                                            <  1:259195/148‑1 (MQ=255)
                                                                                         cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                                                           <  3:475991/149‑1 (MQ=255)
                                                                                         cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttctt                                                            <  2:307250/148‑1 (MQ=255)
                                                                                           tttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCg                                                          <  2:594295/148‑1 (MQ=255)
                                                                                            ttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGtt                                                        <  4:619056/149‑1 (MQ=255)
                                                                                            ttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGtt                                                        <  4:378041/149‑1 (MQ=255)
                                                                                                  cAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATACTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGtt                                                  <  4:549569/149‑1 (MQ=255)
                                                                                                   aGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTg                                                 <  3:7013/149‑1 (MQ=255)
                                                                                                   aGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTg                                                 >  3:41574/1‑149 (MQ=255)
                                                                                                          ggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTa                                          >  3:864183/1‑149 (MQ=255)
                                                                                                           gcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTa                                          >  2:665265/1‑148 (MQ=255)
                                                                                                              gAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGt                                       <  1:534216/148‑1 (MQ=255)
                                                                                                                     tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTcc                               <  3:469944/149‑1 (MQ=255)
                                                                                                                      gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTcc                               <  2:274113/148‑1 (MQ=255)
                                                                                                                              aGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt                       <  1:430953/148‑1 (MQ=255)
                                                                                                                              aGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                      <  3:828041/149‑1 (MQ=255)
                                                                                                                              aGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                      <  4:654137/149‑1 (MQ=255)
                                                                                                                               gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTa                      <  1:502465/148‑1 (MQ=255)
                                                                                                                               gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                     <  4:694208/149‑1 (MQ=255)
                                                                                                                                  gaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgcg                  >  3:244594/6‑149 (MQ=255)
                                                                                                                                  cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgc                   >  1:723432/1‑148 (MQ=255)
                                                                                                                                  cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgc                   >  2:280120/1‑148 (MQ=255)
                                                                                                                                     cTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCCGAATGCTACGCGCa                >  2:289613/1‑148 (MQ=255)
                                                                                                                                     cTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCa                <  1:381546/148‑1 (MQ=255)
                                                                                                                                     aaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCa                >  2:307020/3‑148 (MQ=255)
                                                                                                                                        aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCTAg                                                                 <  1:556310/96‑1 (MQ=255)
                                                                                                                                        aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAg                                                                 >  2:556310/1‑96 (MQ=255)
                                                                                                                                        aGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAata             >  2:764203/1‑148 (MQ=255)
                                                                                                                                         gCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                     <  4:907359/139‑1 (MQ=255)
                                                                                                                                         gCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAc                     >  3:907359/1‑139 (MQ=255)
                                                                                                                                          ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATc           >  1:409060/1‑148 (MQ=255)
                                                                                                                                           cAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCACTATcc          >  1:252113/1‑148 (MQ=255)
                                                                                                                                             gTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgc       >  3:654604/1‑149 (MQ=255)
                                                                                                                                                  tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                             >  1:869282/1‑74 (MQ=255)
                                                                                                                                                  tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                                                                             <  2:869282/74‑1 (MQ=255)
                                                                                                                                                  tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTc  >  3:726463/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTC  >  NC_000913/2814028‑2814322

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