Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 921,575 G→A 5.9% R410H (CGC→CAC)  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,5750GA 5.9% 210.1 / 2.2 67R410H (CGC→CAC) macBABC‑type tripartite efflux pump ATP binding/membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (30/33);  new base A (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGT  >  NC_000913/921426‑921724
                                                                                                                                                     |                                                                                                                                                     
taACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTcgccg                                                                                                                                                       >  1:259805/1‑150 (MQ=255)
    aaTGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAg                                                                                                                                                   >  1:900923/1‑150 (MQ=255)
    aaTGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAg                                                                                                                                                   <  2:259805/150‑1 (MQ=255)
    aaTGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAg                                                                                                                                                   <  1:120306/150‑1 (MQ=255)
      tgttgatgttgCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCt                                                                                                                                                 >  1:1008167/1‑150 (MQ=255)
        ttgatgttgCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCttt                                                                                                                                               <  2:750374/150‑1 (MQ=255)
           atgttgCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTcc                                                                                                                                            >  1:1355918/1‑150 (MQ=255)
             gttgCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTcccc                                                                                                                                          <  2:576525/150‑1 (MQ=255)
              ttgCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTccccc                                                                                                                                         <  2:27355/150‑1 (MQ=255)
               tgctgcCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCa                                                                                                                                        <  2:295681/150‑1 (MQ=255)
                   gcCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCtttt                                                                                                                                              >  2:1427268/1‑140 (MQ=255)
                   gcCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCtttt                                                                                                                                              <  1:1427268/140‑1 (MQ=255)
                   gcCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATaaa                                                                                                                                    >  2:320978/1‑150 (MQ=255)
                   gcCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATaaa                                                                                                                                    >  1:1399412/1‑150 (MQ=255)
                      aGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCa                                                                                                                                 <  2:46330/150‑1 (MQ=255)
                          ccAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGAtg                                                                                                                             >  2:568949/1‑150 (MQ=255)
                          ccAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCACCAGCTTTTCCCCCATAAAGCAGAtg                                                                                                                             >  1:1232213/1‑150 (MQ=255)
                          ccAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCACCAGCTTTTCCCCCATAAAGCAGAtg                                                                                                                             >  1:1227171/1‑150 (MQ=255)
                           cAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGAtgt                                                                                                                            >  2:480694/1‑150 (MQ=255)
                                        ggcgATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTg                                                                                                               >  1:830522/1‑150 (MQ=255)
                                            attattTTAATGTCTATGGCATGACCTTCAGTGAAGGATACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTc                                                                                                           <  1:718777/150‑1 (MQ=255)
                                             ttattTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCCGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCt                                                                                                          <  2:1180543/150‑1 (MQ=255)
                                              tattTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTg                                                                                                         >  1:243323/1‑150 (MQ=255)
                                               attTTAATGTCTATGGCATGACCCTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTgg                                                                                                        >  1:265799/1‑150 (MQ=255)
                                                ttttAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGt                                                                                                       >  1:1010107/1‑150 (MQ=255)
                                                     aTGTCTATGGCATGACCTTCAGCGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGTTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCa                                                                                                  <  2:618885/150‑1 (MQ=255)
                                                        tCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAata                                                                                               <  2:830522/150‑1 (MQ=255)
                                                         cTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAatat                                                                                              >  2:1075148/1‑150 (MQ=255)
                                                          tATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATg                                                                                             <  2:965133/150‑1 (MQ=255)
                                                             ggCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGccc                                                                                          >  2:265838/1‑150 (MQ=255)
                                                                         aGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGtt                                                                                                                        >  2:963546/1‑108 (MQ=255)
                                                                         aGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGtt                                                                                                                        <  1:963546/108‑1 (MQ=255)
                                                                           tGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGACTCTGGCGGGCATATTGCCGACCAGAATCacctc                                                                            <  1:981565/150‑5 (MQ=38)
                                                                            gAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGgt                                                                           >  2:786336/1‑150 (MQ=255)
                                                                             aaGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGgtg                                                                          >  2:74275/1‑150 (MQ=255)
                                                                              aGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTGGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGgtgt                                                                         >  2:575170/1‑150 (MQ=255)
                                                                               ggAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGtgtg                                                                        >  1:778843/1‑150 (MQ=255)
                                                                                 aaaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGc                                                                      <  2:1205207/150‑1 (MQ=255)
                                                                                  aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                     >  1:698540/1‑150 (MQ=255)
                                                                                  aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                     <  1:786336/150‑1 (MQ=255)
                                                                                  aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                     <  1:824789/150‑1 (MQ=255)
                                                                                  aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCACCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                     <  2:1227171/150‑1 (MQ=255)
                                                                                  aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCACCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                     <  2:1232213/150‑1 (MQ=255)
                                                                                      ccTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATaa                                                                                                                                     >  2:1159056/1‑82 (MQ=255)
                                                                                      ccTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATaa                                                                                                                                     <  1:1159056/82‑1 (MQ=255)
                                                                                      ccTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGgaag                                                                 >  1:1124214/1‑150 (MQ=255)
                                                                                        tttAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGgaagaa                                                               >  1:760997/1‑150 (MQ=255)
                                                                                        tttAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGgaagaa                                                               <  1:320978/150‑1 (MQ=255)
                                                                                          tAATCAGGAGCAGCTGAACGGGCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaa                                                             <  1:1261112/150‑1 (MQ=255)
                                                                                            aTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAAc                                                           <  1:1075148/150‑1 (MQ=255)
                                                                                            aTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAAc                                                           <  1:1100070/150‑1 (MQ=255)
                                                                                            aTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAAc                                                           >  1:1227658/1‑150 (MQ=255)
                                                                                            aTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAAc                                                           <  1:824394/150‑1 (MQ=255)
                                                                                                ggAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTc                                                       <  1:265838/150‑1 (MQ=255)
                                                                                                ggAGCACCTGAACGGTCGTGCGCAGGTCCTGGATCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCAAAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTc                                                       >  2:581259/1‑150 (MQ=255)
                                                                                                 gAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGGTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCg                                                      <  2:1399412/150‑1 (MQ=255)
                                                                                                    gagcTGAAGGGTCGTGCGCAGGTCGTGTTTCTCTACAGTATTACTCTCGGCGAGCTTTTCCCGCTTAAAGGATATGTGGTTGGGGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGTCGGAAGAAAAACAGTCGATg                                                   <  2:1087500/149‑1 (MQ=255)
                                                                                                         gAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGTTGTGGCGGAAGAAAAACAGTCGATGTTAgg                                              >  1:604709/1‑150 (MQ=255)
                                                                                                         gAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTgg                                              >  1:720063/1‑150 (MQ=255)
                                                                                                            cGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAg                                           <  2:3763/150‑1 (MQ=255)
                                                                                                                    cgcAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAg                                   >  2:272750/1‑150 (MQ=255)
                                                                                                                       aGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGtgc                                <  2:1106654/150‑1 (MQ=255)
                                                                                                                         gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCGATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGtgctg                              <  2:1418564/150‑1 (MQ=255)
                                                                                                                         gTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGtgctg                              <  2:243323/150‑1 (MQ=255)
                                                                                                                                 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:645215/1‑150 (MQ=255)
                                                                                                                                                     gccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTATCTGGCTACCTTACAGCACGATGt  >  2:650244/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACAATGTTGATGTTGCTGCCAGTGCCAATGGCGTGAGCGGCGATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGT  >  NC_000913/921426‑921724

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: ‑