Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 137,022 G→C 7.0% intergenic (‑105/‑61) yacC ← / → cueO PulS_OutS family protein/multicopper oxidase (laccase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913137,0220GC 7.0% 614.0 / 29.7 171intergenic (‑105/‑61)yacC/cueOPulS_OutS family protein/multicopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (76/83);  new base C (6/6);  total (82/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 179 total aligned reads displayed.

GCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTA  >  NC_000913/136874‑137160
                                                                                                                                                    |                                                                                                                                          
gCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAg                                                                                                                                            <  1:1680295/149‑1 (MQ=255)
     cAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                       <  3:927741/149‑1 (MQ=255)
     cAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                       <  2:951031/149‑1 (MQ=255)
      aGGCTGCCGAATAGCACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                      <  3:63341/149‑1 (MQ=255)
      aGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                      <  2:1054408/149‑1 (MQ=255)
       ggCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                      >  1:1365071/1‑148 (MQ=255)
       ggCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                      <  2:1365071/148‑1 (MQ=255)
         cTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGTCCTTCCCCTTACGGGAAggtc                                                                                                                                   <  3:959143/149‑3 (MQ=255)
            ccGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAg                                                                                                                                      <  1:1292365/143‑1 (MQ=255)
            ccGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAg                                                                                                                                      >  2:1292365/1‑143 (MQ=255)
                      cTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                          <  3:662716/129‑1 (MQ=255)
                         ttCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGa                                                                                                                   <  1:998424/149‑1 (MQ=255)
                         ttCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGa                                                                                                                   <  2:328350/149‑1 (MQ=255)
                             gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGttt                                                                                                                     <  3:1069812/143‑1 (MQ=255)
                             gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGttt                                                                                                                     >  4:1069812/1‑143 (MQ=255)
                              aaaaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTg                                                                                                              <  3:1358007/149‑1 (MQ=255)
                                aaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                            >  2:273133/1‑149 (MQ=255)
                                       ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGtt                                                                                                                      <  4:1737217/132‑1 (MQ=255)
                                       ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTg                                                                                                     >  1:1415212/1‑149 (MQ=255)
                                       ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTg                                                                                                     >  1:1047247/1‑149 (MQ=255)
                                       ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTg                                                                                                     >  2:1027600/1‑149 (MQ=255)
                                          aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAgg                                                                                                                                           <  4:302907/108‑1 (MQ=255)
                                          aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAgg                                                                                                                                           >  3:302907/1‑108 (MQ=255)
                                          aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCtt                                                                                                  >  1:145413/1‑149 (MQ=255)
                                           ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAg                                                                                                                                            <  1:241199/106‑1 (MQ=255)
                                           ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTc                                                                                                          >  2:246939/1‑140 (MQ=255)
                                            tGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                <  4:1346834/149‑1 (MQ=255)
                                               ttCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGaa                                                                                             >  3:699317/1‑149 (MQ=255)
                                                  tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAt                                                                                                                  >  2:81519/1‑125 (MQ=255)
                                                  tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                          >  1:1633208/1‑149 (MQ=255)
                                                   aCCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                         <  2:1253028/149‑1 (MQ=255)
                                                     cAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGa                                                                                       >  4:1023276/1‑149 (MQ=255)
                                                     cAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGa                                                                                       >  2:1766307/1‑149 (MQ=255)
                                                       ggggCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGt                                                                                                             >  4:909331/1‑125 (MQ=255)
                                                       ggggCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGaaa                                                                                     <  3:432031/149‑1 (MQ=255)
                                                          gCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATaa                                                                                  <  4:1284755/149‑1 (MQ=255)
                                                          gCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATaa                                                                                  <  2:1633208/149‑1 (MQ=255)
                                                          gCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATaa                                                                                  <  1:746788/149‑1 (MQ=255)
                                                          gCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATaa                                                                                  <  1:1766307/149‑1 (MQ=255)
                                                             tAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                               <  4:699317/149‑1 (MQ=255)
                                                               aagaagATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                    >  2:1181831/1‑94 (MQ=255)
                                                               aagaagATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                    <  1:1181831/94‑1 (MQ=255)
                                                                 gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAAc                                                                                                                         <  2:1457034/103‑1 (MQ=255)
                                                                 gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAAc                                                                                                                         >  1:1457034/1‑103 (MQ=255)
                                                                 gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggtc                                                                                                                                   >  1:595334/1‑91 (MQ=255)
                                                                 gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggtc                                                                                                                                   >  1:431679/1‑91 (MQ=255)
                                                                 gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggtc                                                                                                                                   <  2:595334/93‑3 (MQ=255)
                                                                    gaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                               <  2:548544/110‑1 (MQ=255)
                                                                    gaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                               >  1:548544/1‑110 (MQ=255)
                                                                    gaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTAt                                                                              <  2:1434323/143‑1 (MQ=255)
                                                                     aTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                                       >  2:306579/1‑149 (MQ=255)
                                                                      tAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTAt                                                                              <  2:717711/141‑1 (MQ=255)
                                                                       aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                <  4:167872/122‑1 (MQ=255)
                                                                       aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                >  3:167872/1‑122 (MQ=255)
                                                                       aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                          >  1:1580898/1‑128 (MQ=255)
                                                                       aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                          <  2:1580898/128‑1 (MQ=255)
                                                                        aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                       <  1:407444/82‑1 (MQ=255)
                                                                        aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                    >  1:1810036/1‑149 (MQ=255)
                                                                        aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                    >  1:308117/1‑149 (MQ=255)
                                                                              aGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                              <  1:396964/149‑1 (MQ=255)
                                                                               gCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCt                                                             >  4:44616/1‑149 (MQ=255)
                                                                                    cTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                                        <  1:273133/149‑1 (MQ=255)
                                                                                          gcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                            <  4:662688/91‑1 (MQ=255)
                                                                                            gagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGt                                                >  2:416951/1‑149 (MQ=255)
                                                                                             agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                               >  3:1850623/1‑149 (MQ=255)
                                                                                             agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                               <  3:44616/149‑1 (MQ=255)
                                                                                                       cATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAttt                                                                                                                >  2:1369178/1‑74 (MQ=255)
                                                                                                       cATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAttt                                                                                                                <  1:1369178/74‑1 (MQ=255)
                                                                                                       cATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGa                                                                  <  3:658740/120‑1 (MQ=255)
                                                                                                       cATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGa                                                                  >  4:658740/1‑120 (MQ=255)
                                                                                                        attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                          <  2:832734/95‑1 (MQ=255)
                                                                                                        attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                          >  1:832734/1‑95 (MQ=255)
                                                                                                        attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                                    >  2:1650437/1‑149 (MQ=255)
                                                                                                          tattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTgg                                  >  2:386994/1‑149 (MQ=255)
                                                                                                            ttGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                >  2:402282/1‑85 (MQ=255)
                                                                                                              gTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTc                              <  2:1415212/149‑1 (MQ=255)
                                                                                                                 gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                   <  3:156295/125‑1 (MQ=255)
                                                                                                                 gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                   >  1:1022801/1‑125 (MQ=255)
                                                                                                                 gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                   >  4:156295/1‑125 (MQ=255)
                                                                                                                 gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                               >  1:1288545/1‑145 (MQ=255)
                                                                                                                    tAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt                        <  2:1830492/149‑1 (MQ=255)
                                                                                                                             tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg               >  4:852927/1‑149 (MQ=255)
                                                                                                                                gCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgt                                      <  4:1303174/123‑1 (MQ=255)
                                                                                                                                   gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                <  1:347026/142‑1 (MQ=255)
                                                                                                                                     ttGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcgcg       >  2:418691/1‑149 (MQ=255)
                                                                                                                                     ttGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcgcg       >  4:515965/1‑149 (MQ=255)
                                                                                                                                        tccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg    >  2:1228932/2‑149 (MQ=255)
                                                                                                                                        tccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTATTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGtgtg                                   >  3:1504218/2‑118 (MQ=255)
                                                                                                                                          cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                <  3:1319539/55‑1 (MQ=255)
                                                                                                                                          cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                >  4:1319539/1‑55 (MQ=255)
                                                                                                                                          cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa  >  4:1258933/1‑149 (MQ=255)
                                                                                                                                           ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGaaa                                                                                     >  1:1013024/1‑65 (MQ=255)
                                                                                                                                           ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGaaa                                                                                     <  2:1013024/65‑1 (MQ=255)
                                                                                                                                           ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                        >  3:887322/1‑110 (MQ=255)
                                                                                                                                            tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                               <  1:893134/118‑1 (MQ=255)
                                                                                                                                            tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                               >  2:893134/1‑118 (MQ=255)
                                                                                                                                              ccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCCTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCa     >  3:1290122/1‑142 (MQ=255)
                                                                                                                                                 tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                >  3:1266383/1‑128 (MQ=255)
                                                                                                                                                 tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg              >  4:1700134/1‑130 (MQ=255)
                                                                                                                                                 tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg              <  3:1700134/130‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTA  >  NC_000913/136874‑137160

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