Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 137,018 G→C 5.3% 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 5.3% 264.0 / 5.2 76intergenic (‑101/‑65)yacC/cueOPulS_OutS family protein/multicopper oxidase (laccase)
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base C (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCG  >  NC_000913/136876‑137166
                                                                                                                                              |                                                                                                                                                    
cATCAGGCTGCCGAATAATACTGTTCTTAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                                >  4:671792/1‑149 (MQ=255)
cATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAggg                                                                                                                                                >  3:200024/1‑149 (MQ=255)
   cAGGCTGCAGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGaa                                                                                                                                             <  2:566952/149‑1 (MQ=255)
      gCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGa                                                                                                                                          <  3:236878/149‑1 (MQ=255)
       cTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCCTTACGGGAAggtc                                                                                                                                         <  3:93150/149‑3 (MQ=255)
        tGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACt                                                                                                                                        <  3:618942/149‑1 (MQ=255)
            gAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                            >  1:687865/1‑141 (MQ=255)
            gAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAg                                                                                                                                            <  2:687865/141‑1 (MQ=255)
             aaTAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCt                                                                                                                                   >  2:290678/1‑149 (MQ=255)
                   acTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGt                                                                                                                             >  3:232909/1‑149 (MQ=255)
                       ttCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGa                                                                                                                         <  3:668283/149‑1 (MQ=255)
                        tCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAt                                                                                                                        <  1:468855/149‑1 (MQ=255)
                          tGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAttt                                                                                                                      <  4:412319/149‑1 (MQ=255)
                           gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGTGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                     >  3:569449/1‑149 (MQ=255)
                           gAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                     >  2:207451/1‑149 (MQ=255)
                                aCGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTc                                                                                                                <  2:77626/149‑1 (MQ=255)
                                   tCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGcc                                                                                                             >  1:532332/1‑149 (MQ=255)
                                   tCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGcc                                                                                                             >  2:79706/1‑149 (MQ=255)
                                        aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTTTCTGCGGCTTGACCTTCCCGTAAGGGGTAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                  <  3:45487/139‑1 (MQ=255)
                                        aTTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGtt                                                                                                                  >  4:45487/1‑139 (MQ=255)
                                         ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCaa                                                                                                                                <  3:47055/124‑1 (MQ=255)
                                         ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCaa                                                                                                                                >  4:47055/1‑124 (MQ=255)
                                         ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                     >  4:422627/1‑135 (MQ=255)
                                         ttGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGAtttt                                                                                                                     <  3:422627/135‑1 (MQ=255)
                                          tGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTaa                                                                                                      <  1:542001/149‑1 (MQ=255)
                                           gCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAg                                                                                                     <  1:78327/149‑1 (MQ=255)
                                             ttCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGTCCTTCCCCTTACGGGAAggtcaa                                                                                                                                       >  3:193232/1‑109 (MQ=255)
                                             ttCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGTCCTTCCCCTTACGGGAAggtcaa                                                                                                                                       <  4:193232/113‑5 (MQ=255)
                                                tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                <  1:207451/149‑1 (MQ=255)
                                                tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                <  1:290678/149‑1 (MQ=255)
                                                tACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATaa                                                                                                <  4:121037/149‑1 (MQ=255)
                                                  ccAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAgg                                                                                              >  3:633326/1‑149 (MQ=255)
                                                    aGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                               >  3:528708/1‑146 (MQ=255)
                                                    aGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAg                                                                                               <  4:528708/146‑1 (MQ=255)
                                                     ggggCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGaaa                                                                                           >  1:22717/1‑149 (MQ=255)
                                                          ttAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGGAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACt                                                                                      >  4:653307/1‑149 (MQ=255)
                                                           tAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                                     <  1:567620/149‑1 (MQ=255)
                                                           tAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTa                                                                                     <  2:532332/149‑1 (MQ=255)
                                                              agaagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGc                                                                                  >  2:465716/1‑149 (MQ=255)
                                                                aagaTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCaa                                                                                >  4:433567/1‑149 (MQ=255)
                                                                   aTAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                                             >  3:390188/1‑149 (MQ=255)
                                                                    tAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtc                                                                            >  4:316540/1‑149 (MQ=255)
                                                                      aaCGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                                          <  4:232909/149‑1 (MQ=255)
                                                                         gTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCt                                                                                                         >  2:2678/1‑115 (MQ=255)
                                                                         gTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCt                                                                                                         <  1:2678/115‑1 (MQ=255)
                                                                         gTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAt                                                                       >  2:262692/1‑149 (MQ=255)
                                                                          ttAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAtt                                                                      <  1:465716/149‑1 (MQ=255)
                                                                            aGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTTCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                                    <  1:95000/149‑1 (MQ=255)
                                                                                cacTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTc                                                                                                                >  3:666769/1‑101 (MQ=255)
                                                                                cacTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTc                                                                                                                <  4:666769/101‑1 (MQ=255)
                                                                                   tGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAGGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAat                                                             >  1:44155/1‑149 (MQ=255)
                                                                                     cgagcgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAatat                                                           <  2:22717/149‑1 (MQ=255)
                                                                                         cgagGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCg                                                       <  1:262692/149‑1 (MQ=255)
                                                                                                  gacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              >  3:77076/1‑149 (MQ=255)
                                                                                                  gacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              <  4:77076/149‑1 (MQ=255)
                                                                                                      attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                                          >  1:506124/1‑149 (MQ=255)
                                                                                                       ttattGTCGATTAAATTGTGTCCGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGtgtg                                         <  4:569449/149‑1 (MQ=255)
                                                                                                          ttGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCt                                      >  1:41067/1‑149 (MQ=255)
                                                                                                            gTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTc                                    <  2:385843/149‑1 (MQ=255)
                                                                                                             tCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCg                                   >  1:158610/1‑149 (MQ=255)
                                                                                                               gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                         >  3:487478/1‑125 (MQ=255)
                                                                                                               gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                                         <  4:487478/125‑1 (MQ=255)
                                                                                                                     aTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc                           <  2:158610/149‑1 (MQ=255)
                                                                                                                     aTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc                           >  3:650702/1‑149 (MQ=255)
                                                                                                                      ttGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCg                          >  1:667761/1‑149 (MQ=255)
                                                                                                                         tgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgt                                            >  3:178814/1‑128 (MQ=255)
                                                                                                                         tgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgt                                            <  4:178814/128‑1 (MQ=255)
                                                                                                                          gtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt                      <  1:79706/149‑1 (MQ=255)
                                                                                                                              gCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAg                  >  1:222583/1‑149 (MQ=255)
                                                                                                                                ggCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGcc                >  2:195669/1‑149 (MQ=255)
                                                                                                                                 gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              <  2:39567/118‑1 (MQ=255)
                                                                                                                                 gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                                              >  1:39567/1‑118 (MQ=255)
                                                                                                                                     gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCa           <  2:222583/149‑1 (MQ=255)
                                                                                                                                      tccTTCCCCTTACGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg          >  4:19143/2‑149 (MQ=255)
                                                                                                                                      aCCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAg          >  1:499519/1‑149 (MQ=255)
                                                                                                                                        cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTa        >  3:25318/1‑149 (MQ=255)
                                                                                                                                              gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCg  >  1:132092/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CATCAGGCTGCCGAATAACACTGTTCTGAAAAACGTCTTCATTGCTTCTACCAGGGGCTTAAAGAAGATAAACGTTAGCACACTGCGAGCGAGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGTGGAGCCGCGCAGTATTTGCG  >  NC_000913/136876‑137166

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