Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 921,617 T→C 5.1% I424T (ATT→ACT)  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,6170TC 5.1% 242.2 / 2.5 78I424T (ATT→ACT) macBABC‑type tripartite efflux pump ATP binding/membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

ATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTG  >  NC_000913/921470‑921766
                                                                                                                                                   |                                                                                                                                                     
attattTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTc                                                                                                                                                     <  2:1367613/150‑1 (MQ=255)
attattTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTc                                                                                                                                                     <  1:928912/150‑1 (MQ=255)
attattTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGACGTGGTTGGCGAGGTGATTc                                                                                                                                                     <  2:1444774/150‑1 (MQ=255)
   attTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTgg                                                                                                                                                  <  1:828726/150‑1 (MQ=255)
    ttttAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGt                                                                                                                                                 >  2:1102270/1‑150 (MQ=255)
    ttttAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGt                                                                                                                                                 <  2:1166064/150‑1 (MQ=255)
     tatAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGATCGGTCGTGCGCAGGTCGTGGTTCACGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTc                                                                                                                                                >  1:886365/3‑150 (MQ=255)
         aTGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCa                                                                                                                                            >  2:541067/1‑150 (MQ=255)
         aTGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCa                                                                                                                                            >  1:954315/1‑150 (MQ=255)
         aTGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCa                                                                                                                                            <  1:811589/150‑1 (MQ=255)
          tGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCaa                                                                                                                                           <  2:1437574/150‑1 (MQ=255)
          tGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCaa                                                                                                                                           <  2:1220317/150‑1 (MQ=255)
           gTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAat                                                                                                                                          <  1:1535889/150‑1 (MQ=255)
              tATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGACAGCTTTTCACCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATg                                                                                                                                       >  2:1015383/1‑150 (MQ=255)
               aTGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGc                                                                                                                                      >  1:833025/1‑150 (MQ=255)
                   cATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGc                                                                                                                                  >  1:1023487/1‑150 (MQ=255)
                     tGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTc                                                                                                                                                >  2:100543/1‑134 (MQ=255)
                     tGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTc                                                                                                                                                <  1:100543/134‑1 (MQ=255)
                        ccTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCagag                                                                                                                             >  1:953665/1‑150 (MQ=255)
                        ccTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCagag                                                                                                                             >  1:700556/1‑150 (MQ=255)
                          ttCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTc                                                                                                                           >  2:568622/1‑150 (MQ=255)
                           tCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGCTGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCa                                                                                                                          >  2:266005/1‑150 (MQ=255)
                             aGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGGCAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGACTCTGGCGGGCATATTGCCGACCAGAATCacc                                                                                                                        >  2:131321/1‑148 (MQ=255)
                             aGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCAtt                                                                                                                        <  2:1039695/150‑1 (MQ=255)
                              gTGAAGGAAACACCTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTgg                                                                                                                      >  1:1016171/1‑150 (MQ=255)
                               tGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGGCAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGACTCTGGCGGGCATATTGCCGACCAGAATCacctc                                                                                                                      <  1:131321/150‑5 (MQ=37)
                                 aaGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAATCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGgtg                                                                                                                    <  1:187890/150‑1 (MQ=255)
                                  aGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCACCCATAAAGCAGATGTTGTTGGCGAGTTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGgtgt                                                                                                                   <  2:343572/150‑1 (MQ=255)
                                   ggAAACACCTTTAATCAGGAGCAGCTGAACAGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGtgtg                                                                                                                  <  2:273079/150‑1 (MQ=255)
                                     aaaCACCTTTAATCAGGACCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGc                                                                                                                <  2:1508076/150‑1 (MQ=255)
                                      aaCACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCg                                                                                                               >  1:924082/1‑150 (MQ=255)
                                        cacCTTTAATCAGTAGCATCTTAACGGTGGTGCGCAGGTCGTGGTTGTCGACAGTAATAGTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCTAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCTga                                                                                                             <  2:408936/150‑1 (MQ=255)
                                        cacCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGga                                                                                                             >  1:906113/1‑150 (MQ=255)
                                        cacCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGACTCTGGTCGGCAATATGCCCGCCagag                                                                                                                             >  2:48123/1‑134 (MQ=255)
                                        cacCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGACTCTGGTCGGCAATATGCCCGCCagag                                                                                                                             <  1:48123/134‑1 (MQ=255)
                                              tAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaa                                                                                                       <  2:65362/150‑1 (MQ=255)
                                               aaTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGaaaaa                                                                                                      >  2:74087/1‑150 (MQ=255)
                                                aTCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAAc                                                                                                     <  1:120997/150‑1 (MQ=255)
                                                 tCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGGAGAAAAACa                                                                                                    <  2:1038902/150‑1 (MQ=255)
                                                  cAGGAGCAGCTGAACGGTCGTGAGCAGGTCGTGGTTCTCCACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAg                                                                                                   >  2:1144874/1‑150 (MQ=255)
                                                   aGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGt                                                                                                  >  1:1525814/1‑150 (MQ=255)
                                                    ggAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTc                                                                                                 <  2:77617/150‑1 (MQ=255)
                                                     gAGCAGCTGAACGTTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGACGCCAGCTTTTCCGCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGTTGTGGCGTAAGAAAAACAGTCg                                                                                                <  2:1160674/150‑1 (MQ=255)
                                                      agcagcTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGa                                                                                               <  1:258870/150‑1 (MQ=255)
                                                         agcTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAACTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGt                                                                                            >  2:403160/1‑150 (MQ=255)
                                                            tGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTg                                                                                         <  2:1525814/150‑1 (MQ=255)
                                                            tGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTg                                                                                         >  2:1176123/1‑150 (MQ=255)
                                                              aaCGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGt                                                                                       <  1:561664/150‑1 (MQ=255)
                                                                       gcgcAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTaaa                                                                              <  1:929631/150‑1 (MQ=255)
                                                                         gcAGGTCGTGGTTCTCGACAGTCATACTCGCCGCCAGCTTTTCGCCCATAAAGCAGCTGTGGTTGGCGCGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGAGGTTTGGTAGCAGTAAAgg                                                                            >  2:1267307/1‑149 (MQ=255)
                                                                                       tcGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGc                                                              <  2:1023487/150‑1 (MQ=255)
                                                                                           cAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTAcc                                                          >  2:1056104/1‑150 (MQ=255)
                                                                                            aGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCt                                                         >  1:1166518/1‑150 (MQ=255)
                                                                                            aGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCt                                                         >  1:1172037/1‑150 (MQ=255)
                                                                                            aGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCt                                                         >  2:1432099/1‑150 (MQ=255)
                                                                                             gTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGGGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCtt                                                        <  1:92932/150‑1 (MQ=255)
                                                                                                  aCTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGGGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGc                                                   <  2:722160/150‑1 (MQ=255)
                                                                                                    tCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCAc                                                 <  1:1169879/150‑1 (MQ=255)
                                                                                                          ccAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTc                                           >  2:1259052/1‑150 (MQ=255)
                                                                                                                    ccccATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGtt                                 >  2:1266121/1‑150 (MQ=255)
                                                                                                                     cccATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTa                                >  2:1134384/1‑150 (MQ=255)
                                                                                                                        aTAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATgg                             >  1:1056224/1‑150 (MQ=255)
                                                                                                                         tAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATggg                            <  1:1056104/150‑1 (MQ=255)
                                                                                                                                   gtgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTgg                  <  1:1036957/150‑1 (MQ=255)
                                                                                                                                    tgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGc                 <  2:1172037/150‑1 (MQ=255)
                                                                                                                                    tgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGc                 <  2:1166518/150‑1 (MQ=255)
                                                                                                                                    tgGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGGCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGc                 <  1:1400037/150‑1 (MQ=255)
                                                                                                                                        tGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTaa             >  1:465381/1‑150 (MQ=255)
                                                                                                                                        tGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTaa             >  1:851216/1‑150 (MQ=255)
                                                                                                                                         ggCGATGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAAc            <  1:1015383/150‑1 (MQ=255)
                                                                                                                                         ggCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAAc            <  2:924082/150‑1 (MQ=255)
                                                                                                                                          gCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACt           <  1:568622/150‑1 (MQ=255)
                                                                                                                                                tGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGGTGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTa     >  1:1111979/1‑150 (MQ=255)
                                                                                                                                                 taTTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGGGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTAc    <  2:600785/149‑1 (MQ=255)
                                                                                                                                                 gATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTAc    >  1:864773/1‑150 (MQ=255)
                                                                                                                                                 gATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTAc    >  1:563471/1‑150 (MQ=255)
                                                                                                                                                  aTTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACt   >  2:679610/1‑150 (MQ=255)
                                                                                                                                                  aTTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACt   <  1:1259052/150‑1 (MQ=255)
                                                                                                                                                   ttCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTg  <  1:266005/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTATTTTAATGTCTATGGCATGACCTTCAGTGAAGGAAACACCTTTAATCAGGAGCAGCTGAACGGTCGTGCGCAGGTCGTGGTTCTCGACAGTAATACTCGCCGCCAGCTTTTCCCCCATAAAGCAGATGTGGTTGGCGAGGTGATTCTGGTCGGCAATATGCCCGCCAGAGTCATTGGTGTGGCGGAAGAAAAACAGTCGATGTTTGGTAGCAGTAAAGTGCTGCGTGTCTGGCTACCTTACAGCACGATGTCCGGGCGAGTTATGGGCCAGTCGTGGCTTAACTCCATTACTG  >  NC_000913/921470‑921766

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: