Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 921,700 G→A 5.7% V452I (GTC→ATC)  macB → ABC‑type tripartite efflux pump ATP binding/membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913921,7000GA 5.7% 223.5 / 2.2 70V452I (GTC→ATC) macBABC‑type tripartite efflux pump ATP binding/membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base A (2/2);  total (33/37)
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.41e-01

TCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAA  >  NC_000913/921555‑921845
                                                                                                                                                 |                                                                                                                                                 
tctcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTg                                                                                                                                               <  2:1008167/150‑1 (MQ=255)
    gACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTa                                                                                                                                           >  1:1303259/1‑150 (MQ=255)
        gTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCtt                                                                                                                                       >  1:902577/1‑150 (MQ=255)
               tCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCAc                                                                                                                                <  1:505518/150‑1 (MQ=255)
                cgccgcCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGTCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACg                                                                                                                               >  2:717410/1‑150 (MQ=255)
                cgccgcCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACg                                                                                                                               >  1:1291910/1‑150 (MQ=255)
                    gccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGt                                                                                                                           >  2:367463/1‑150 (MQ=255)
                    gccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTATCTGGCTACCTTACAGCACGATGt                                                                                                                           >  2:650244/1‑150 (MQ=255)
                    gccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTATCTGGCTACCTTACAGCACGATGt                                                                                                                           >  2:645215/1‑150 (MQ=255)
                     ccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTc                                                                                                                          >  1:95494/1‑150 (MQ=255)
                     ccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTc                                                                                                                          <  2:468261/150‑1 (MQ=255)
                      cAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTcc                                                                                                                         >  1:729288/1‑150 (MQ=255)
                           tttCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCg                                                                                                                    <  2:778843/150‑1 (MQ=255)
                             tCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAg                                                                                                                  <  2:760997/150‑1 (MQ=255)
                                  cATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCGAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCGGCGTGGCTGGCTACCTTACAGCACGATGTGCGGGCGAGTTATg                                                                                                             <  2:900923/150‑1 (MQ=255)
                                    tAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCTATTTTTGGTTGGAGTAAAGTGCTGCGTGTCTGGCTTCCTTACAGCACGATGTCTGGTCGAGTTATggg                                                                                                           <  2:774472/150‑1 (MQ=255)
                                    tAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATggg                                                                                                           >  1:141610/1‑150 (MQ=255)
                                    tAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATggg                                                                                                           >  2:288238/1‑150 (MQ=255)
                                      aaGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCgg                                                                                                                       >  1:526320/1‑136 (MQ=255)
                                      aaGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCgg                                                                                                                       <  2:526320/136‑1 (MQ=255)
                                      aaGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGcc                                                                                                         <  2:1227658/150‑1 (MQ=255)
                                       aGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCa                                                                                                        >  2:727050/1‑150 (MQ=255)
                                         cAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGt                                                                                                      >  1:214866/1‑150 (MQ=255)
                                          aGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGTTGTCCGGGCGAGTTATGGGCCAGTc                                                                                                     <  2:265799/150‑1 (MQ=255)
                                             tgtgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTg                                                                                                  <  2:1010107/150‑1 (MQ=255)
                                              gtgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTgg                                                                                                 <  2:1135329/150‑1 (MQ=255)
                                                        aGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCa                                                                                       <  2:95494/150‑1 (MQ=255)
                                                         ggTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCAt                                                                                      >  1:1016843/1‑150 (MQ=255)
                                                          gTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCGAGTCGTAGCTTAACTCCAtt                                                                                     <  2:698540/150‑1 (MQ=255)
                                                            gATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTAc                                                                                   <  2:230948/150‑1 (MQ=255)
                                                              ttCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTg                                                                                 >  1:842796/1‑150 (MQ=255)
                                                                 tGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTATCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCa                                                                              <  1:650244/150‑1 (MQ=255)
                                                                 tGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTATCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCa                                                                              <  1:645215/150‑1 (MQ=255)
                                                                  ggTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAg                                                                             >  2:526326/1‑150 (MQ=255)
                                                                  ggTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAg                                                                             >  1:363554/1‑150 (MQ=255)
                                                                     cGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGt                                                                          >  1:1153256/1‑150 (MQ=255)
                                                                     cGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGt                                                                          <  1:727050/150‑1 (MQ=255)
                                                                      ggCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTg                                                                         <  1:581259/150‑1 (MQ=255)
                                                                      ggCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTg                                                                         <  1:568949/150‑1 (MQ=255)
                                                                      ggCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTg                                                                         <  1:74275/150‑1 (MQ=255)
                                                                             tGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAaagaag                                                                  <  2:720063/150‑1 (MQ=255)
                                                                             tGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAaagaag                                                                  >  2:440078/1‑150 (MQ=255)
                                                                             tGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAaagaag                                                                  >  2:435193/1‑150 (MQ=255)
                                                                             tGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGGGAaagaag                                                                  <  2:902577/150‑1 (MQ=255)
                                                                                ccGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGAt                                                               >  2:792615/1‑150 (MQ=255)
                                                                                 cGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCAGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGAtt                                                              >  2:860612/1‑150 (MQ=255)
                                                                                   ccAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTg                                                            <  2:141610/150‑1 (MQ=255)
                                                                                            ttGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCg                                                   >  1:595897/1‑150 (MQ=255)
                                                                                               gtgtGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAgg                                                >  2:66969/1‑150 (MQ=255)
                                                                                               gtgtGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCGTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAgg                                                <  2:1355918/150‑1 (MQ=255)
                                                                                                tgtgGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGc                                               >  2:95210/1‑150 (MQ=255)
                                                                                                 gtgGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCg                                              <  1:526326/150‑1 (MQ=255)
                                                                                                   ggCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGa                                            >  1:203434/1‑150 (MQ=255)
                                                                                                     cGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAAc                                          <  2:363554/150‑1 (MQ=255)
                                                                                                        aagaaAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGc                                       <  1:367463/150‑1 (MQ=255)
                                                                                                          taaaaaCAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCaa                                     <  2:1124214/149‑1 (MQ=255)
                                                                                                            aaaaCAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACt                                   >  2:537794/1‑150 (MQ=255)
                                                                                                              aaCAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCa                                 <  2:595897/150‑1 (MQ=255)
                                                                                                                     gATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGttt                          <  1:13852/150‑1 (MQ=255)
                                                                                                                       tGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTAc                        >  1:705631/1‑150 (MQ=255)
                                                                                                                         tttGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACtt                      >  2:838984/1‑150 (MQ=255)
                                                                                                                            ggTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCa                   <  1:792615/150‑1 (MQ=255)
                                                                                                                             gTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCg                                                                                                                        <  2:819792/48‑1 (MQ=255)
                                                                                                                             gTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCg                                                                                                                        >  1:819792/1‑48 (MQ=255)
                                                                                                                             gTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCAc                  <  1:435193/150‑1 (MQ=255)
                                                                                                                             gTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCAc                  <  1:440078/150‑1 (MQ=255)
                                                                                                                             gTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCAc                  <  1:288238/150‑1 (MQ=255)
                                                                                                                              tAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACt                 >  1:896335/1‑150 (MQ=255)
                                                                                                                                           ctgcGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAaag    <  1:95210/150‑1 (MQ=255)
                                                                                                                                             gcGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAaagaa  >  1:720279/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTGTCAGGGTGAAAGAAGGATTTGACAGCGCCGAGGCGGAACAGCAACTCACGCGTTTACTTTCACTGCGCCACGGAAAGAA  >  NC_000913/921555‑921845

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