Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 137,018 G→C 6.7% intergenic (‑101/‑65) 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,0180GC 6.7% 325.7 / 11.6 105intergenic (‑101/‑65)yacC/cueOPulS_OutS family protein/multicopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (50/47);  new base C (5/2);  total (56/49)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 105 total aligned reads displayed.

ACGGCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTG  >  NC_000913/136871‑137164
                                                                                                                                                   |                                                                                                                                                  
aCGGCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTa                                                                                                                                                  <  2:107759/150‑1 (MQ=255)
    ccATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                              <  2:796210/150‑1 (MQ=255)
    ccATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                              >  1:1164868/1‑150 (MQ=255)
    ccATCAGGCTGCCGAATAACACTGTGCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                              <  1:1592110/150‑1 (MQ=255)
     cATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAgggg                                                                                                                                             <  1:1184843/150‑1 (MQ=255)
         aGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAgg                                                                                                                                         <  2:761093/150‑1 (MQ=255)
          ggCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                        <  1:88451/150‑1 (MQ=255)
           gCTGTCGAATAACACTGTTCTGAAAAACTTCTTAATAGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggtc                                                                                                                                       <  2:134738/150‑3 (MQ=255)
                cGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGc                                                                                                                                  >  1:883344/1‑150 (MQ=255)
                cGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGc                                                                                                                                  <  1:499105/150‑1 (MQ=255)
                cGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGc                                                                                                                                  <  2:1249777/150‑1 (MQ=255)
                cGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGc                                                                                                                                  >  2:1264548/1‑150 (MQ=255)
                 gAATAACACTGTTTTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCt                                                                                                                                 >  2:1158276/1‑150 (MQ=255)
                   aTAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCa                                                                                                                               <  2:260211/150‑1 (MQ=255)
                    tAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCaa                                                                                                                              >  1:1404087/1‑150 (MQ=255)
                     aaCACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAAc                                                                                                                             >  1:596476/1‑150 (MQ=255)
                      acacTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACg                                                                                                                            >  1:1164458/1‑150 (MQ=255)
                       cacTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGt                                                                                                                           >  2:1572352/1‑150 (MQ=255)
                         cTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCAATATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGttt                                                                                                                         >  2:724031/1‑150 (MQ=255)
                          tGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTg                                                                                                                        <  1:101103/150‑1 (MQ=255)
                             tCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtt                                                                                                                     <  2:222207/150‑1 (MQ=255)
                               tGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                   >  1:1429996/1‑150 (MQ=255)
                                gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTg                                                                                                                  >  2:1593172/1‑150 (MQ=255)
                                  aaaaCGTCTTCATTGCTTCTACCAGGGGTTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                <  2:448403/150‑1 (MQ=255)
                                  aaaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCTCACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                <  1:68458/150‑1 (MQ=255)
                                  aaaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                <  1:626179/150‑1 (MQ=255)
                                   aaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGttt                                                                                                               >  1:684383/1‑150 (MQ=255)
                                      cGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGc                                                                                                            >  1:1680923/1‑150 (MQ=255)
                                       gTCTTCATTGCTTCTACCAGGGGCTTCAAGAAGATAAACGTTAGCTCACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGAGTCTGCGGCTGGACCTTCCGGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGcc                                                                                                           <  1:1090364/150‑1 (MQ=255)
                                       gTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGcc                                                                                                           <  2:685355/150‑1 (MQ=255)
                                           tCATTGCTTCTACCAGGGGCTTAAAGACTATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCt                                                                                                       <  2:253261/150‑1 (MQ=255)
                                            cATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCtt                                                                                                      >  2:1462378/1‑150 (MQ=255)
                                            cATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCtt                                                                                                      >  2:732122/1‑150 (MQ=255)
                                            cATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCtt                                                                                                      <  2:1128541/150‑1 (MQ=255)
                                              ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                    >  1:1342344/1‑150 (MQ=255)
                                               tGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAg                                                                                                   <  1:1572352/150‑1 (MQ=255)
                                                gCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGa                                                                                                  >  2:553471/1‑150 (MQ=255)
                                                   tCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCGCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATa                                                                                               <  2:1052719/150‑1 (MQ=255)
                                                       ccAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGa                                                                                           <  1:125926/150‑1 (MQ=255)
                                                       ccAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGa                                                                                           <  1:47201/150‑1 (MQ=255)
                                                         aGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGATCGAGGCGACGACCATTATTTTCGATTAAATTTTGTAAGCGGCTTGACATTCCCGTAAGGGCAAGGACTATGCTCAAAGGTTGATTTTGTTTCGCGTGCATAACAATAAGGaaa                                                                                         >  1:115855/1‑150 (MQ=255)
                                                          ggggCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAAt                                                                                        <  1:732122/150‑1 (MQ=255)
                                                             gCTTAAATAAGATAAAAGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAAc                                                                                     <  1:78926/150‑1 (MQ=255)
                                                             gCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAAc                                                                                     >  1:1027084/1‑150 (MQ=255)
                                                                 aaaGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAATTAACTAtt                                                                                 >  1:1249866/1‑149 (MQ=255)
                                                                 aaaGAAGATAAACGTTAGAACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCTTAAGGGGAAGGACTATTCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATg                                                                                 >  1:1256530/1‑150 (MQ=255)
                                                                  aagaagATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGc                                                                                >  2:221062/1‑150 (MQ=255)
                                                                   agaagaTAACCGTTAGCACACTGCGTGAGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCa                                                                               <  1:341275/150‑1 (MQ=255)
                                                                      agaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcg                                                                            >  1:364381/1‑150 (MQ=255)
                                                                         tAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcg                                                                         >  2:1303339/1‑150 (MQ=255)
                                                                         tAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcg                                                                         <  1:1593172/150‑1 (MQ=255)
                                                                             cGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAt                                                                     >  1:571676/1‑150 (MQ=255)
                                                                                tAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                                  <  2:883344/150‑1 (MQ=255)
                                                                                  gCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCtt                                                                >  2:1161565/1‑150 (MQ=255)
                                                                                  gCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCtt                                                                >  2:526867/1‑150 (MQ=255)
                                                                                  gCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCtt                                                                >  2:1073079/1‑150 (MQ=255)
                                                                                   cacaCTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTa                                                               <  2:684383/150‑1 (MQ=255)
                                                                                   cacaCTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTa                                                               >  2:513574/1‑150 (MQ=255)
                                                                                   cacaCTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTa                                                               <  1:1131999/150‑1 (MQ=255)
                                                                                    acacTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaa                                                              <  2:1027084/150‑1 (MQ=255)
                                                                                     cacTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                                             <  1:283038/150‑1 (MQ=255)
                                                                                       cTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAat                                                           <  2:1342344/150‑1 (MQ=255)
                                                                                       cTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAat                                                           <  1:1303339/150‑1 (MQ=255)
                                                                                       cTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAat                                                           >  2:1157124/1‑150 (MQ=255)
                                                                                         gcgagcgaGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAatat                                                         >  1:931081/1‑150 (MQ=255)
                                                                                           gagcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                       <  2:1164458/150‑1 (MQ=255)
                                                                                             gcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCg                                                     <  2:1164868/150‑1 (MQ=255)
                                                                                               gagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                                   <  1:1642017/150‑1 (MQ=255)
                                                                                                    gacgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTCAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCt                                              <  1:1468392/150‑1 (MQ=255)
                                                                                                    gacgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCt                                              >  1:1613553/1‑150 (MQ=255)
                                                                                                          cATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                                        >  2:429099/1‑150 (MQ=255)
                                                                                                            ttattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTgg                                      <  1:1264548/150‑1 (MQ=255)
                                                                                                             tattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGc                                     <  2:1680923/150‑1 (MQ=255)
                                                                                                              attGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCt                                    <  1:744441/150‑1 (MQ=255)
                                                                                                                 gTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCg                                 >  2:537458/1‑150 (MQ=255)
                                                                                                                   cGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                                             <  2:1577514/120‑1 (MQ=255)
                                                                                                                   cGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                                             >  1:1577514/1‑120 (MQ=255)
                                                                                                                    gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGATCAACGTTTTATTTTTTTTCGCCAGCTTAAGAAAAAGGAAATAAATATGCAACGTCGTGATTTCTTAAAATATTAAGTAGCGCTGGGTGTGGCTTCGGCt                              >  1:1084115/1‑150 (MQ=255)
                                                                                                                      ttAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTTATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt                            <  1:429099/150‑1 (MQ=255)
                                                                                                                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc                          >  1:1485759/1‑150 (MQ=255)
                                                                                                                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTAGTGATTTCTTAAAATATTACGTCGCGCTGGGTGTGGCTTCGGCTTTGc                          >  1:1715543/1‑150 (MQ=255)
                                                                                                                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTAAGCTCAACGTTTGATTTTGTTTAGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGTCtttta                          >  2:850770/1‑148 (MQ=255)
                                                                                                                          aTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCg                        <  1:1462378/150‑1 (MQ=255)
                                                                                                                             gtgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGACGcgg                     >  1:1072414/1‑148 (MQ=255)
                                                                                                                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                  <  1:1018883/150‑1 (MQ=255)
                                                                                                                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTgg                  >  1:612357/1‑150 (MQ=255)
                                                                                                                                     ggCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcg             <  1:553471/150‑1 (MQ=255)
                                                                                                                                     gcatagtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcg             >  1:900319/8‑150 (MQ=255)
                                                                                                                                      gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTATGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcgc            >  1:529010/1‑150 (MQ=255)
                                                                                                                                       cTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCcgcg           <  1:448635/150‑1 (MQ=255)
                                                                                                                                          gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg        >  1:528139/3‑150 (MQ=255)
                                                                                                                                          gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg        >  2:759144/3‑150 (MQ=255)
                                                                                                                                          gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg        >  2:798890/3‑150 (MQ=255)
                                                                                                                                          gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg        >  1:708693/3‑150 (MQ=255)
                                                                                                                                          gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg        >  1:1409383/1‑150 (MQ=255)
                                                                                                                                            ccTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATTAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa      <  1:1323489/150‑1 (MQ=255)
                                                                                                                                            ccTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa      >  1:243943/1‑150 (MQ=255)
                                                                                                                                            ccTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa      >  1:33272/1‑150 (MQ=255)
                                                                                                                                            ccTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa      >  1:426034/1‑150 (MQ=255)
                                                                                                                                                cccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTg  >  1:1626120/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGGCCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTG  >  NC_000913/136871‑137164

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