Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,174 T→A 9.5% 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,1740TA 9.5% 303.0 / 18.5 105intergenic (+20/+44)emrB/luxSmultidrug efflux system protein/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (42/53);  new base A (4/6);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 110 total aligned reads displayed.

AGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTC  >  NC_000913/2814027‑2814322
                                                                                                                                                   |                                                                                                                                                    
aGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                                                                     >  2:254553/1‑149 (MQ=255)
 ggggCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGTATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                                    <  2:1025145/149‑1 (MQ=255)
      tGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAGACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACt                                                                                                                                               <  1:362481/149‑1 (MQ=255)
       gATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTg                                                                                                                                              >  1:137047/1‑149 (MQ=255)
        attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTgg                                                                                                                                             >  2:593869/1‑149 (MQ=255)
        attattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                             <  1:394334/149‑1 (MQ=255)
         ttattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGTTGCTGGTGTGGTTTGCTAAACCGCCAGTTGGCGCAGGGGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGc                                                                                                                                            <  1:1428028/149‑1 (MQ=255)
            tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttt                                                                                                                                         >  1:872765/1‑149 (MQ=255)
            tttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGACCTGGCtta                                                                                                                                         <  1:1355844/149‑2 (MQ=255)
             ttCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCttag                                                                                                                                        <  1:463049/149‑3 (MQ=255)
              tCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACAAAGCCAGTTCAAATGAACTGGCttagt                                                                                                                                       <  1:913924/149‑4 (MQ=255)
               ccGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttttt                                                                                                                                      <  2:397689/149‑1 (MQ=255)
                 gCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCa                                                                                                                                    <  1:1469372/149‑1 (MQ=255)
                  ccACTGCGATCTTCTGGATGTCCGCCGTGATATTCCGCGCCCTGCCGGGGCGTGTGTGGTTTGCTACACCGCCATTTGGCGCAGGTGGCGGCGTAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCttagttgta                                                                                                                                   <  1:112949/149‑8 (MQ=255)
                  ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTGaa                                                                                                                                   >  1:302124/1‑149 (MQ=255)
                   cAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaat                                                                                                                                  <  2:610436/149‑1 (MQ=255)
                        agagCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaatt                                                                                                                             <  1:17950/149‑1 (MQ=255)
                        agaTCTTCTTGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGACAGTTCATTTGAACTGGCTTTTTTCaattaatt                                                                                                                             <  1:792956/149‑1 (MQ=255)
                        agaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaatt                                                                                                                             <  1:914436/149‑1 (MQ=255)
                        agaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaatt                                                                                                                             <  1:1037287/149‑1 (MQ=255)
                         gaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtg                                                                                                                            >  1:311556/1‑149 (MQ=255)
                          aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                                                                           >  1:554731/1‑149 (MQ=255)
                          aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                                                                           >  2:701252/1‑149 (MQ=255)
                          aTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                                                                           >  2:1517752/1‑149 (MQ=255)
                           tcttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACCGGCTTTTTTCAATTAATtgtg                                                                                                                          >  2:88105/1‑149 (MQ=255)
                             ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                                                                        <  2:1366382/149‑1 (MQ=255)
                             ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                                                                        >  2:817173/1‑149 (MQ=255)
                             ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                                                                        >  1:1173703/1‑149 (MQ=255)
                             ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                                                                        <  2:942874/149‑1 (MQ=255)
                                tGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGAt                                                                                                                     >  2:423382/1‑149 (MQ=255)
                                 ggATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATa                                                                                                                    <  2:652684/149‑1 (MQ=255)
                                 ggATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATa                                                                                                                    <  1:1231147/149‑1 (MQ=255)
                                 ggATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATa                                                                                                                    >  2:188883/1‑149 (MQ=255)
                                  gATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAg                                                                                                                   >  1:1064859/1‑149 (MQ=255)
                                  gATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAg                                                                                                                   >  1:555554/1‑149 (MQ=255)
                                      tCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                                                                               <  1:1582605/149‑1 (MQ=255)
                                        aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                                                             <  1:1371859/149‑1 (MQ=255)
                                        aGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                                                             <  2:684048/149‑1 (MQ=255)
                                                tatTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAt                                                                                                     <  1:1484235/149‑1 (MQ=255)
                                                  ttCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgt                                                                                                   >  2:274796/1‑149 (MQ=255)
                                                  ttCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgt                                                                                                   <  2:1465461/149‑1 (MQ=255)
                                                  ttCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGCGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgt                                                                                                   >  2:1215968/1‑149 (MQ=255)
                                                   tCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgtg                                                                                                  <  1:254553/149‑1 (MQ=255)
                                                     cTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                                >  1:1005851/1‑149 (MQ=255)
                                                     cTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                                <  1:1248073/149‑1 (MQ=255)
                                                            tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTg                                                                                         <  2:757324/149‑1 (MQ=255)
                                                              ctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCa                                                                                       >  2:729872/1‑149 (MQ=255)
                                                              ctgCGGCTGGTGTGGTTTGCTATACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCa                                                                                       <  1:162022/149‑1 (MQ=255)
                                                                ggggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                                                                    <  2:114848/134‑1 (MQ=255)
                                                                ggggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtg                                                                                                    >  1:114848/1‑134 (MQ=255)
                                                                  ggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACtt                                                                                   <  1:917410/149‑1 (MQ=255)
                                                                  ggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACtt                                                                                   <  1:1359576/149‑1 (MQ=255)
                                                                  ggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACtt                                                                                   <  2:833034/149‑1 (MQ=255)
                                                                   gCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACtta                                                                                  >  2:1480649/1‑148 (MQ=255)
                                                                     tGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtctc                                                                                >  2:1566218/1‑149 (MQ=255)
                                                                      ggTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGTTAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTctct                                                                               <  2:782242/149‑1 (MQ=255)
                                                                      ggTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGCGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTctct                                                                               <  1:1215968/149‑1 (MQ=255)
                                                                         gtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAGTTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTc                                                                            <  2:980178/149‑1 (MQ=255)
                                                                                gCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTg                                                                     >  1:1296750/1‑149 (MQ=255)
                                                                                 cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                                    >  1:388259/1‑149 (MQ=255)
                                                                                  tAAACCGCCATTTTGCGCCGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTTGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                                                                   <  1:793328/149‑1 (MQ=255)
                                                                                  tAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                                                                   >  1:1024800/1‑149 (MQ=255)
                                                                                   aaaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                                                                  >  2:1035079/1‑149 (MQ=255)
                                                                                       cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc                                                              >  1:1124270/1‑149 (MQ=255)
                                                                                          cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                                                           <  2:311556/149‑1 (MQ=255)
                                                                                          cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                                                           <  1:1491077/149‑1 (MQ=255)
                                                                                            tttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGt                                                         <  1:593869/149‑1 (MQ=255)
                                                                                                 cgcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg                                                    >  2:1054506/1‑149 (MQ=255)
                                                                                                     ggTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGa                                                <  2:689367/149‑1 (MQ=255)
                                                                                                      gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                               <  1:274796/149‑1 (MQ=255)
                                                                                                      gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                               <  2:302124/149‑1 (MQ=255)
                                                                                                      gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                               >  1:1303971/1‑149 (MQ=255)
                                                                                                      gTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGGGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGAt                                               <  2:1134484/149‑1 (MQ=255)
                                                                                                       tGGCGGCGGAGGCGGTGCACACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATg                                              <  1:499634/149‑1 (MQ=255)
                                                                                                          cggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGt                                           <  2:872765/149‑1 (MQ=255)
                                                                                                            gcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTAc                                         <  1:423382/149‑1 (MQ=255)
                                                                                                                aGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCa                                     <  1:429533/149‑1 (MQ=255)
                                                                                                                 ggCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCAc                                    <  2:388259/149‑1 (MQ=255)
                                                                                                                     gTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTc                                <  1:1018395/149‑1 (MQ=255)
                                                                                                                      tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTcc                               <  1:817173/149‑1 (MQ=255)
                                                                                                                       gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCa                              <  2:98713/149‑1 (MQ=255)
                                                                                                                       gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCa                              <  2:784602/149‑1 (MQ=255)
                                                                                                                       gcgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCa                              <  1:582529/149‑1 (MQ=255)
                                                                                                                          cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGaa                           <  2:1056975/149‑1 (MQ=255)
                                                                                                                          cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTATTCGTTGCTGTTGATGCGTACGTCACGTTCCAGaa                           >  1:679715/1‑149 (MQ=255)
                                                                                                                             tAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGc                        >  2:507481/1‑149 (MQ=255)
                                                                                                                                 ttgaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcg                    >  1:466849/8‑149 (MQ=255)
                                                                                                                                   gaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgcg                  >  2:664295/6‑149 (MQ=255)
                                                                                                                                    aaCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTATTTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACgcgc                 >  1:605131/1‑149 (MQ=255)
                                                                                                                                    aaCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCCACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATTCGTACGTCACGTTCCAGAATGCTACgcgc                 >  2:672634/1‑149 (MQ=255)
                                                                                                                                    aaCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATAGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACgcgc                 >  1:1506950/1‑149 (MQ=255)
                                                                                                                                       aaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAat              >  1:1010334/2‑149 (MQ=255)
                                                                                                                                        aaGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAata             >  2:824760/1‑149 (MQ=255)
                                                                                                                                         aGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTTCGTCACGTTCCAGAATGCTACGCGCAatat            >  1:246835/1‑149 (MQ=255)
                                                                                                                                         aGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAatat            >  1:1600226/1‑149 (MQ=255)
                                                                                                                                          gCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATc           <  2:1064859/149‑1 (MQ=255)
                                                                                                                                           ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATcc          >  1:1127460/1‑149 (MQ=255)
                                                                                                                                               ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTgcg      >  2:1153687/1‑149 (MQ=255)
                                                                                                                                                   tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTc  <  1:615441/149‑1 (MQ=255)
                                                                                                                                                   tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTc  <  1:1566218/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGGGCTGATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTC  >  NC_000913/2814027‑2814322

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