Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 137,022 G→C 9.4% 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 9.4% 218.8 / 12.4 64intergenic (‑105/‑61)yacC/cueOPulS_OutS family protein/multicopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base C (2/4);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

CCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCG  >  NC_000913/136875‑137166
                                                                                                                                                   |                                                                                                                                                
ccATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAgg                                                                                                                                                 <  4:83212/149‑1 (MQ=255)
 cATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                                <  1:55410/149‑1 (MQ=255)
    cAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                             >  4:195180/1‑149 (MQ=255)
     aGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAg                                                                                                                                            <  3:357559/149‑1 (MQ=255)
     aGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAg                                                                                                                                            <  1:329592/149‑1 (MQ=255)
       gCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                          <  3:53179/149‑1 (MQ=255)
       gCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                          >  3:32149/1‑149 (MQ=255)
                 aaCACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCaa                                                                                                                                >  3:470196/1‑149 (MQ=255)
                    acTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGt                                                                                                                             >  3:135784/1‑149 (MQ=255)
                            gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                     <  4:271882/149‑1 (MQ=255)
                               aaaCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                  <  1:56379/149‑1 (MQ=255)
                                      ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                            <  4:202739/116‑1 (MQ=255)
                                      ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                            >  3:202739/1‑116 (MQ=255)
                                      ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggt                                                                                                                                          >  1:225043/1‑117 (MQ=255)
                                      ttCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggt                                                                                                                                          <  2:225043/118‑2 (MQ=255)
                                         aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAAc                                                                                                                               <  1:467888/126‑1 (MQ=255)
                                         aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAAc                                                                                                                               >  2:467888/1‑126 (MQ=255)
                                          ttGCTTCTACCAGGGGCTTCAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGTGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTa                                                                                                       >  3:281335/1‑149 (MQ=255)
                                          ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                             >  1:147243/1‑111 (MQ=255)
                                          ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                             <  2:147243/111‑1 (MQ=255)
                                           tGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                      <  4:135784/149‑1 (MQ=255)
                                              ttCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGaa                                                                                                   >  2:17656/1‑149 (MQ=255)
                                                 tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                <  3:441134/149‑1 (MQ=255)
                                                  aCCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                               <  3:52854/149‑1 (MQ=255)
                                                   ccAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAgg                                                                                              >  1:441978/1‑149 (MQ=255)
                                                    cAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                <  3:411789/146‑1 (MQ=255)
                                                    cAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                >  4:411789/1‑146 (MQ=255)
                                                           ttAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACt                                                                                      >  2:144795/1‑149 (MQ=255)
                                                           ttAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACt                                                                                      >  4:222817/1‑149 (MQ=255)
                                                            tAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                                     <  1:17656/149‑1 (MQ=255)
                                                              aagaagATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATg                                                                                   >  2:240029/1‑149 (MQ=255)
                                                              aagaagATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATg                                                                                   <  1:226243/149‑1 (MQ=255)
                                                               agaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGc                                                                                  >  3:118777/1‑149 (MQ=255)
                                                               agaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGAATGCGGCTTGAGCTACCCGTAAGGGGAAGGACTATGATCTTCGTTTGATTTTGTTTCGCCTTCTTAAGAATAAGGAAATAACTATGc                                                                                  >  2:319920/1‑149 (MQ=255)
                                                                gaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCa                                                                                 >  4:350370/1‑149 (MQ=255)
                                                                    aTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                                     <  1:332994/141‑1 (MQ=255)
                                                                    aTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                                     >  2:332994/1‑141 (MQ=255)
                                                                      aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                               >  1:27950/1‑129 (MQ=255)
                                                                      aaaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                               <  2:27950/129‑1 (MQ=255)
                                                                       aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                          <  1:319920/149‑1 (MQ=255)
                                                                       aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                          >  1:317424/1‑149 (MQ=255)
                                                                       aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                          <  3:285420/149‑1 (MQ=255)
                                                                                  acTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                                               <  4:118777/149‑1 (MQ=255)
                                                                                   cTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                                              <  1:393393/149‑1 (MQ=255)
                                                                                        agcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                         <  4:470196/149‑1 (MQ=255)
                                                                                         gcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTcc                                                        >  4:61034/1‑149 (MQ=255)
                                                                                                       attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                                          <  4:311039/149‑1 (MQ=255)
                                                                                                       attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                                          <  4:32149/149‑1 (MQ=255)
                                                                                                        ttattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              <  1:425039/144‑1 (MQ=255)
                                                                                                        ttattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              >  2:425039/1‑144 (MQ=255)
                                                                                                                 aTTAAATTTTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                      >  3:437942/1‑79 (MQ=255)
                                                                                                                 aTTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                      <  4:437942/79‑1 (MQ=255)
                                                                                                                  ttAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt                               >  4:365797/1‑149 (MQ=255)
                                                                                                                   tAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt                              <  2:162412/149‑1 (MQ=255)
                                                                                                                       ttGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                  >  1:184266/1‑61 (MQ=255)
                                                                                                                       ttGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                  <  2:184266/61‑1 (MQ=255)
                                                                                                                       ttGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCg                          >  2:61879/1‑149 (MQ=255)
                                                                                                                           gtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                      <  1:144795/149‑1 (MQ=255)
                                                                                                                            tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgtg                     <  3:61034/149‑1 (MQ=255)
                                                                                                                                 ggCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGcc                <  1:240029/149‑1 (MQ=255)
                                                                                                                                      gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTAt                                                                                    <  4:177803/74‑1 (MQ=38)
                                                                                                                                      gtccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTAt                                                                                    >  3:177803/3‑76 (MQ=38)
                                                                                                                                       aCCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              <  4:177700/113‑1 (MQ=255)
                                                                                                                                       aCCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              >  3:177700/1‑113 (MQ=255)
                                                                                                                                         cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa        >  2:34871/1‑149 (MQ=255)
                                                                                                                                               gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCg  >  4:170233/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCG  >  NC_000913/136875‑137166

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